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Abstract 


In  recent  years  such  progress  has  been  reported  in  the  fabrication  of 
•ultilayer  reflectors  for  x-ray  and  XUV  radiation  (1A  -  600A) . 

The  diaracteristic  optical  properties  of  Materials  at  these 
wavelengths  (dielectric  constants  are  aonplex  and  approach  unity),  allow 
one  to  obtain  solutions  to  Maxwell's  equations  for  a  quasi-periodic 
reflector  whose  layers  May  contain  arbitrary  index  gradients. 

This  solution  can  be  formulated  as  a  difference  equation  that 
propagates  the  Miplitude  reflectance  across  cad)  layer  pair.  The 
difference  equation  resembles  the  Airy  svnmaticn  for  single  layers,  but 
has  a  simpler  Ricatti  form. 

From  the  difference  equation  one  can  derive  design  criteria  for 
maximization  of  nultilayer  reflectivity.  These  criteria  provide  guidance 
in  the  selection  of  appropriate  multilayer  materials,  and  have  been  used 
to  derive  approximate  scaling  laws  for  multilayer  reflecting  properties. 

The  difference  aquation  forms  the  basis  for  a  non-perturbative 
analysis  of  wltilayer  reflectivity  in  the  presence  of  randan  thickness 
errors  (including  the  residual  loss  in  raflectivity  that  remains  when 
reflectance  monitoring  ia  used  to  ocnpensate  for  thickness  errors  during 
multilayer  fabrication) .  Under  certain  circunstanoes,  the  difference 
equation  can  be  used  to  snalyze  the  effect  of  interfacial  roughness  on 
multilayer  reflectivity.  Accurate  closed-form  solutions  to  such 
stodhastic  problems  can  be  found  by  neglecting  higher  order  powers  in  the 
incoherent  ccwponent  of  multilayer  reflectivity. 
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She  reflecting  properties  of  x-ray  multilayers  may  contain 
qualitative  signatures  that  correspond  to  these  different  kinds  of 
structural  defects. 

Our  theoretical  results  indicate  that  significant  constraints  on 
efficiency  and  field  of  view  are  involved  in  using  ulti layers  to  extend 
optical  technology  to  the  x-ray  regia*,  bangles  discussed  include 
•icrosoopes  operating  at  short  x-ray  wavelengths  (1&  -  2A) ,  share 
multilayers  can  provide  a  useful  level  of  spectral  selectivity,  and 
resonant  cavities  for  projected  x-ray  lasers  (50&  -  200&) ,  share  the 
coupling  of  the  intracavity  field  to  the  amplifier  can  be  strongly 
increased  if  the  cavity  configuration  is  baaed  on  two  multilayers 
that  are  tuned  to  reflect  at  normal  incidence. 
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Chapter  I  Introduction 


Section  1-1  -  Literature  Review 

Part  A)  X-Ray  Multilayer  Reflectors 

In  recent  years  such  progress  has  been  reported  in  the  fabrication  of 
reflecting  multilayers  having  period  lengths  in  the  100A  regime  (Haelbich 
and  Kunz,  1976;  Haelbich,  Segmuller,  and  Spiller,  1979;  Barbee  and 
Keith,  1979;  Spiller  et  al.,  1980;  Barbee,  1981;  Spiller,  1982b;  Underwood 
and  Barbee,  1982).  Such  devices  may  be  regarded  either  as 
short-wavelength  multilayer  coatings,  or  as  artificial  crystals  having 
large  spacings  between  the  diffracting  "planes".  (There  are  few  true 
crystals  having  spacings  larger  than  about  12A.) 

These  devices  should  prove  useful  for  the  spectral  region  below 
—'600A  (Spiller,  1972,  1976).  In  this  part  of  the  spectrin  all  materials 
have  complex  dielectric  constants  that  depart  only  slightly  from  unity; 
i.e.  materials  are  absorbing  and  are  incapable  of  providing  a  high 
single-surface  reflectivity. 

In  this  work  we  present  a  theoretical  analysis  of  the  reflecting 
properties  of  x-ray  multilayers.  In  the  present  chapter  we  review  the 
field,  and  also  discuss  related  topics  (and  associated  terminology')  frcm 
other  fields,  such  as  x-ray  diffraction  theory. 
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According  to  Barbee  (1982),  the  first  attenpts  to  fabricate  jr-ray 
multilayers  (generally  unsuccessful)  were  made  by  Koeppe  (1929)  and 
Deubner  (1930). 

More  recently,  Dinklage  (1967)  succeeded  in  fabricating 
iron/magnesiuti  reflectors  having  a  reasonably  stable  reflecting  power 
(roan  tenperature  half-life  about  one  year).  Dinklage  has  not  reported  an 
absolute  reflectivity  for  these  structures. 

Spiller  (1972,  1976)  recognized  that  in  an  absorbing  multilayer  it  is 
not  necessarily  desirable  to  minimize  the  absorption  of  both  the  high  and 
the  low  index  layers.  By  sulking  the  thickness  of  a  strongly  absorbing 
high  index  layer  substantially  thinner  than  one  quarter-wave,  one 
decreases  the  absorption  of  the  layer  sore  rapidly  than  one  does  the 
reflectivity,  allowing  in  principle  for  a  more  than  compensating  increase 
in  the  total  number  of  reflecting  layers.  larger  total  reflectivities  are 
then  obtained  when  a  dense,  high  2  material  is  used  in  high  index  layers 
that  are  alternated  with  low  absorption  (pacing  layers. 

Haelbich  and  Kunz  (1976)  achieved  1.21  peak  reflectivity  at  15QA  (IS* 
incidence  to  the  normal)  from  a  9  layer  structure  of  gold  and  carbon  that 
was  designed  according  to  the  above  considerations.  Carbon  was  chosen  as 
a  spacer  material  because  of  its  low  absorption,  and  because  it  does  not 
diffuse  readily  with  metals. 

Haelbich,  Segmuller,  and  Spiller  (1979)  used  a  rheniviv'tungsten  alloy 
(ReW)  as  a  high  index  layer  and  carbon  as  a  low  index  layer  to  obtain  9.5* 
peak  reflectivity  (15*  angle  of  incidence)  with  7  layers  st  200A.  The  ReK 
alloy  was  found  to  provide  ■soother  layers  than  did  other  high  density 
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Material*  such  as  gold. 

In  these  early  efforts,  the  number  of  participating  layers  that  the 
structures  could  contain  ms  limited  by  an  inability  to  precisely  control 
the  layer  thicknesses,  and  so  avoid  the  aocvnulaticn  of  dephasing  due  to 
random  errors  (Baelbieh  and  Kunz,  1976}  Baelbich,  Segnuller,  and 
Spiller,  1979).  In  order  to  reach  absorption-limited  reflectivities, 
x-ray  multilayers  must  contain  a  larger  mmfeer  of  layer  pairs  than  are 
required  in  optical  multilayers,  because  at  x-ray  wavelengths  the 
interaction  of  an  individual  unit  cell  with  the  incident  bean  is  weak. 

(As  a  rule  of  thumb,  a  multilayer  should  contain  more  than  (600A  /  X  )* 
pairs  of  layers  in  order  to  achieve  absorption-limited  reflectivity  at 
normal  incidence;  see  sec.  XI-3.) 

In  the  above  experiments,  the  deposited  materials  were  transferred  to 
the  substrates  by  vacuum  evaporation.  Barbee  and  Keith  (1979) •  and 
Barbee  (1982),  report  the  fabrication  of  tungsten/carbon  multilayers  with 
periods  as  short  as  15A  using  a  vacuim  sputtering  technique. 

Sophisticated  deposition  apparatus  and  procedures  aH*ar  to  be  the  key  to 
their  achievement  of  high  quality  diffracting  structures. 

Spiller  et  al.  (1980)  report  the  i*q>lementation  of  an  in-situ 
reflectance  monitoring  technique  in  which  the  interference  oscillations  of 
•  reflected  x-ray  probe  beam  are  monitored  during  multilayer  fabrication. 
Such  a  system  directly  monitors  the  phasing  between  the  partially 
deposited  qpper  layer  and  the  preoeeding  stack  (in  contrast  to  the  usual 
crystal  oscillator  which  monitors  the  (mass)  thickness  of  the  upper  layer 
alone).  In  essence,  the  reflectance  monitoring  technique  allows  the 
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reflection  axponent  from  the  upper  layer  to  be  brought  into  phase  with 
the  pceceeding  stack,  despite  any  thickness  errors  that  this  stack  any 
contain. 

At  present  Spiller's  prop  at  XBM  and  Barbee's  group  at  Stanford  aay 
be  regarded  as  the  dominant  centers  for  research  in  this  field.  X-ray 
multilayers  have  also  recently  been  fabricated  by  Soviet  researchers  using 
laser  evaporation  of  the  materials  bo  be  deposited  (Gaponov  et  si.,  1981). 

The  most  commonly  used  materials  in  multilayers  to  date  are  carbon 
for  the  low  index  layer  and  tungsten  or  rhenium/tungsten  alloy  for  the 
high  index  layer.  The  optical  constants  of  rhenius  (Z  ■  75)  and  tungsten 
(Z  •  74)  are  generally  comparable,  so  the  nominal  performance  of  the  two 
materials  should  be  similar  in  the  absence  of  defects.  We  will  generally 
use  tungsten/carbon  multilayers  when  presenting  nisnerical  examples.  In 
doing  so  we  use  preliminary  versions  of  the  optical  constants  data  in 
Henke  et  al.  (1982) ,  generously  provided  to  us  by  the  authors. 

Vinogradov  and  Zeldovich  (1977)  present  an  analytic  treatment  of 
reflection  from  periodic  x-ray  multilayers.  They  use  (normal  incidence) 
coupled  wave  equations  whose  unknowns  are  the  envelope  functions  of 
high-frequency  electromagnetic  field  oscillations.  (The  "carrier 
frequency”  is  taken  to  be  the  spatial  periodicity  of  the  multilayer 
structure.)  Such  coupled-wave  equations  apply  to  multilayer  reflectivity 
when  the  coupling  constant  is  weak,  as  it  is  in  the  x-ray  case. 

Vinogradov  and  Zeldovich  obtain  expressions  for  the  reflectivity  of  s 
periodic  multilayer  containing  an  arbitrary  lumber  of  layers.  They  also 
determine  ttrn  condition  which  optimizes  the  ratio  in  thickness  of  the  two 


different  layers,  and  finally  they  derive  what  is  in  effect  a  hybrid 
optimization  condition  that  gives  the  aptlansn  total  thickness  of  the 
period,  when  the  thickness  ratio  is  also  optimized.  (The  two  relations 
together  thus  implicitly  give  a  refractive  correction  to  Bragg's  law.) 


Part  B)  Crystal  Diffraction 

Zn  aany  respects,  raflaction  from  x-ray  multilayers  is  similar  to 
diffraction  from  crystals  in  tha  Bragg  node  (atonic  plants  parallel  to  the 
crystal  surface) ;  this  is  particularly  true  in  the  case  of  periodic 
multilayers. 

Most  crystals  have  period  lengths  (or  "d-spacings")  that  are  a  few 
angstroms  or  less  between  adjacent  diffracting  planes,  so  that  the  longest 
wavelength  that  can  be  reflected  from  most  crystals  is  less  than  ten 
angstroms.  (Itie  longest  wavelength  that  can  be  reflected  is  approximately 
twice  the  period  length,  which  is  referred  to  as  the  “2d- spacing*,  the 
basic  period  of  a  crystal  (as  a  physical  entity)  is  usually  referred  to  as 
the  “unit  cell".) 

The  optical  constants  of  materials  in  the  soft  x-ray  region 
(5A  -  125A)  and  the  extreme  vacuum  ultraviolet  (125A  -  600A)  have  a 
qualitative  similarity  to  the  optical  constants  of  materials  at 
conventional  x-ray  wavelengths  (0.5A  to  SA) . 

In  the  high  frequency  (but  non-relativistic)  regime  where  the 
radiation-induced  electron  oscillations  are  limited  primarily  by  electron 
inertia  rather  than  by  coupling  with  the  atom  (“free-electron  regime"), 
the  real  part  of  the  index  of  refraction  (less  unity)  scales  as  X2  , 
while  the  imaginary  part  scales  as  X4  (Heitler,  1954,  p.208),  so  that  at 
short  x-ray  wavelengths  absorption  becomes  snail  oanpared  to  dispersion. 

Since  the  real  part  of  the  index  of  refraction  is  less  than  one  for 
electromagnetic  frequencies  above  the  plaana  frequency,  "total  external 
reflection"  of  x-rays  can  occur  when  the  angle  of  incidence  to  a  single 


reflecting  surface  is  less  than  the  critical  angle  given  by 
%  m  J  t  ( 1 -  Re  0»))  radians,  with  n  the  index  of  refraction  of  the 
substrate  (Janes,  1965,  p.171).  X-ray  reflectivities  that  approach  unity 
can  therefore  be  obtained  at  glancing  angles  (typically  —  89*  to  the 
normal)  when  absorption  is  snail.  Total  external  reflection  of  short 
wavelength  x-rays  can  therefore  be  quite  efficient. 

In  the  soft  x-ray  region  the  imaginary  part  of  the  index  of 
refraction  can  became  comparable  in  magnitude  to  the  real  part  (less 
unity).  The  effective  number  of  free  electrons  per  atom  (known  as  the 
'atomic  scattering  factor”)  Bust  then  be  treated  as  a  complex  quantity 
whose  dependence  on  wavelength  is  irregular  due  to  atomic  resonances. 
(However,  it  appears  from  the  data  in  (Henke,  et  al. ,  1982)  that  for  meet 
materials  the  imaginary  part  of  the  index  scales  very  broadly  as  X  in 
the  soft  x-ray  region.) 

The  mo6t  rigorous  theory  of  diffraction  from  crystals  is  generally 
considered  to  be  the  dynamical  theory  of  Bwald  and  von  Laue  (Batterman  and 
Cole,  1964;  James,  1965,  p.413) .  In  the  Ewald  -  von  Laue  theory  the 
physical  structure  of  the  crystal  is  represented  by  s  complex  dielectric 
constant  that  is  spatially  periodic,  and  that  has  as  interaction  with 
radiation  given  by  Maxwell's  material  equations. 

The  dynamical  theory  can  be  used  to  treat  multiple  reflection 
processes  in  which  an  incident  beam  interacts  resonantly  with  more  than 
one  set  of  atomic  planes  at  once,  but  in  the  most  carman  case  only  one 
pair  of  structural  spatial  frequency  earponents  (positive  and  negative) , 
plus  a  DC  component,  are  considered  to  be  active  at  one  time. 


X-l-8 


Fran  a  fully  classical  point  of  view,  the  polarizability  can  bt  taken 
to  be  proportional  to  a  local  electron  density,  but  aider  a  sanl  classical 
•odel  the  polarizability  Must  be  calculated  quantize  Mechanically.  At 
short  wavelengths  where  the  dipole  approximation  breaks  down,  the 
dielectric  constant  will  therefore  depend  on  the  scattering  angle.  This 
effect  is  anall  if  the  2d- spacing  is  larger  than  an  individual  a tan. 

When  the  surface  of  a  thick  crystal  is  cut  parallel  to  the 
diffracting  planes,  the  reflectivity  predicted  by  the  Bwald  -  von  Laue 
thoery  is  the  sane  as  that  obtained  with  the  earlier  theory  of  Darwin  and 
Prins  (Janes,  1965,  p.429).  Zn  essence,  the  Darwin-Prins  solution  is 
obtained  by  requiring  that  the  change  in  the  forward  and 
backward-traveling  amplitudes,  as  each  cell  is  traversed,  be  consistent 
with  the  cell  reflectance  and  transmittance. 
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part  C)  Thin  Film  Formalins 

Kogelnik  (1976)  has  derived  (from  coupled  wave  equations)  a 
differential  equation  that  prcpagates  the  wrplitude  reflectanoe  through  a 
quasi-per iodic  structure,  this  differential  equation  is  similar  to  the 
difference  equation  with  which  we  treat  x-ray  multilayer  reflectivity  in 
Chapter  XI. 

He  derive  our  difference  equation  frcn  a  characteristic  matrix 
solution  for  the  fields  within  the  unit  cell  of  an  x-ray  multilayer;  this 
unit  aell  may  contain  an  arbitrary  one-dimensional  index  gradient.  Such  a 
solution  exists  (in  principle)  for  any  median  in  which  the  index  of 
refraction  varies  only  in  one  dimension  (Born  and  Wolf,  1976,  p.51).  The 
characteristic  matrix  solution  for  homogeneous  layers  is  cotnonly  used  in 
optical  thin-film  calculations. 

Although  this  work  centers  on  the  analysis  of  structures  with  a 
one-dimensional  index  variation,  we  summarize  in  Chapter  XV  the  results  of 
a  preliminary  analysis  w>  have  made  of  the  effect  of  interfacial  roughness 
on  multilayer  reflectivity. 

This  analysis  is  based  on  that  of  Eastman  (1978) ,  who  showed  that  a 
one-dimensional  formalism  can  be  used  to  treat  Interfacial  roughness  when 
the  roughness  has  a  sufficiently  gradual  variation  within  the  layers. 
Eastman’s  methods  have  been  summarized  and  extended  by  Carniglia  (1981). 

Shellan  et  al.  (1978)  use  a  coupled  wave  formalism  to  make  a 
perturbation  treatment  of  the  effect  of  randan  thickness  errors  in 
dielectric  multilayers.  Our  non-perturbative  analysis  of  randan  thickness 
errors  in  *-rsy  multilayers  is  presented  in  sec.  XX-5. 


_  i 


Section  1-2  -  Outline  of  Text 


This  work  contains  four  chapters,  each  of  which  is  divided  into 
sections,  and  parts  of  sections.  The  present  chapter  reviews  the  field. 
Chapter  II  presents  our  theoretical  analysis  of  x-ray  multilayer 
reflectors.  Chapter  III  uses  this  theory  to  discuss  applications  for  x-ray 
multilayers,  and  Chapter  IV  stsamarizes  cur  results. 

The  text  uses  a  common  amsaeratian  system  to  number  p tapes,  aquations, 
figures  and  tables.  M  an  example,  the  first  equation  of  section  1  in 
Chapter  II  is  lumbered  "eq.  II-l-l",  the  second  is  eq.  II-1-2,  and  so  on 
through  the  final  equation  of  the  section  (eq.  1 1-1-25! .  The  first 
equation  of  the  next  section  is  then  eq.  1 1 -2-1,  and  so  forth.  The  same 
system  is  used  to  nunber  pages,  figures,  and  tables. 

In  sec.  II-l  we  introduce  the  basic  formalism  for  our  treatment  of 
reflection  from  x-ray  multilayers.  Like  the  Ewald  -  von  Laue  theory  of 
diffraction  from  crystals,  we  begin  with  a  physical  description  of  the 
multilayer  in  terms  of  a  spatially  varying  oonplex  dielectric  constant. 

Ne  assisne,  as  is  nstural  for  multilayers,  that  the  dielectric  constant 
varies  only  in  one  dimension.  This  ass  nation  has  an  analogy  with  the 
Bweld  -  sen  Laue  theory,  in  that  we,  in  effect,  neglect  s  portion  of  that 
set  of  spatial  frequencies  which  are  alao  neglected  in  the 
Ewsld  -  van  Laue  theory  when  only  one  reflection  is  regarded  as  being 
active  (sac  above);  for  axm^le,  In  the  case  of  a  periodic  multilayer,  we 
in  effect,  neglect  «iy  Bpatiel  frequencies  that  are  not  parallel  to  the 
active  pair  of  spatial  frequencies  (assumed  to  be  normal  to  the 
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substrate) . 

However,  unlike  the  dynamical  theory,  we  do  not  require  that  the 
■ulti layer  be  periodic;  in  effect  we  allow  each  of  the  active  spatial 
frequencies  to  be  replaoed  by  a  narrow  continuum  of  spatial  frequencies. 

Our  formalism  does  not  introduoe  such  a  continues  explicitly,  but 
instead  makes  direct  use  of  the  structural  properties  of  the  discrete 
cells.  Ibis  is  accomplished  with  a  difference  equation  that  propagates 
the  amplitude  reflectance  from  cell  to  cell  (derived  in  sec.  II— 1) . 

In  sec.  11-2  we  consider  the  reflectivity  of  ideal,  defect-free 
multilayers;  our  main  interest  is  in  the  optimization  of  the  layer 
thicknesses  and  materials. 

Sec.  II-3  considers  the  angular  and  spectral  scaling  of  the 
reflection  properties  of  tungsten/carbon  multilayers. 

Sec.  11-4  considers  the  possible  utility  of  aperiodic  reflector 
designs,  with  particular  emphasis  on  the  layer  by  layer  optimization 
scheme  of  Carniglia  and  hpfel  (1980).  The  in-si tu  reflectance  monitoring 
technique  reported  by  Spiller  et  al.  (1980)  lends  itself  naturally  to  this 
kind  of  design  procedure. 

Sec.  XI-5  considers  the  effect  of  random  errors  in  the  layer 
thicknesses.  We  make  the  problem  tractable  by  neglecting  higher  pcwers  in 
a  quantity  that  may  be  regarded  as  the  incoherent  reflectivity  of  the 
degraded  structure.  In  this  way  the  analysis  is  able  to  treat  a  full 
range  of  error  magnitudes,  from  those  which  cause  only  a  slight  decrease 
in  reflectivity,  to  those  large  enough  to  substantially  degrade  the 
reflectivity. 
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Spiller  et  al.  (1980)  have  verified  experimentally  that  in-aitu 
reflectance  monitoring  strongly  reduces  the  degradation  in  reflectivity 
caused  by  thickness  errors.  We  model  this  effect  quantitatively  in 
sec.  11-6.  However,  due  to  the  complicated  interaction  involved  in  the 
monitoring  process,  the  results  of  sec.  II-8  should  probably  not  be 
regarded  as  strongly  quantitative  predictions. 

In  the  first  section  of  Chapter  m  we  discuss  points  of  general 
relevance  concerning  the  application  of  x-ray  multilayers.  In  sec.  II 1-2 
we  discuss  the  possibility  of  using  x-ray  multilayers  to  form  cavity 
mirrors  for  projected  x-ray  lasers  ( X  ~  100A) .  Sec.  III-3  considers  the 
utility  of  multilayer  coatings  in  short  wavelength  x-ray  microscopes 
(  X  ~  1A) . 

The  first  section  of  Chapter  IV  svznnarizes  the  results  of  the 
preceding  chapters.  Sec.  TV-2  discusses  possible  avenues  for  future 
research,  and  also  simnar ires  a  preliminary  investigation  of  the  effect  of 
inter facial  roughness  on  multilayer  reflectivity.  In  terms  of  the  above 
discussion,  the  analysis  of  sec.  IV- 2  may  be  said  to  permit  the  continuum 
of  spatial  frequencies  in  the  multilayer  structure  to  have  a  narrow  extent 
in  the  non-normal  direction. 


Chapter  II  theoretical  Analysis  of  Multilayer  Reflection 


Section  II-l  -  Analytic  Formalism 

Part  A)  Introduction 

Me  have  seen  in  the  preceeding  chapter  that  the  dynamical  theory  of 
x-ray  diffraction  is  based  on  the  assumption  of  a  spatially-varying,  complex 
dielectric  constant  which  has  an  interaction  with  radiation  defined  by 
Maxwell's  material  equations.  Despite  the  seemingly  ad-hoc  nature  of  this 
assumption,  the  dynamical  theory  is  generally  considered  to  be  the  most 
rigorous  theory  of  x-ray  diffraction,  and  has  been  found  to  yield  fairly 
reliable  results  with  crystalline  structures  having  angstrom  periodicities 
(Batterman  and  Cole  1964) . 

Because  such  assumptions  permit  the  use  of  classical  techniques  of 
electromagnetic  theory,  we  will  also  use  them  in  our  treatment  of  x-ray 
multilayer  reflectors.  However,  in  our  case  the  effective  dielectric 
constant  need  not  be  assisted  to  be  periodic,  but  must,  on  the  other  hand,  be 
assumed  to  be  a  purely  one-dimensional  function  varying  only  in  the 
direction  perpendicular  to  the  multilayer  surface. 

Our  formalism  will  tie  in  closely  with  traditional  methods  of  thin- film 
optics.  It  will  be  applicable  to  a  broader  range  of  problems  involving 
x-ray  multilayers  than  are  other  x-ray  formalisms,  and  will  provide  analytic 
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solutions  in  the  x-ray  regime  to  problms  that  require  asnbersame  mnerical 
calculations  if  traditional  thin- film  Methods  are  used. 

Our  formalism  will  contain  the  structural  properties  of  the  different 
cells  explicitly,  whereas  the  Ewald-von  Laue  or  coupled-wave  formalisms  are 
based  cn  the  Fourier  components  of  the  diffracting  structures.  Fourier 
exponents  are  somewhat  inconvenient  to  use  when  analyzing  aperiodic 
multilayers. 

An  additional  advantage  of  our  formalism  is  that  it  proceeds  via  a 
well-characterized  series  of  assumptions  from  standard  results  of 
electromagnetic  theory.  For  this  reason  our  solutions  include  eertain 
higher-order  terms  that  must  be  included  in  order  to  obtain  even  a 
lowest-order  analysis  of  certain  problems  of  interest.  Such  problems 
include  the  determination  of  the  layer  thicknesses  in  multilayers  made  by 
in-situ  reflectance  monitoring  (sec.  11-4)  and  the  determination  of  the 
effect  on  reflectivity  of  random  thickness  errors  in  multilayer  structures 
(sec.  II— 5) . 


( 
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Part  B)  Notation,  Preliminary  Result s,  and  Characterization  of 
Approximations 

The  core  of  our  formalism  is  an  analytic  solution  for  the 
electromagnetic  fields  inside  the  unit  aell  of  an  x-ray  multilayer 
that  has  an  arbitrary  gradient  structure.  In  the  optical  regime,  such  a 
solution  cannot  be  obtained  analytically  with  arbitrary  gradients  (see 
below) ;  we  are  able  to  obtain  analytic  solutions  in  the  x-ray  case  by 
treating  the  decrement  in  the  dielectric  constant  as  a  snail  quantity. 

It  has  been  shown  (Born  and  Wolf,  1975,  p.51)  that  when  the  dielectric 
constant  has  a  purely  cne-dimemensicnal  spatial  variation  along  a  direction 
z,  Maxwell's  partial  differential  equations  can  be  separated  into  ordinary 
differential  equations  involving  variables  D(z)  and  V(z),  where  U  and  V  are 
defined  in  terms  of  the  electromagnetic  fields  through  the  relations 


H*  * 


for  the  P  case  and 


uw,,k.!i" 

» sin  0  y 
i  lr#  sin  By 


EM  s  U  (*)e 
Hy  *  V  (*)e 


(i-i-i) 


(1-1-2) 


for  the  S  case.  (We  define  the  P  case  to  be  that  in  which  the  magnetic 
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fitld  H  points  slang  the  x-axis,  and  the  S  case  to  be  that  in  which  the  E 
field  points  along  x.) 

Here  the  parameter  k^  is  defined  by 

(1-1-3) 


with  X.  the  vacuum  wavelength.  0  is  the  vacuun  angle  of  incidence  to  the 
normal.  In  most  cases  we  will  folia*  the  usual  optics  convention  where 
angles  of  incidence  are  specified  to  the  normal. 

The  exponents  in  the  y  oscillation  factors  have  been  put  in  a 
trigonometric  form  that  enables  them  to  be  fit  to  boundary  conditions  that 
are  independent  of  x  and  that  apply  in  planes  of  constant  2,  i.e.  our 
solutions  are  aj^ropriate  to  a  Bragg  geometry  in  which  the  y-2  plane  is  the 
plane  of  incidence  (See  fig.  11-1-1). 

The  P  case  is  somewhat  more  interesting  than  the  S  case,  and  is 
disscussed  in  somewhat  less  detail  by  Born  and  Wolf  (1975,  p.54) .  For  this 
reason  we  will  briefly  describe  the  steps  by  which  the  solutions  for  the 
fields  in  P  polarization  are  put  into  characteristic  matrix  form. 

As  is  shown  by  Born  and  Wolf,  the  wnplitudes  u  and  v  satisfy  the 
ordinary  differential  equations 
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and 
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k.  €  (*) 


dU 

dz 


(i-j-j) 


in  the  P  case. 

We  consider  next  the  unit  aell  of  an  x-ray  multilayer  reflector  that 
extends  from  z-zt  to  *»zt  (with  z4  ^  zA  ).  As  indicated  in 
fig.  II -1-1  we  will  define  the  direction  of  increasing  t  to  be  towards  the 
substrate  from  the  incident  side.  Let  (z)  ,V4  (z);  Ux  (2)  ,V,  (2)  be  two 
sets  of  solutions  for  the  transverse  fields  and  E y  that  satisfy 


Ui  (*i)  *  0,  Vi  (*i)  *=  1 
(*i)  =  1,  V*  (*x)  =  0 


U-l-6) 


Then  following  well-known  procedures  for  linear  ordinary  differential 
equations,  we  have  the  following  unique  relation  connecting  general  fields 
at  2^  with  those  at  Zj  : 

/-E^n  /  v,  (*,)  v,  r*,>\  /-e,  c.n 

V  H,  (*,) }  l  U,  (*,)  U,  (**>./  \  H,  (*t) ) 


'xrS  — 
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(The  minus  signs  in  the  field  vectors  are  a  result  of  the  definitions  in 
•q.  11-1-1) .  This  linear  relation  obtains  because  aqs.  11-1-4,5  are 
homogeneous . 

The  determinant  of  the  above  matrix  must  be  one  according  to  the  ss*e 
mathematical  arguments  as  are  used  in  the  S  case  by  Born  and  Wolf  (1975,  p.56) . 
(Note  that  their  supplemental  physical  arguments  concerning  conservation  of 
energy  cannot  be  applied  in  the  case  of  an  absorbing  structure) . 

Since  the  determinant  is  unity  w  can  write  the  following  inverse  to 
eq.  11-1-7: 


"  By  (Z,Y 


-VXUX)\  /-E,  UAV 


or,  in  more  usual  form: 


Uilz^  \  H, 


(l-l-F) 


h  V-x  (*J  \  /  E,  W 


,hx  (*4) J  VujC**)  vtW/ 


(a-i-3) 


This  is  what  Born  and  Wolf  call  the  characteristic  matrix  solution  for  the 
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fields  in  •  stratified  Medium;  in  the  case  of  a  unit  aell  in  an  x-ray 
reflector  it  expresses  the  fields  an  the  incident  side  of  the  cell  in  terms 
of  those  cn  the  substrate  side.  In  order  to  more  conveniently  treat 
reflection  at  non-normal  angles  of  incidence  we  will  refer  to  a  slightly 
different  Matrix  as  the  charcteristic  matrix  (eg.  11-1-14  below).  In  this 
way  we  trill  automatically  incorporate  in  our  equations  certain  well-known 
rules  of  thin-film  optics  that  reduoe  the  case  of  off -normal- incidence 
reflection  to  the  normal  incidence  case  (Baumeister,  1962). 

We  will  subscript  all  quantitities  that  appear  on  the  right  side  of  the 
characteristic  matrix  equation  with  the  oontnun  subscript  K.  Thus  the  fields 
at  the  substrate  interface  of  a  particular  cell  in  a  multilayer  will  have 
the  same  subscripts  as  the  structural  parameters  of  the  cell,  while  the 
subscripts  of  the  fields  at  the  incident  interface  will  be  larger  by  one. 

Generally,  we  will  use  the  indicial  label  K  to  refer  to  an  arbitrary 
cell  in  a  multilayer,  and  the  label  J  to  refer  to  the  last  or  total  nutiber 
of  cells.  Thus,  K  is  equal  to  1  at  the  substrate  and  rises  to  J+l  at  the 
upper  interface  (see  fig.  II-l-l).  (Note  that  given  the  matrix  inversion 
performed  between  eqs.  II-1-7  and  8,  we  have  chosen  to  have  K  increase  in 
the  direction  of  decreasing  z.) 


f 
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Bg.  II-1-4  can  aasily  be  transformed  into  a  Schrodinger  equation  by 
eliminating  the  J  U/ dz  term.  Since  the  Schrodinger  equation  con  only 

be  solved  analytically  in  a  few  isolated  oases  (Killians,  1M2) ,  eq.  11-1-4 
cannot  be  solved  analytically  in  the  ease  of  an  arbitrary  gradient. 

However,  in  Appendix  1,  we  show  how  straightforward  perturbation 
techniques  can  be  used  to  solve  eqs.  11-1-4,5  under  the  as sanction  that  the 
parameter  A  ,  defined  by 


tk  (z)  *  ~  (e  (*)  - 1 )  (jl-i-k >; 


is  small  compared  to  one.  Since  A  may  be  thought  of  as  the  unit  decrement 
to  the  index  of  refraction,  this  perturbation  solution  is  appropriate  to  the 
x-ray  regime  where  the  indieies  of  refraction  of  all  materials  approach 
unity. 

The  aell  structure  is  considered  to  be  specified  by  the  function  A (z) 
in  the  range  za  £  z  4=.  *fc  • 

Without  loss  of  generality  we  can  position  the  z*D  origin  at  the 
midpoint  of  the  cell.  He  then  specify  the  oell  thickness  in  terms  of  a 
dimensionless  Bragg  detuning  parameter  ^  : 


f  m  n%  -  ka  cos  &  (zt-zj  (l-i-11) 

z.  w  7t  -  cos  6  d 


with  d  the  oell  thickness. 
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Xn  aq.  XX >1-11  r  m  is  the  order  of  dif friction.  Under  the 
approximations  to  which  mb  will  work  (see  below) ,  s  resonance  in  mth  order 
may  unambiguously  be  defined  as  the  situation  in  which  9  i»  wall* 

A  single  unit  aell  can  be  treated  with  a  perturbation  approach  because 
the  interaction  of  the  single  cell  with  the  Incident  x-ray  bean  nust  be 
weak. 

We  will  generally  be  interested  in  multilayers  operating  sufficiently 
near  resonance  that  the  reflectivity  of  the  total  structure  is  significantly 
higher  than  the  reflectivity  of  the  individual  cells.  For  this  reason  it 
will  generally  not  be  appropriate  to  treat  the  entire  structure  using  a 
perturbation  approach. 

As  an  alternative,  we  oould,  of  course,  obtain  the  character istic 
matrix  solution  for  the  entire  structure  by  multiplying  the  matrix  solutions 
for  the  ireSividual  cells  together.  This  standard  procedure  is  based  on  the 
continuity  of  the  transverse  field  opponents  across  the  cell  interfaces 
Such  a  procedure  may  be  thought  of  as  a  numerical  integration  of 
«qs.  11-1-4,5  with  the  thickness  of  each  unit  cell  serving  as  the  step  size 
(thus  the  requirement  that  the  change  introduced  by  each  unit  cell  be 
•mall) . 

Rather  than  sultiplying  the  individual  oell  matricies  together  to  find 
the  overall  solution,  we  will,  instead,  derive  from  the  individual  cell 
solution  a  difference  equation  that  propagates  the  «rplitude  reflectivity 
from  cell  to  aell.  The  Imown  (and  usually  negligable)  reflectivity  of  the 
substrate  then  serves  as  s  boundary  condition  from  which  to  determine  the 
reflectivity  of  the  entire  structure. 


-k. 
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It*  derivations  of  the  avatrix  aolutien  and  difference  equation  are 
presented  in  A^endicies  1  and  2;  before  presenting  these  results  we  will 
first  daaracterize  wore  precisely  the  approximations  used. 

In  terms  of  our  notation,  the  usual  theories  of  x-ray  diffraction 
generally  include  oily  first  order  terms  in  the  decrement  A  and  the 
detuning  parameter  f  .  Neglect  of  high  orders  in  A  is  generally  well 
justified  in  the  re-ray  regime. 

Neglect  of  higher  orders  in  ^  is  justified  when  the  structure  is 
periodic  and  operating  near  resonance  (near  enough  that  the  reflectivity  of 
the  overall  structure  is  significantly  higher  than  that  of  the  individual 
cell) . 

This  may  be  seen  by  the  following  argument.  The  number  of  cells 
participating  in  the  reflection  can  be  at  most  of  order  1/A  due  to 
extinction  of  the  beam.  Here  A  may  be  considered  to  be  an  average  of  the 
decrement  over  typical  values  for  the  eell.  (If  extinction  were  due  only  to 
absorption,  we  could  consider  oily  the  imaginary  part  of  A  here,  but  at 
shorter  wavelengths,  depletion  by  back  reflection  is  also  significant.) 

Success  in  obtaining  high  reflectivity  will  depend  on  the  total  phase 
detuning  in  the  stack  of  participating  oells;  this  detuning  will  be  of  order 
$/&  since  <p  is  essentially  the  detuning  per  cell.  For  operation  near 
resonance,  the  total  detuning  cannot  be  large  compared  to  cane  radian,  so 
that  in  cases  of  greatest  interest,  ?/A  fc  1.  Given  that  A  1 .  a-, 
expansion  that  is  first  order  in  f  as  well  as  in  A  is  adequate. 

Higher  order  terms  are  necessary  to  treat  kinds  of  aperiodicaties. 


t 
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The  first  are  aperiodicities  of  a  stochastic  nature;  In  auch  eases  a 
randcn  error  in  at  cne  point  in  the  structure  night  be  canceled  out  by 
that  at  another,  so  that  larger  RMS  detunings  within  the  individual  cells 
can  be  tolerated  than  is  possible  in  the  periodic  case.  Specifically,  since 
we  would  expect  the  total  RMS  phase  error  front  the  3  cells  that  participate 
in  the  reflection  to  be  of  order  ^  ^  J  ^  ,  which  is  of  order 

[<**>/a]*  •  we  would  expect  that  in  the  cases  of  greatest  interest: 


A 


(1-1-12) 


In  fact,  we  will  see  in  sec.  11-5  that  if  the  expectation  value  of  is 
equal  to  its  unperturbed  value  (i.e.  its  value  in  the  periodic  case),  the 
lowest  order  effect  of  randan  errors  in  ^  on  the  reflectivity  is  through 
terms  of  order  ?  . 

The  second  kind  of  multilayer  structure  whose  performance  is  strongly 
dependent  on  higher  order  terms  is  that  which  results  when  irrsitu 
reflectance  monitoring  is  used.  In  such  multilayers  the  ^  value  for  the 
Kth  pair  of  layers  is  intended  to  be  that  value  which  maximizes  the 
reflectivity  RK44  of  the  entire  stack  of  cells  1  ...  K  (see  sec.  JI-4). 

Determination  of  this  ?  value  mnounts  to  setting  a  derivative  of 
intensity  reflectance  with  respect  to  equal  to  zero.  We  will  see  in 
sec.  11-4  that  the  lowest  order  terms  in  such  a  derivative  are  of  order  A 
(i.e.  A1  rather  than  A°  )•  Terms  of  order  unity  are  missing  because 


( 
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the  lowest  order  effect  of  a  change  in  f  having  magnitude  S  •  is  to 
multiply  the  reflected  arplitude  by  exp (-2i  SfK )  leaving  the  intensity 
reflectance  unchanged.  To  obtain  terms  of  order  A  in  a  differentiation 
with  respect  bo  f  ,  terms  of  order  ?  •  A  must  be  included  initially. 

In  order  to  treat  both  kinds  of  aperiodicity,  our  analysis  will  include 
all  orders  in  the  parameter  9  »  and  also  terms  of  order  9 '  A  • 

As  in  the  traditional  x-ray  diffraction  theories,  we  will  include  only 
first  order  terms  in  A  .  Such  an  approximation  breaks  down  at  grazing 
angles  of  incidence.  As  the  angle  approaches  the  surface  the  path  traversed 
through  each  cell  increases  for  two  reasons;  first  because  the  x-ray  photons 
traverse  the  cell  at  more  oblique  angles,  and  second  because  the  cell 
thicknesses  have  to  be  increased  as  the  angle  decreases  in  order  to  maintain 
a  resonant  reflection  from  the  multilayer.  These  two  factors  imply  a 
quadratic  increase  in  extinction  per  layer  with  angle;  therefore  when 

(f-e)*  <  |a  |  d-i-13) 


a  regime  has  been  reached  in  which  the  beam  has  a  strong  interaction  within 
one  single  cell,  and  our  formalism  is  no  longer  valid. 

Our  formalism  will  depart  from  the  common  practice  in  x-ray  diffraction 
theory  in  that  we  will  include  the  trigonometrical  factors  that  are 
necessary  for  the  theory  to  be  applied  near  normal  incidence. 


i 


Xn  the  reminder  of  the  text  we  will  follow  the  standard  practice  of 


using  primes  and  double  primes  to  denote  real  and  imaginary  parts, 
respectively. 
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Part  C)  Characteristic  Matrix  Solution  for  the  X-Ray  Regime 
In  Appendix  1  we  show  that  subject  to  the  approximations  discussed 
above,  the  solution  for  the  P  polarisation  field  solitudes  04  2  ,  Vi  2  is 


where 


j  (-i)mcos  tx-P* 


ICM)"1  5in  tK  4  rK) 


(l-i-14) 


k  P  (B)  I 

rKB  -  J  dz  A(t)cos(ik#to5  6x) 
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k  P  (&)  f 
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r  cos  ZB  P  Polarization 

P(B)  *  1 


1 


S  Polarization 
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Bq.  XI-1-14  represents  •  straightforward  perturbation  solution  to 
aq.  II-1-4;  terns  of  order  A*  are  neglected,  and  the  vaeutse  (i.e.  A  ■  0) 
solution  for  0  is  substituted  into  all  terns  that  are  first  order  in  A* 
Bq.  II-1-4  then  becanes  a  driven  harmonic  oscillator  equation. 

In  eq.  XI-1-14,  i/  and  are  the  field  oenponenta  that  are  normal  to 
the  k-vector  in  vicuna;  from  fig.  XX-1-2  we  see  that  these  are  defined  by 


t  E  -  *  /  cor  6 

*  ■ 


(Again  we  note  that  0  is  the  vacuum  angle  of  incidence) .  Thus  if  we  follow 
the  usual  sign  convention  of  thin-film  optics  in  which  the  wplitude 
reflectivity  £  is  real  and  positive  if  the  reflected  zy  is  in  phase  with 
the  incident  (for  P  polarisation),  we  see  from  fig.  II-1-2: 


t  *  (14  f)  A 
*  -  <1-  ?)  A 


(1*1-  17) 


where  A  is  some  K-dependent  field  wplitude. 

He  show  in  Appendix  1  that  eq.  XI-1-14  is  also  valid  for  the  case  of 
S  polarisation  (in  which  E  ia  polarised  along  the  x-axia)  if  we  define 


V  E  E* 


(1-1-U) 
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COORDINATE  SYSTEM  FOR 
CHARACTERISTIC  MATRIX  SOLUTION 


In  all  figures  p  is  real  and  positive.  Field  vectors  are  shown  as  solid 
arrows,  propagation  directions  as  dotted  arrows. 

•  Figure  a  shows  S  polarization  at  normal  Incidence. 

•  Figure  b  shows  P  polarization  at  normal  Incidence.  The 

coordinate  system  of  fig.  a  has  been  rotated  90° 
relative  to  the  field  vectors. 

•  Figure  c  shows  S  polarization  at  non-normal  Incidence. 

•  Figure  d  shows  P  polarization  at  non-normal  Incidence. 
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COORDINATE  SYSTEM  FOR 
CHARACTERISTIC  MATRIX  SOLUTION 
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Part  D)  Difference  Equation  for  Amplitude  Reflectance 

As  disaissed  in  sec.  n-l-B  above.  choose  to  convert  the  matrix 
solution  for  the  unit  cell  to  a  difference  aquation  in  ^  before  proceeding 
to  consider  the  reflectivity  of  the  multilayer  as  a  whole. 

The  derivation  in  Appendix  2  follows  a  well-known  procedure,  except 
that  in  the  course  of  the  derivation  we  make  truncated  expansions  consistent 
with  the  approximations  discussed  in  see.  2I-1-B  above. 

In  brief,  eq.  II-1-17  is  substituted  into  aq.  II-1-14,  and  the 
anplitude  A  is  cancelled  from  the  resulting  pair  of  aquations.  Terms  of 
order  A*  are  then  neglected. 

The  resulting  equation  that  propagates  the  amplitude  reflectivity 
across  the  Kth  cell  is: 

(1-J-2C) 

Eq.  II-1-20  is  analogous  to  the  well-kncam  Airy  recursion  formula  that 
is  used  to  propagate  the  reflectivity  frcrn  single-layer  to  single- layer  in 
optical  multilayers. 

In  fact,  if  we  near  employ  a  phenomenological  Airy  argument  based  on  the 
usual  smmation  of  partial  reflection  and  transmission  ccnponents  within  the 
aingle-oell,  we  can  physically  interpret  eq.  II-1-20  by  comparing  it  with 
the  result  of  the  Airy  sumaticn. 

The  Airy  autmaticn  argument  is  summarized  in  fig.  II- 1-3;  the 
phenomenological  derivation  of  the  partial  reflection  and  transnission 
caponents  has  been  omitted,  but  argunents  that  are  essentially  the  same  are 
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used  in  the  Darwin-Prins  theory  of  *-ray  diffraction  (see  Janes,  1965,  p.62). 
As  shewn  in  the  figure,  the  first  ter*  an  the  right  side  of  eq.  II-1-20 
represents  s  double  transmission  through  cell  K  with  an  interaediate 
reflection  from  the  stack  of  preoeeding  oells  F-l,  K-2,  ...  1.  second 
term  represents  a  reflection  from  oell  F  back  into  the  direction  of  the 
incident  medium.  The  third  term  represents  a  multiple  reflection  process 
consisting  of  reflection  from  the  preoeeding  stack  of  oells,  reflection  back 
tevards  the  substrate  from  oell  K,  and  finally  a  second  reflection  from  the 
preoeeding  stack  back  into  the  direction  of  the  incident  medium.  Higher 
order  nultiple  reflections  are  of  order  A*  and  are  neglected. 

The  comparison  thus  implies  (as  can  be  shown  rigorously  from  the  matrix 
solution  of  eq.  II-1-14)  that  the  reflectance  of  the  Kth  oell  is 
(-1)"  (iTjj  -p^e*’**  from  the  incident  side  and  (-1)"1  P,  from 

the  substrate  side,  and  that  the  transmittance  of  the  cell  is  -e  ~  . 

The  analogy  between  eq.  II-1-20  and  the  single-layer  Airy  formula  can 
be  formalised  in  another  way  if  the  cell  is  centx osymetr ic ,  i.e.  if  it 
contains  a  central  plane  of  syrnnetry.  In  this  case  the  method  of  equivalent 
parameters  can  be  employed  (Herpin  1947,  Knittl  1976).  If  the  cell  is 
centrosymmetric,  the  parameter  p  defined  in  eq.  II-1-15  is  zero,  and 
eq.  II-1-14  becomes,  to  first  order  in  f  and  A 


( 
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PROPAGATION  EQUATION  FOR  AMPLITUDE 
REFLECTANCE 


•I  (r  -  Ip)  t*2V 


-I  (r  + 


The  amplitude  recursion  equation  is: 

PKO  *  e*2'*  PK  -  C»r  -p)  •*"  -  (Ir  «•  p)  e'3u  p* 


Figure  11-1-3 
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in  the  simplest  case  where  m»l.  ZC  this  aetrix  equation  is  ooRpared  to  the 
usual  (non  x-ray)  characteristic  matrix  solution  for  single  homogeneous  layers 
(Born  snd  Wolf,  1975,  p.58)  with  6*0: 


(l-l  -22) 


we  see  that  the  two  matricies  can  be  formally  equated  if  we  make  the 
identifications 


and  pe  are  known  as  the  equivalent  parameters  of  the  unit  cell. 
(Tl,  is  the  equivalent  index  and  the  equivalent  phase  thickness  of  the 
cell) .  There  is  s  straightforward  analogy  between  the  unit  cell  of  an  x-ray 
reflector  snd  s  single  homogeneous  layer  that  has  an  index  of  refraction 
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nc  and  a  phase  thickness  .  Similarly,  a  reflector  with  J  cells 
corresponds  to  a  stack  of  J  equivalent  layers,  which  is,  in  affect,  a  single 
layer  having  index  n#  and  phase  thickness  J  .  this  coalescence  of  the 
J  equivalent  layers  into  cne  is  permissible  because  the  internal  equivalent 
layers  have  no  interaction  with  cne  another  (since  their  refractive  indicies 
are  equal),  the  equivalent  parameters  can  therefore  be  considered  to 
incorporate  implicitly  the  effect  of  multiple  reflections  within  the 
structure. 

Bq.  II-1-20  is  also  analogous  to  the  propagation  equation  that  Kogelnik 
derives  (1976)  from  standard  coupled  wave  equations,  this  may  be  seen  by 
setting 


tK..  -X,  *  vf  •  AK  ■* 

■  +  °^) 


(1-1-24; 


so  as  to  obtain  the  differential  equation 


«  .  2ii?-(ir-p)-(ir4p)^40(A,;  +  0(9A)+0(^2; 

1  ^ 


(1-1-25; 


is  also  of 


where  we  have  indicated  that  the  neglected  tern  of  order 
order  ? 2  ,  f- A  ,  and  A*.  In  treating  stochastic  problems,  a  difference 
equation  has  the  important  advantage  over  differential  equations  of  having 
terms  whose  statistical  properties  can  be  determined  quite  easily:  this  is 
because  the  terms  of  the  difference  equation  have  a  cne-to-cne 
correspondence  with  the  properties  of  the  individual  cells. 

Since  the  layer  thicknesses  are  proportional  to  the  order  of 
diffraction  m,  the  absorption  is  least  in  first  order  (for  given  wavelength 
and  angle) ,  and  the  first  order  reflectivity  is  therefore  the  largest. 

Thus,  although  our  formalism  includes  all  orders  of  diffraction,  our 
discussion  will  generally  involve  only  the  first  order  case. 
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Section  13-2  -  Defect-Free  Multilayers 

Fart  A)  Reflectivity  of  Periodic  Multilayers 

The  simplest  ease  that  can  be  treated  with  the  formaline  developed  in 
*^e  previous  section  is  that  of  a  perfectly  periodic  multilayer  reflector. 
Periodic  x-ray  multilayers  have  been  treated  by  several  authors  in  recent 
years  (Vinogradov  and  Zeldovich  1977,  Lee  1981,  and  Underwood  and  Barbee 
1981).  Periodic  *-ray  multilayers  are  quite  analogous  in  their  diffracting 
properties  to  crystals  cut  in  the  Bragg  geometry. 

Me  will  present  brief  derivations  of  the  reflecting  properties  of 
periodic  *-ray  sultilayers;  the  derivations  will  tend  to  be  quite  sketrfiy 
except  where  previous  work  has  been  extended. 

As  disaissed  in  sec.  II-l-B,  periodic  multilayers  can  be  treated  with 
an  analysis  that  is  first  order  in  f  and  A  .  Eq.  II-1-20  becomes 


$*♦1  *  (‘r«-  P*)  ”  <irK+  P*)  Cl -2-1) 


The  reflectivity  of  the  periodic  multilayer  increases  almost 
monotonically  as  more  wit  cells  are  added,  until  a  steady-state  regime  is 
reached.  In  the  steady-state  regime,  we  can  set 


«... '  *  * 
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•o  that  ag.  11-2-1  becomes  a  quadratic  aquaticr  to  be  solved  for  The  two 
solutions  are 


C  -  (1-2-3) 

'  r-  i  p 


where  two  roots  result  frcn  the  sign  ambiguity  in  the  quantity 


5 


(1-2-4) 


We  new  discuss  the  sign  choice  that  must  be  made  here.  The  sign  choice 
that  %e  establish  will  be  applicable  in  most  of  the  remaining  calculations 
of  the  text.  The  few  exceptions  will  be  explicitly  identified. 

The  parameter  S  has  already  appeared  in  aq.  II-1-23;  in  eg.  11-2-4  we 
make  the  mnbiguity  in  sign  *jg>licit,  and  generalize  the  definition  to  include 
the  non- centre* ynrnetr ic  case. 

In  order  for  the  equivalent  phase  parameter  to  imply  exponential 
attenuation  instead  of  exponential  anplif ication,  we  must  choose  that  root 
in  eq.  II-1-23  which  causes  the  imaginary  part  of  S  to  be  positive.  In  the 
rsmainder  of  the  text,  the  square  roots  that  occur  in  expressions  of  the 
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form  /  a.1  -  ^ 1  will  by  definition  be  chosen  to  be  on  the  branch 
corresponding  to  positive  imaginary  parts. 

He  will  ncv  show  that  under  this  convention  it  is  the  upper  sign  rather 
than  the  lower  sign  that  must  be  used  in  eq  11-2-4  as  well  as  in 
eq.  11-1-23. 

If  we  multiply  the  two  possible  solutions  together  (i.e.  the  two 
solutions  corresponding  to  different  evaluations  of  the  sign  anbiguity  in 
eq.  11-2-4) ,  we  find 


c.c,  =  ( iilEEIEli  )  .  (  zlt1-  r»-  - 1  )  r+if 

m  'i  v  r-ip  /  V  r-ip  /  r-ip 
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bo  that  if  we  represent  the  true  physical  solution  with  an  unsubscripted  ^  : 

9l  x  (AL±)\  ( 13.Z±)  ( JLtil  .  rT\j>\  s  _  r*  ip 
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or 
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In  order  to  have  R  <  1,  we  require 


R -  -Tjp  Re(tj*)  <  0  U-*-» 


We  will  show  later  in  this  eecticn  that  at  the  wavelength  of  peak 
reflectivity, 


i 
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with  d>0  and  with  the  upper  and  lexer  sign  in  eq.  II-2-9  representing  the 
sane  sign  ambiguity  as  in  eq.  11-2-4. 

Thus,  from  eq.  II-1-15,  we  require  at  peak  reflectivity 


;  y*  <  o 
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for  R  ^1,  with  J4."  defined  to  be  the  imaginary  part  of  the  parameter  jj- 
(defined  in  eq.  11-1-15). 


Since  all  Materials  absorb  in  the  x-ray  regiae,  jl  >  0,  ao  that  in 
•gs.  II -2-4, 10  we  suat  choose  the  gpper  sign  when  the  Multilayer  operates  at 
the  Bragg  condition. 

Further,  we  trill  show  later  in  this  section  that  the  imaginary  part  of 
$  can  never  dtange  sign  from  its  value  at  the  Bragg  condition.  Therefore, 
we  will  always  choose  the  root  in  the  definit'on  of  S  to  be  the  root  with 
the  positive  imaginary  part. 

In  this  case  we  have 


S<i 

r-ip 


r+»p 
*- 1 


(l-2-il) 


and  also  the  useful  result 


p*  =  -  -J-VP  -i*t 

K  T-ip  X  -  t 


(a -2-12; 


(This  is  eq.  II-2-6  with  the  anbiguity  in  the  sign  of  t  new  resolved) . 
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In  Appendix  3  w  solve  aq.  11-2-1  in  the  ease  ef  periodic  multilayers 
that  do  not  contain  sufficiently  many  layer  pairs  to  reach  the  steady-state 
regime.  The  solution  (integration  of  eq.  11-1-25  with  constant  coefficients) 

is 


i  -  £  (■??£) 


(l  -  2-13) 


where  ^  is  the  steady-state  reflectivity  given  by  eq.  11-2-11.  When  the 
siultilayer  is  oentrosynjnetxic  (as  in  the  bilayer  ease  of  fig.  11-2-1) , 


.  _  ^(l-e2;Wy~n) 

i  _lii  1 1-1) 


(11-2-14} 


The  resartolanoe  of  eq.  11-2-14  to  an  Airy  sunmation  is  readily 
understood  in  terms  of  the  equivalent  index  analogy  (sec.  II-l-D) . 

The  most  important  oentrosyrrmetric  case  is  that  of  the  bilayer 
structure  shown  in  fig.  11-2-1.  This  figure  incorporates  the  ideal 
assimpticns  of  perfectly  sharp  interfaces.  Under  this  assumption  the 
structural  parameters  defined  in  eq.  II -1-15  beoane 

9  m  0 
rK 

r  «  (  A  -  A  )  Sift  3  JSC.2  0 

«  •«  l  r 

14,  m.  x.  A  +  fl  (  A  “  A  )  Mt  0 

f  K  t  H  L 

t  m  k  us  B  K  (J-2-1 5) 

r*,  *  •  * 
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The  parameter  pM  K  ,  the  phase  thickness  of  the  high  index  layer,  will 
be  used  ao  often  that  we  will  establish  a  convention  in  which  the  subscript  H  is 
made  implicit.  The  symbol  will  thus  by  definition  be  the  phase 
thickness  of  the  Kth  high  index  layer. 

The  high  index  layer  will  be  by  definition  that  layer  with  the  larger 
absorption;  generally  the  real  part  of  the  decrement  of  the  high  index  layer 
will  also  be  larger  in  magnitude  than  the  real  decrement  of  the  1cm  index 
layer. 


Jt 
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Part  B)  Structural  Optimization  of  Periodic  Multilayers 
We  now  turn  fran  the  calculation  off  reflectivity  to  the  associated 
optimization  problem,  i.e.  the  problem  off  determining  the  particular 
structure  that  optimizes  the  performance  of  a  multilayer.  Here  we  will 
mostly  be  concerned  with  maximizing  the  peak  reflectivity  of  the  structure. 
This  would  be  particularly  appropriate  in  such  applications  as  the  design  of 
cavity  mirrors  for  potential  x-ray  lasers;  in  other  applications  one  might 
ultimately  prefer  more  complicated  criteria. 

The  reflectivity  of  a  multilayer  mirror  operating  near  resonance  can 
for  practical  purposes  always  be  considered  to  be  an  increasing  function  of 
the  nisnber  off  layer  pairs  in  the  structure;  thus  in  order  to  maximize  the 
reflectivity  one  would  want  to  deposit  sufficiently  many  layer  pairs  that 
the  multilayer  operates  in  the  steady-state  regime.  (In  sane  cases  the 
reflectivity  is  found  to  undergo  very  small  oscillations  about  the 
steady-state  reflectivity  as  the  nvxnber  of  layer  pairs  approaches  the  regime 
where  radiation  no  longer  reaches  the  substrate) . 

Further,  in  many  applications  polarization  effects  can  be  neglected. 
This  includes  applications  at  normal  incidence  (x-ray  laser  mirrors, 
normal- incidenoe  microscopy) ,  grazing  incidence  applications  (all 
high-energy  applications) ,  reflections  taking  place  near  45  degrees  (where 
the  P  reflectivity  is  zero) ,  and  most  of  the  applications  that  involve 
synchrotron  radiation  (which  is  linearly  polarized) . 
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In  such  cases  w  Maximize  the  reflectivity  by  Maximizing  the  magnitude 
squared  of  eq.  21-2-11. 

We  first  consider  the  preplan  of  choosing  the  optimum  unit  call 
thickness  far  the  nultilayer,  which  from  eq.  Xl-l-U  is  equivalent  to 
finding  an  optimum  value  for  the  paraneter  f  . 

The  well-known  requirement  for  Maximizing  the  reflectivity  of  a 
dielectric  multilayer  (i.e.  a  multilayer  with  real  A(z))  ,  is  that  the  unit 
cell  thickness  be  chosen  in  accordance  with  Bragg's  law,  as  corrected  for 
dispersion.  By  this  criterion  the  optical  thickness  of  the  unit  cell  should 
be  one  half-wavelength,  so  that 

I  it  Re  (1  +AW)  LOS  e  (11-2-16) 

Period 

Here  0  (z)  is  the  angle  of  refraction  at  a  particular  depth  in  the 
structure;  this  is  defined  by 


(  1  ♦  4U))  Sift  6  (*)  *  sin  9  (1-2-17) 
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Then 
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or  from  eq.  II -1-11; 
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This  x-ray  version  of  the  dispersion-corrected  Bragg’s  lav  does  not 
oorrectly  produce  the  position  of  maximum  reflectivity  in  the  case  where  A 
is  aenplex.  The  appropriate  condition  in  this  case  was  first  found  in  the 


'  -«* 
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context  o f  diffraction  from  oentrosyranetric  crystals  by  Millar  (1935) . 

Me  will  briefly  present  a  stabler  derivation  that  applies  in  the 
non-oentr  osyimetr ic  case  as  veil,  and  that  provides  physical  insight  into 
the  origin  of  the  absorption  correction  to  Bragg's  law. 

Since  ,  we  can  say  that  if 


2.Re(iii]  s  ±  il!  +  J_  it 

z  Ke  U  f  +  ^  4v 


1  il 

ft  ’  Jv 
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(1-2-21) 


where  v  is  sane  structural  paramet  r  of  the  multilayer,  then  the 
reflectivity  is  maximized  with  respect  to  v.  In  this  section  we  will  use  a 
dot  to  represent  d  /dv. 

In  the  present  case  we  optimize  with  respect  to  the  parameter  f  ,  so 
that  t  »  1  to  first  order  in  f  according  to  eg.  II-1-15.  (As  discussed 
above,  in  the  case  of  periodic  multilayers  we  need  only  wort  to  within  first 
order  in  the  par meter  9  ,  and  can  neglect  terms  of  order  ?  ’  A . ) 

Also, 


Y  «p«/t«0  +  0(9’  A)  fl-2-22) 
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Since 


Re(f)  « 
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our  condition  for  maximizing  R  is  that 
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(1-2-25) 
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If  ve  cxnpare  eq.  II-1-23  with  aq.  21-2-26,  we  see  that  our  correction 
to  Bragg’s  law  is  a  requirement  that  the  real  part  of  the  equivalent  phase 
thickness  of  the  cell  be  7L .  The  difference  here  from  the  usual  requirement 
that  the  real  part  of  the  optical  phase  thickness  be  71  is  due  to  the 
presence  of  multiple  reflections. 

This  may  be  seen  by  manipulating  eq.  11-2-11  to  obtain 

5  -  -*-*  (r-ip) 

so  that 

s  * 

S 

to  within  order  A  . 

The  term  -  ft'1*  represents  the  phase  oscillation  imrelved  in  transmission 
through  the  cell.  The  term  ♦  p  ( l f  +  p)  also  contributes  to  the  overall 

phase  oscillation  across  the  cell  (which  is  )»  this  term  represents  a  multiple 

reflection  process  consisting  of  transmission  through  the  cell,  reflection 
from  all  succeeding  cells,  and  reflection  back  into  the  initial  direction. 
(Higher-order  sultiple  reflections  are  of  order  A*).  A  phase  change  occurs 
during  the  multiple  reflections,  and  if  n  is  complex,  the  overall 


-e*'4  (1-  ip  (r- ip;) 
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oscillation  of  the  field  across  the  aell  is  shifted. 

For  brevity's  sake  we  have  described  the  calculation  of  this  affect  as 
a  correction  for  absorption;  hmever  we  note  that  the  imaginary  part  of  A 
causes  a  shift  in  the  position  of  peak  reflectivity  through  its  affect  on 
the  phase  of  ^  ,  rather  than  by  any  extinction  mechanism  (see  below) . 

We  can  convert  the  resonanoe  condition  of  eg.  II-2-26  to  a 
computationally  more  useful  form  as  follows.  The  requirement  that  Re($ )*0 
can  only  be  satisfied  if  is  a  negative  real  number.  He  show  in 
Appendix  4  that  when  f  is  chosen  to  make  4*  real,  4*  is  automatically 
made  negative.  This  implies  a  result  stated  above,  that  the  imaginary  part 
of  S  can  never  change  sign  from  its  sign  at  the  Bragg  condition.  Such  a 
sign  change  would  require  the  existence  of  a  f  value  at  which  w  was  both 
real  and  positive. 

For  Im(  S1)  to  be  zero,  we  must  have  from  eq.  1 1-2-4 
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To  first  order  in  L  ,  we  csn  evaluate  eq.  11-1-15  (the  defining  equations 
for  the  parameters  in  our  solution  for  for*  )  at  Btm  . 

In  terms  of  a  twvelength  or  ^spacing  shift  away  fran  the  uneorrected 
Bragg  value  defined  by  2d  COS  0#*  X  ,  we  can  set 

*  -  *  *  -  Jl  ~~  (1-2'Sl) 


To  convert  the  phase  shift  to  an  angular  shift,  we  have  from 


eq.  11-1-11 
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where 
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We  can  set 
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Terms  of  order  A0Xcos  0#  have  been  included  in  eq.  11-2-34  because  the 
linear  term  becomes  very  anall  near  normal  incidence. 

Solving  for  A0  gives 


a6 


+ 


2  y-t 

TL 


txn.  Q 
0 
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which  has  physical  solutions  so  long  as  the  radicand  is  non-negative,  or,  to 
a  good  approximation,  so  long  as  0^ ^ .  At  B^m ,  the  true 
anqle  of  maximum  reflectivity  6  is  at  0*.  _ 
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(1-2-36) 


In  eq.  11-2-36  the  first  term  on  the  right  is  the  dispersion  correction  and 
the  second  is  the  absorption  correction. 

He  show  in  Appendix  5  that  the  absorption  correction  must  always  be 
smaller  than  the  dispersion  correction  except  in  regions  of  very  strong 
anomalous  dispersion,  but  there  are  many  cases  in  which  the  two  are 
ccnparable  in  magnitude. 
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In  fig.  II-2-2  we  shcx  plotted  the  ratio  of  absorption  correction  to 
net  correction  in  a  particular  exanple;  specifically,  the  ratio  plotted  is 


(1-2-37) 


with  6X  the  true  Bragg  angle  as  corrected  for  both  absorption  and 
dispersion,  and  Bt  the  Bragg  angle  as  corrected  by  eq.  11-2-20  for 
dispersion  only.  This  example  exhibits  the  cornonest  ease  in  which  the 
absorption  and  dispersion  corrections  are  in  opposite  directions.  (Note 
that  in  this  figure  the  incidence  angle  is  measured  relative  to  the 
surface.) 

Fig.  II-2-2  shows  three  regimes  in  which  the  relative  effect  of 
absorption  is  fairly  large. 

First,  near  normal  incidence,  eq.  II-l-ll  shows  that  the  correspondence 
between  angular  shifts  and  shifts  in  phase  thickness  be ocrnes  strongly 
non-linear  due  to  the  cosine  factor.  Jn  the  normal  incidence  regime,  a 
relatively  small  absorption  correction  to  the  phase  thickness  can  therefore 
cause  a  large  decrease  in  the  total  angular  shift.  Both  the  angular 
absorption  correction  and  the  angular  dispersion  correction  become  large  in 
absolute  terms  as  the  angle  of  incidence  approaches  0*. 

The  plot  shows  two  other  regions  at  which  the  relative  importance  of 
the  absorption  correction  becomes  strong;  at  wavelengths  just  above  the 
car boo- K  edge,  and  at  short  wavelengths.  In  each  case  the  absorption 
correction  becomes  strong  because  multiple  reflections  beecme  strong,  and 
multiple  reflections  are  strong  precisely  where  the  absorption  imposes  a 
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COMPARISON  OF  ABSORPTION  CORRECTION 
TO  DISPERSION  CORRECTION 


ANCLE  OF  PEAK  REFLECTIVITY 
(to  turtocc) 


Figure  1 1-2-2 
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weaker  limit  cr  the  reflectivity,  (in  .  *onent  we  will  coraider  the 
limiting  case  where  the  absorption  goes  to  uro) . 

This  is  why  the  term  "absorption  correction*  must  not  be  taken  too 
literally;  as  discussed  above,  the  imaginary  part  of  A  shifts  the  resonance 
angle  through  the  mechanism  of  phase  changes  that  occur  during  multiple 
reflections,  rather  than  ty  any  mechanism  Inwlving  extinction  of  the 
fields. 

In  the  limit  of  sero  absorption,  the  steady-atate  reflectivity  of  the 
structure  approaches  unity  throughout  the  high  reflectance  sane  (often 
called  the  stqp-hand) .  According  to  the  above  analysis,  in  the  limit  of  low 
absorption  the  reflectivity  approaches  closest  to  one  at  an  angle  quite  far 
front  the  aenter  of  the  stop-band,  because  with  low  absorption  multiple 
reflections  became  quite  strong,  which  in  turn  causes  the  ratio  defined  in 
eq.  II-2-37  to  become  quite  large.  However  with  very  low  absorption  the 
reflectivity  at  0,  departs  only  slightly  further  from  unity  than  the 
reflectivity  at  ©x,  and  indeed  the  reflectivity  throughout  the  stopband  is 
almost  constant. 

In  contrast,  in  the  soft  x-rsy  regime  where  absorption  is  fairly 
substantial,  peak  reflectivities  generally  do  not  approach  very  close  to  one 
and  the  absorption  correction  manifests  itself  as  a  distinct  shift  in  the 
position  of  peak  reflectivity  (aee  fig.  II-2-3) . 
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REFLECTANCE  OF  W/C  MULTILAYER  - 
UNPOLARIZED  LIGHT 
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Figure  II- 
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Once  the  cell  thickness  of  the  structure  is  chosen  so  as  to  optimize 
the  reflectivity,  one  may  wish  to  optimize  the  structure  in  other  degrees  of 
freedom;  for  exanple  in  the  usual  case  of  a  bilayer  reflector  erte  would  like 
to  optimize  the  ratio  in  thickness  of  the  two  types  of  layers. 

In  order  to  maximize  the  reflectivity  at  the  Bragg  angle,  (which  we 
will  refer  to  as  maximization  of  peak  reflectivity) ,  we  use  the  steady-state 
version  of  eq.  II -2-1 


-Cir-p)  -  (ir  +  f»)  *  0 


(1-2-38) 


Dividing  by  ^  and  differentiating  with  respect  to  sane  secondary  structural 
parameter  (any  parameter  other  than  ?  )  we  obtain 
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+  (r-ip)^-  (r  +  ip)  X  =  0  (1-2-39) 
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Fran  eq.  IT-2-11  we  have 


•  •  * 

(T-  i  p)^-  (r+  ip)  *  -i-  +  =  -ZSJ  U-2-40) 

To  maximize  the  peak  reflectivity  we  require  that  Re ( — -  )*0;  further, 

* 

since  we  are  operating  at  peak  reflectivity,  $  is  pure  imaginary,  so  that 
lm(2  S  4-1-0. 

Thus  if  we  take  the  imaginary  part  of  eq.  JI-2-39  we  obtain,  using 
eq.  XX -2-11,  the  condition  for  maximizing  the  peak  reflectivity 

1,1  [“"?=£  <*♦*>♦  Inf  »-«]  *  o 

(1-2-41) 

If,  as  in  the  case  of  a  bilayer  reflector,  the  cell  is  oentro6ynmetxic 
in  such  a  way  that  p  is  zero  along  with  p,  eq.  II-2-41  achieves  the  sinple 
form: 


I  m 


■  0 


(1-2-42) 


XI-2-24 


He  now  aonsider  optimization  with  respect  to  the  parameter  p  used  in 
characterizing  a  bilayer  reflector  (eq.  11-2-15).  The  terms  of  eq.  XX-2-41 
are  of  order  A  ,  so  that  to  obtain  optimisation  accurate  overall  to  first 
order  in  A  ,  we  need  only  differentiate  with  respect  to  ^  in  "zeroth"  order 
for  our  optimization. 

For  the  same  reason,  we  can  aonsider  optimization  of  ^  to  yield  an 
optimum  ratio  of  high-index  layer  thickness  to  low-index  layer  thickness 
given  fay 


g  k  cos  6  dM  _  dH 
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Upon  carrying  out  the  differentiation  in  eq.  XI-2-42  we  obtain  the 
optimization  condition  first  obtained  by  Vinogradov  and  Zeldovich  (1977) 
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jBq.  u-2-44  is  valid  for  either  polarization  at  any  angle  of  incidence. 


We  Rentier)  that  one  should  avoid  being  mislead  by  the  dissuasion 
following  sq.  23b  in  Vinogradov  and  Zeldovich  (1977) ,  which  sight  suggest 
that  eg.  II-2-44  can  be  treated  as  a  primary  optimization  condition. 

Bq.  u-2-44  is  a  secondary  optimization  condition  that  obtains  only  if  the 
total  cell  thickness  has  also  been  optimized. 

Bq.  11-2-44  is  a  tranoendental  equation  with  no  analytic  solution,  but 
in  Appendix  6  we  exhibit  a  fast  numerical  prooedure  that  can  be  used  to 
solve  it.  the  initial  analytic  seed  that  is  used  in  this  prooedure  may  well 
be  accurate  enough  for  most  purposes. 

In  Appendix  7  we  show  that  the  peak  reflectivity  (i.e.  the  releetivity 
at  the  Bragg  condition)  is  given  by: 
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Part  C)  Optimum  Materials  Choices 

Bq.  11-2-45  in  conjunction  with  the  fast  algorithm  of  Appendix  6  makes 
it  possible  to  very  rapidly  calculate  the  maximan  reflectivity  Obtainable 
from  a  given  pair  of  materials  reflecting  at  a  particular  x-ray  wavelength. 

It  therefore  becomes  possible  to  have  a  ocrputer  program  search 
efficiently  through  all  possible  pairs  of  sultilayer  materials  and  so 
determine  the  optimum  reflectivity  attainable  at  the  particular  wavelength 
under  consideration. 

Fig.  11-2-4  and  table  XI-2-1  show  the  results  of  such  a  materials 
search  at  each  of  125  wavelengths  in  the  soft  x-ray  region  from  6.2A  to 

124A.  A  and  a  were  calculated  using  atomic  scattering  factors  cnnpiled 

H  i 

by  Henke  et  al.  (1982). 

The  calculated  reflectivities  of  the  optimum  materials  pairs  are  shown 
plotted  in  red.  For  comparison,  fig.  11-2-5  shows  the  reflectivity 
obtainable  from  the  usual  materials  choice  of  tungsten  and  carbon.  The 
reflectivity  of  the  new  materials  combinations  are  surprisingly  high, 
particularly  at  the  longer  wavelengths.  In  general  the  high  reflectivities 
are  a  consequence  of  atomic  resonances;  the  anomalous  dispersion  in  the 
spectral  vicinities  of  such  resonances  cause  layers  made  with  the  elements 
to  yield  high  reflectivities,  primarily  because  their  unit  decrements  have 
anomalously  small  imaginary  parts. 

Fig.  II-2-4  also  shows  plotted  in  green  the  optimum  fraction  of  the 
unit  cell  thickness  that  should  be  occupied  by  the  high  index  layer. 


t 


XX-2-28 


Optimum  X-Ray  Multilayer  Parameters 

MINIMUM  l  ALLOWED  IS  5 

WAVE- 


LENGTH 

R 

N 

GAMMA 

ELEMENT  -  8 

ELEMENT  -  L 

( 

AT  ) 

(  DH  ) 

( 

NORMAL  ) 

( - ) 

(INCIDENCE) 

(DH  *  DL) 

124.00 

.8447 

9.9 

.3568 

RHODIUM 

SILICON 

121.04 

.7932 

15.4 

.3627 

RHODIUM 

STRONTIUM 

118.15 

.7992 

15.8 

.3785 

RHODIUM 

STRONTIUM 

115.33 

.7897 

16.8 

.3850 

RHODIUM 

STRONTIUM 

112.58 

.7758 

18.0 

.3862 

RHODIUM 

STRONTIUM 

109.89 

.7667 

19.1 

.3869 

RHODIUM 

STRONTIUM 

107.27 

.7505 

19.9 

.3785 

RUTHENIUM 

STRONTIUM 

104.71 

.7392 

20.9 

.3800 

RUTHENIUM 

STRONTIUM 

102.21 

.7280 

21.9 

.3814 

RUTHENIUM 

STRONTIUM 

99.77 

.7171 

23.0 

.3829 

RUTHENIUM 

STRONTIUM 

97.39 

.7165 

23.6 

.3899 

SILVER 

STRONTIUM 

95.06 

.7176 

24.8 

.4049 

SILVER 

STRONTIUM 

92.79 

.6966 

26.4 

.4099 

SILVER 

STRONTIUM 

90.58 

.6693 

27.7 

.4205 

SILVER 

STRONTIUM 

88.42 

.6313 

29.1 

.4344 

SILVER 

STRONTIUM 

86.31 

.6218 

42.6 

.3573 

RUTHENIUM 

BORON 

84.25 

.6210 

45.4 

.3503 

RUTHENIUM 

BORON 

82.23 

.6201 

48.2 

.3456 

RUTHENIUM 

BORON 

80.27 

.6208 

51.3 

.3388 

RUTHENIUM 

BORON 

78.36 

.6239 

54.5 

.3320 

RUTHENIUM 

BORON 

76.49 

.6287 

57,7 

.3261 

RUTHENIUM 

BORON 

74.66 

.6341 

60.9 

.3204 

RUTHENIUM 

BORON 

72.88 

.6414 

63.8 

.3146 

RUTHENIUM 

BORON 

71.14 

.6504 

66.2 

.3100 

RUTHENIUM 

BORON 

69.44 

.6774 

79.9 

.3797 

LANTHUNUM 

BORON 

67.78 

.7211 

75.3 

.3905 

LANTHUNUM 

BORON 

66.16 

.9200 

52.4 

.4033 

LANTHUNUM 

BORON 

64.59 

.5446 

99.7 

.4933 

CALCIUM 

LANTHUNUM 

63.04 

.5097 

108.6 

.4924 

LANTHUNUM 

CALCIUM 

61.54 

.4932 

73.8 

.3703 

RUTHENIUM 

CALCIUM 

60.07 

.4879 

78.6 

.3668 

RUTHENIUM 

CALCIUM 

58.64 

.4827 

93.6 

.3634 

RUTHENIUM 

CALCIUM 

57.24 

.4778 

99.0 

.3598 

RUTHENIUM 

CALCIUM 

55.87 

.4732 

94.8 

.3563 

RUTHENIUM 

CALCIUM 

54.54 

.4676 

101.4 

.3521 

RUTHENIUM 

CALCIUM 

53.23 

.4687 

161.0 

.3011 

RUTHENIUM 

CARBON 

51.96 

.4707 

172.1 

.2956 

RUTHENIUM 

CARBON 

50.72 

.4797 

155.5 

.2204 

COBALT 

CARBON 

49.51 

.4943 

176.9 

.2623 

CHROMIUM 

CARBON 

49.33 

.5109 

183.2 

.2610 

CHROMIUM 

CARBON 

47.18 

.5291 

188.5 

.2596 

CHROMIUM 

CARBON 

Table  11-2-1 
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MINIMUM  X  ALLOWED  IS  5 


HAVE- 


LENGTH 

R 

■ 

CAMKA 

ELEMENT  -  ■ 

ELEMENT  -  L 

(  AT  ) 

(  DB  ) 

(  NORMAL  ) 

(INCIDENCE) 

(DH  4  DL) 

46.05 

.5530 

191.5 

.2561 

CHROMIUM 

CARBON 

44.95 

.5830 

189.3 

.2568 

CHROMIUM 

CARBON 

43.06 

.6625 

158.3 

.2556 

CHROMIUM 

CARBON 

42.03 

.4987 

157.2 

.3122 

CHROMIUM 

CALCIUM 

41  •  01 

.5072 

163.9 

.3126 

CHROMIUM 

CALCIUM 

40.01 

.5159 

170.6 

.3129 

CHROMIUM 

CALCIUM 

39.04 

.5256 

177.5 

.3134 

CHROMIUM 

CALCIUM 

30.09 

.5358 

103.8 

.3137 

CHROMIUM 

CALCIUM 

37.96 

.5489 

108.9 

.3141 

CHROMIUM 

CALCIUM 

37.05 

.5663 

191.6 

.3145 

CHROMIUM 

CALCIUM 

36.17 

.5924 

108.9 

.3149 

CHROMIUM 

CALCIUM 

35.30 

.4174 

164.9 

.3257 

COBALT 

SCANDIUM 

34 .46 

.4302 

171.0 

.3260 

COBALT 

SCANDIUM 

33.64 

.4458 

176.3 

.3263 

COBALT 

SCANDIUM 

32.04 

.4645 

180.1 

.3267 

COBALT 

SCANDIUM 

32.05 

.4904 

100.2 

.3271 

COBALT 

SCANDIUM 

31.29 

.5392 

198.7 

.3790 

CHROMIUM 

SCANDIUM 

30.54 

.3595 

101.5 

.3653 

NICKEL 

BARIUM 

29.81 

.3708 

198.1 

.3596 

NICKEL 

TITANIUM 

29.10 

.3969 

200.5 

.3599 

NICKEL 

TITANIUM 

20.41 

.5022 

144.4 

.4967 

CHROMIUM 

CADMIUM 

27.73 

.5053 

184.3 

.3745 

COBALT 

TITANIUM 

27.07 

.3960 

274.1 

.3405 

NICKEL 

MAGNE5.UM 

26.42 

.4058 

291.0 

.3473 

NICKEL 

MAGNESIUM 

25.79 

.4153 

308.9 

.3462 

NICKEL 

MAGNESIUM 

25.17 

.4254 

328.1 

.3451 

NICKEL 

MAGNESIUM 

24.57 

.4750 

209.3 

.3992 

NICKEL 

VANADIUM 

23.99 

.4443 

370.2 

.3429 

NICKEL 

MAGNESIUM 

23.41 

.4537 

393.6 

.3418 

NICKEL 

MAGNESIUM 

22.05 

.4625 

418.3 

.3409 

NICKEL 

MAGNESIUM 

22.31 

.4970 

216.3 

.4319 

NICKEL 

TELLURIUM 

21.78 

.7359 

250.8 

.4003 

TELLURIUM 

MAGNESIUM 

21.26 

.4605 

504.9 

.3377 

NICKEL 

MAGNESIUM 

20.75 

.5824 

522.0 

.4019 

CHROMIUM 

MAGNESIUM 

20.25 

.5038 

575.0 

.3356 

NICKEL 

MAGNESIUM 

19.77 

.5109 

614.5 

.3346 

NICKEL 

MAGNESIUM 

19.30 

.5173 

657.9 

.3336 

NICKEL 

MAGNESIUM 

10.04 

.5236 

707.6 

.3322 

NICKEL 

MAGNESIUM 

10.39 

.5291 

761.6 

.3308 

NICKEL 

MAGNESIUM 

17.95 

.5321 

020.4 

.3302 

NICKEL 

MAGNESIUM 

17.52 

.6007 

800.3 

.3609 

IRON 

MAGNESIUM 

17.10 

.5338 

969.6 

.3286 

NICKEL 

MAGNESIUM 

16.69 

.5320 

1061.5 

.3276 

NICKEL 

MAGNESIUM 

16.30 

.5295 

1126.9 

.3140 

COPPER 

MAGNESIUM 

15.91 

.6136 

1244.9 

.3971 

BARIUM 

MAGNESIUM 

15.53 

.5446 

906.2 

.1821 

RHENIUM 

MAGNESIUM 

15.16 

.5549 

1046.9 

.1007 

RHENIUM 

MAGNESIUM 

14.79 

.5648 

1111.2 

.1795 

RHENIUM 

MAGNESIUM 

14.44 

.5743 

1179.1 

.1782 

RHENIUM 

MAGNESIUM 

14.10 

.5039 

1251.4 

.1770 

RHENIUM 

MAGNESIUM 

13.76 

.5935 

1320.6 

.1757 

RHENIUM 

MAGNESIUM 

13.43 

.6027 

1410.4 

.1745 

RHENIUM 

MAGNESIUM 

13.11 

.6117 

1496.3 

.1734 

RHENIUM 

MAGNESIUM 

12.00 

.6200 

1507.6 

.1722 

RHENIUM 

MAGNESIUM 

Table  11-2-1 
(continued ) 
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MINIMUM  Z  ALLOWED  IS  5 


WAVE- 


LENGTH 

R 

M 

GAMMA 

ELB4ENT  -  B 

ELEMENT  -  L 

(  AT  ) 

(  DH  ) 

(  NORMAL  ) 

( - ) 

(INCIDENCE) 

(DB  4  DL) 

12.49 

.6298 

1685.1 

.1711 

RHENIUM 

MAGNESIUM 

12.19 

.6386 

1787.5 

.1700 

RHENIUM 

MAGNESIUM 

11.90 

.6475 

1895.0 

.1690 

RHENIUM 

MAGNESIUM 

11.62 

.6565 

2010.0 

.1680 

RHENIUM 

MAGNESIUM 

11.34 

.6650 

2132.0 

.1671 

RHENIUM 

MAGNESIUM 

11.07 

.6731 

2260.9 

.1661 

RHENIUM 

MAGNESIUM 

10.81 

.6816 

2398.4 

.1651 

RHENIUM 

MAGNESIUM 

10. 55 

.6896 

2543.6 

.1642 

RHENIUM 

MAGNESIUM 

10.30 

.6978 

2695.9 

.1632 

RHENIUM 

MAGNESIUM 

10. OS 

.7059 

2854.0 

.1623 

RHENIUM 

MAGNESIUM 

9.81 

.7148 

3008.6 

.1614 

RHENIUM 

MAGNESIUM 

9.5B 

.7652 

5512.9 

.4435 

BORON 

MAGNESIUM 

9.35 

.6502 

2747.2 

.1900 

RHENIUM 

BORON 

9.13 

.6581 

2932.2 

.1688 

RHENIUM 

BORON 

8.91 

.6659 

3134.8 

.1875 

RHENIUM 

BORON 

8.70 

.6733 

3352.0 

.1863 

RHENIUM 

BORON 

8.49 

.6802 

3584.7 

.1852 

RHENIUM 

BORON 

8.29 

.6869 

3840.5 

.1840 

RHENIUM 

BORON 

8.09 

.6941 

3472.8 

.1934 

RHENIUM 

ALUMINUM 

7.89 

.6991 

4427.4 

.1815 

RHENIUM 

BORON 

7.71 

.7077 

4862.6 

.2011 

RUTHENIUM 

BORON 

7.52 

.7166 

5156.5 

.2003 

RUTHENIUM 

BORON 

7.34 

.7252 

5469.1 

.1996 

RUTHENIUM 

BORON 

7.17 

.7338 

5802.7 

.1988 

RUTHENIUM 

BORON 

7.00 

.7420 

6157.1 

.1961 

RUTHENIUM 

BORON 

6.83 

.7500 

6535.3 

.1973 

RUTHENIUM 

BORON 

6.67 

.7578 

6941.6 

.1965 

RUTHENIUM 

BORON 

6.51 

.7652 

7374.8 

.1958 

RUTHENIUM 

BORON 

6.35 

.7724 

7836.9 

.1950 

RUTHENIUM 

BORON 

6.20 

.7795 

8327.8 

.1943 

RUTHENIUM 

BORON 

Table  11-2-1 
(Continued) 


REFLECTIVITY  ATTAINABLE  FROM  TUNGSTEN- 
CARBON  MULTILAYERS 


WAVELENGTH  (A) 


Adjacent  pairs  of  dotted  lines  in  the  figure  dmaroate  the  regions  in 
which  some  one  pair  of  aaterials  has  proved  to  be  optimum;  atonic  resonances 
cause  the  curves  to  be  highly  discontinuous  across  these  lines.  Fig.  II-2-4 
shore  plotted  in  blue  the  bench-mark  value  N  that  the  mmber  of  layer  pairs 
J  should  exceed  in  order  to  approach  the  steady-state  regime.  N  is 
calculated  at  normal  incidence;  the  nanber  of  layer  pairs  required  scales  as 
cos1  6  for  fixed  X  (see  sec.  11-3) .  in  order  to  have  the  layers  remain 
of  practical  thickness  (say  at  least  a  few  Angstroms) ,  it  is  desirable  to  go 
to  reflection  angles  away  from  normal  incidence  if  the  wavelength  falls 
below  about  30 A.  Even  for  wavelengths  as  long  as  50 X,  we  would  expect  the 
reflectivities  to  be  degraded  significantly  from  the  displayed  values  at 
angles  near  normal  incidence  (see  chapter  III) . 

There  are  of  course  many  other  qualifications  to  be  made  about  the 
results  of  this  optimization  study,  ihe  optical  constants  are  presumably 
less  reliable  in  the  regions  of  strong  anomalous  dispersion  tnat  acranonly 
occur  in  the  elemental  selections  of  table  1 1 -2-1;  this  is  particularly  true 
in  the  case  of  elements  that  are  somewhat  esoteric.  In  fact,  optimization 
searches  such  as  this  will  tend  automatically  to  have  a  statistical  bias  in 
favor  of  erroneous  data  points  in  which  the  error  happens  to  result  in  a 
higher  calculated  performance. 

Further,  no  effort  was  made  to  assess  the  individual  suitability  or 
joint  conpatibility  of  the  ctosen  materials  in  terms  of  either  layer  or 
interfacial  quality  vd  stability;  even  freedom  f ran  chemical  reactivity  was 
not  considered.  The  only  requirements  on  the  elements  examined  were  that 
they  be  nan-radioactive  and  have  melting  points  greater  than  270* C  (that  of 
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bismuth) .  In  Addition  beryllim  mbs  excluded  for  its  toxicity;  unlike  ether 
toxic  elements,  beryllian  Mould  have  been  frequently  selected  by  the  Beard; 
program  had  it  been  included,  due  to  its  low  atonic  mrfcer. 

However,  Host  of  the  materials  selections  in  table  11-2-1  seen  at  least 
to  be  fairly  reasonable,  with  the  exception  of  two  fairly  oaron  selections 
for  the  lew  index  layers,  calcium  and  strontion,  which  Bight  be  difficult  to 
fabricate  in  thin  high  quality  layers  of  the  kind  required  in  x-ray 
multilayers  (Spiller  1982a) . 

In  sunmary,  the  results  of  this  materials  search  should  be  treated  as 
suggestive  until  the  materials  selections  can  actually  be  tried.  For  this 
reason  most  of  the  text  will  deal  with  the  more  carman  materials  choice  of 
tungsten  and  carbon. 

We  have  also  carried  out  a  multilayer  materials  search  using  a  modified 
optimization  program  that  seeks  to  maximize  integrated  reflectivity  (or 
collection  solid  angle) ,  rather  than  peak  reflectivity.  In  the  revised 
program  we  have  modified  the  output  routines  in  order  to  have  printed  out  a 
number  of  possible  materials  pairs  for  each  wavelength. 

An  abbreviated  tabulation  of  these  new  results  is  given  in  Appendix  16. 


Section  11-3  -  Sealing  of  Miltilayer  Reflection  Properties 

In  this  section  we  consider  the  approximate  angular  and  spectral 
scaling  of  the  reflection  properties  of  tungsten- car  bon  multilayers.  Our 
intention  will  be  to  determine  simple  dependencies  that  will  be  accurate  to 
within  a  factor  of  two  or  so  throughout  the  soft  x-ray  region. 

Fig.  II-2-5  shows  the  maximum  reflectivity  attainable  fran  tungsten  and 
carbon  as  a  function  of  wavelength.  The  results  are  for  S  polarization;  in 
the  case  of  P  polarization  the  reflectivity  of  an  individual  unit  cell,  and 
so,  approximately,  the  reflectivity  of  the  entire  structure,  is  reduced  by  a 
factor  cos2  2  0  (see  eqs.  II-1-15  and  11-3-1). 

Fig.  11-3-1  shows  the  values  used  to  obtain  the  reflectivities  in 
fig.  11-2-5.  (Fig  II-3-1  is  essentially  a  plot  of  the  solution  to 
eg.  11-2-44  based  on  the  data  in  (Henke,  et  al.  1982)).  The  optimum  thicknesses 
of  the  high  index  tungsten  layers  become  small  in  regions  of  low  absorption, 
i.e.  p>  is  snail  at  short  wavelengths  and  at  wavelengths  just  above  the 
carbon-K  edge.  In  fact,  Spiller  (1976)  has  shewn  that,  in  principle,  unit 
reflectivities  can  be  obtained  in  the  limit  of  no  absorption  in  the  low 
index  layers  (but  with  finite  absorption  in  the  high  index  material)  if  the 
high  index  layers  are  made  infinitely  thin. 


SCALING  OF  OPTIMUM  THICKNESS 
RATIO  PARAMETER 
(Tungsten-Carbon  Multilayer) 


WAVELENGTH  (A) 


X346 


Figure  11-3-1 
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The  reflectivity  of  the  tungsten-carbon  combination  is  usually  fairly 
■mall.  When  the  reflectivity  is  wall,  the  fourth  term  on  the  right-hand 
side  of  eq.  II-2-1  is  wall  because  it  represents  the  effect  of  Multiple 
reflections.  When  eg.  II-2-1  is  linearized  by  neglecting  this  term,  it  has 
the  steady-state  solution  ^  S£  -  ,  or 


R  « 


M.’2) 


(1 -3-1) 


This  Lorentzian  intensity  profile  was  first  found  by  Henke  (1982a).  If 
we  neglect  the  2  term  in  the  denominator  of  our  solution  for  multilayer 
reflectivity  outside  the  steady-state  regime  (eq.  II-2-14) ,  we  obtain 
another  result  of  Henke's  (1982b): 


j+ i 


(  1  -  e 


■US  (J-l)y 


(I -3-2) 


From  these  approximate  expressions  we  can  derive  simplified  formulas 
for  multilayer  properties,  whose  wavelength  and  angular  scaling  is  then 
relatively  easy  to  determine. 

For  example  f ran  eq.  II-3-1  we  see  that  the  fWHM  of  the  reflection 
profile  in  radian  vrits  is  approximately 


S<f  s  2/  (1-3-3) 

FWHH 


Also,  according  to  eq.  II-3-2  the  nunber  of  layer  pairs  N  required  for 
eq.  11-2-14  to  approach  the  Darwin-Prins  limit  of  eq.  11-2-11  is 
approximately ($*) *.  Further,  since  eq.  11-3-2  applies  under  the  condition 
that  ^is  snail,  frcn  eq.  11-3-1  we  see  that  the  parameter  r  can  be 
considered  to  be  snail  compared  to  the  parameter  t  ,  to  the  extent  that 
eq.  II-3-1  is  applicable.  In  that  case,  the  nunber  of  layer  pairs  N 
required  in  the  multilayers  is  given  approximately  by 

A/  ■  [  I  m  J  (y-M)1-  V2  ] 

£  1/ajl'  (*-3-0 

We  can  determine  the  approximate  wavelength  scaling  of  the  optical 
constants  of  tungsten  and  carbon  from  the  sending  plots  of  figs.  II-3-2, 
which  are  based  on  a  preliminary  version  of  the  data  in  (Henke,  et 
al.  1982). 

We  find 

(6A<\<  114a)  A*  *  ([ 1 7 < to* } A^J )  +  i (U 6 *  10* ‘] A^5} ) 

(6A<K<44A)  Ac  •  i(UD'10*]K*.)) 

(44A<\<il4A)  A.  *  CC 6.1  *ldT]  A**?)  *  *  il4.0  *  10*IC]  ) 
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WAVELENGTH  SCALING  OF 
UNIT  DECREMENT  -  TUNGSTEN  AND  CARBON 


WAVELENGTH  (A)  WAVELENGTH  (A) 


pc  =  2.0  gm/em3 
p  =15.0gm/cm3 

Results  are  from  a  preliminary  version  of  a  compilafion 
by  Henke  (1982a). 


Figure  1 1-3-2 


i  ^ f 


n-a-€ 


The  first  quantity  dependent  on  these  indicies  that  we  will  consider  is 
the  optimum  thickness  ratio  p .  Sinoe  lm(  Ac  )  tends  to  be  ■nail  compared 
to  Im(  A  )/  we  can  consider  6.  .  to  be  less  than  or  of  order  unity, 

W  * 

and  to  a  rough  approximation,  can  set 


(1-3-6; 


Eq.  11-2*44  then  be  cones  for  a  tungsten-carbon  multilayer: 


JCA„ 


A*  -  A* 

•  t 


(1-3-7) 


Then 


(i  -a-#; 


and  for  X  <  44A,  we  have  from  eq.  11-3-5 


1.86  «  fO’JXa  i 

l-J.iO’V' 


(1-3-7) 


n-3-7 


n-3-e 


Thus,  our  overall  result  for  the  range  6A  <  X  <  114A  is 


r  o.i5 

•.7 

A  * 

.  0.04 


(\<  44  A)* 

(K>44A)  . 


(1-3-13) 


This  power-law  fit  is  plotted  in  fig.  IX-3-1.  As  the  figure  shows,  the 
dependence  of  Po,i  0,1  X  is  not  very  strong.  Por  many  purposes,  it  is 
sufficiently  accurate  to  take  6  «  1. 

I  Opt 
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With  this  approximation,  we  nw  consider  the  approximate  scaling  of  the 
multilayer  resolution  and  required  number  of  layer  pairs  N.  We  first 
consider  the  dependence  of  these  two  quantities  on  d- spacing. 

If  we  take  ^  ■  1,  we  have  from  eq.  II-2-15 

aS  -  [  71  i  •  C  A*  -  L*  )]  •  ddc1© 

01-3-14} 

s  dec* 0  *  6.6  a  10  *  X *',f  Sec* 8 

w 


Taking  \ 


*  30A,  we  can,  to  a  fairly  good  approximation,  set 


ju"  s  6.6  *  to *  A*  Cx‘#  ,r)  dec1© 

-  4.0  *  10~‘ dec1© 


(1-3-15) 


We  note  that  eq.  II-3-15  essentially  states  that  the  spatially  averaged 
imaginary  part  of  the  index  of  refraction  scales  as  A  .  Thus,  the 
well-known  scaling  law  for  the  hard  x-ray  regime,  which  states  that  the 
absorption  coefficient  M  s  scales  as  X1  ,  is  not  applicable  in  the 

soft  x-rsy  regime.  The  results  in  (Henke,  et  si.,  1982)  indicate  that  the 
Imaginary  parts  of  the  refractive  indicies  of  most  materials  scale  very 

roughly  as  X*  in  the  soft  x-ray  regime. 


•>-  ‘ 


No?  if  we  denote  the  PVWM  resolution  of  the  multilayer  in  phase  units 

•s  •  we  have  after  incorporating  the  approximately  quadratic 

prfmh 

wavelength  seeding  into  eq  n-3-3 


5®  *  Zm"  «  8.0  *  id*\X  StC&  m  f.O  •  id4  (U,»$ 

TFWMH  ‘A) 


(AY 
(TL-3-U) 


Since 


dtp 


2K.d  cos  & 
k2 


d  A 


(1-3-17) 


we  have  that  the  spectral  resolution  scales  as 


s\ 


FWHH  „ 

X 


s9 


— ^  =  2.S  »  lo'4  ( Zd(fa)2 


(1-3-18) 


Prom  eg.  11-3-4,  the  number  of  layer  pairs  N  required  to  achieve 


auximum  reflectivity  is  given  by 


Thus,  i\  /X  and  N  are  ap^arcocimately  independent  of  X  and  6  for  a 
given  d.  In  general,  properties  of  multilayer  reflection  tend  to  depend 
■ore  strongly  on  the  spacing  d  than  on  X  or  6  . 

We  now  consider  the  scaling  of  the  angular  resolution  whose 
relation  to  $  a  is 


*  — -  Sin  6  SB  (u -3-20) 

FtfMM  X  fwM* 


If  ^  *  1,  we  can  use  eq.  II-3-16 

-a  Wei*)* 

S8tuuii  *  Z.S  tiO  - (inreutiiuis)  (1-3'Zl) 

wwHn  tan  v 

In  crystallographic  terminology,  J0  is  essentially  a 

rWnPf 

rocking-curve  width,  and  the  scaling  of  rocking  curve  width  with  6  given  by 
eq.  I I -3-21  is  similar  to  the  dependence  one  would  expect  to  find  in 
crystals.  (Canpare,  for  example,  fig.  13  in  (Burek,  1976)  with  fig.  II-3-3  of 
this  work.  The  point-by-point  calculation  in  fig.  II-3-3  is  a  nunerical 
solution  for  the  angular  FVJHM  of  the  magnitude  squared  of  eq.  II-2-11 . ) 


ANGULAR  WIDTH  (FWHM  In  minutes) 
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ANGULAR  WIDTH  OF  REFLECTION  PROFILE 
VERSUS  ANGLE  OF  INCIDENCE 
(Tungsten-Carbon  Multilayer) 
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U_E  ™ 


XM9 


Figure  11-3-3 
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Bw»vtr ,  crystals  are  available  only  with  fixed  2d- spacing*;  the  aet  of 
crystal  2d-spacings  available  in  the  soft  x-ray  regime  is  particularly 
limited.  It  is  therefore  of  interest  to  consider  the  effect  of  a  variation 
in  2d- spacing  with  wavelength  held  constant;  this  structural  variation 
represents  a  degree  of  freedom  not  allowed  in  crystalline  reflection. 

By  setting  in  aq. 1-3-21 

u  -  STS’ 

we  have 


This  scaling  law  implies  a  large  acceptance  angle  in  the  limits 
6  »>  0*  and  fl  — >90*.  (Bj.  II-3-23  has  a  singularity  at  the  limit 
points;  our  theory  does  not  apply  in  the  grazing  regime,  and  eq.  II-3-23 
must  be  modified  to  apply  in  the  normal  incidence  regime) . 

The  large  acceptance  angles  at  the  two  extremes  can  be  explained 
physically  by  the  following  arguments.  (In  fact,  we  will  see  that  a  certain 
symmetry  exists  between  the  t vk>  regions) . 


* 


II-3-14 


Pirst,  in  the  grazing  regime,  oq.  II-l-ll  isplies  a  linear  relationship 
between  the  wavelength  and  the  angle  {  to  the  surface  (  £  >  10*-  8  ) .  If  we 
set  for  the  grazing  regiaw 


s£tnJSL  se  c  un.o  se  *  -il 
x  u 


-ii. 


then  the  linearity  is  reflected  explicitly  in  the  relation 


(1-3-25) 


Further,  the  width  of  the  high-reflectivity  stopband  in  phase  units  is 
of  order  A/  1  ;  this  is  because  A/  is  the  nunber  of  layer  pairs 
participating  in  the  reflection  process.  Thus  the  resolution  of  the 
multilayer  is  governed  by  the  familiar  rule 


so  that 


(1-3-26) 


(1-3 ’2V 
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The  absorption  within  each  aell  is  given  by  the  product  of  the 
absorption  coefficient  with  the  path  length  through  the  cell,  so  that  the 
scaling  of  the  absorption  per  aell  with  angle  oust  be  of  the  fora 


4*  A*  d  _  A* 

A  cos  6  cm*  6 


(ff-3-2f) 


# 

Here  A  represents  the  spatially  averaged  imaginary  part  of  the  index 
of  refraction. 

Since  the  absorption  per  cell  is  approximately  the  reciprocal  of  the 
nianber  of  effective  layer  pairs  N, 


N  ~ 


(1-3-11) 


and  so  tie  have 


ALALJL 

X  ~  t,  ~  I* 


(JL-3-3D) 


11-3-16 


Therefore  SX/X  end  both  rise  as  f,  — >  0,  but  iX/X 

rises  more  rapidly. 

As  discussed  in  sec.  II-l-B,  our  formalism  is  invalid  if 


*  S  I /S' I 


(1-3 -31) 


In  the  near-gra2ing  regime  where  our  formalism  is  still  valid,  the 
dominant  factor  in  the  scaling  of  sultilayer  reflection  properties  is  the 
quadratic  dependence  of  the  absorption  per  cell  on  (  .  The  dependence  is 
quadratic  because  as  (  is  decreased,  the  photon  path  beccmes  more  oblique, 
and  also  d  must  be  increased. 

While  N  depends  quadratically  cm  angle  in  the  regime  near  grazing 
incidence,  it  becomes  almost  independent  of  angle  as  the  normal  incidence 
regime  is  approached,  because  the  paths  traversed  through  the  layers  beocne 

insensitive  to  angular  variations. 

-  • 

For  8  <■  45  ,  we  can  therefore  set 


X! 


A 


(1-3-32) 
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Similarly,  in  the  normal  incidence  regime  the  linear  dependence  e£  K 
on  8  disappears;  from  eg.  11-1*11 


5  (cos  6) 
cos  8 


~  ten  6  S 6 


(1-3-33) 


assuming  that  we  are  outside  the  regime  very  near  normal  incidence 
(i.e.  outside  the  regime  6  <  J JJL  ).  For  6  <  45*  we  can  set 


tan  6  —  Q  so  that 


—  8  56 


(1-3-34) 


Using 


i  • 

T  ~  77  ~  A 


(1-3-35) 


we  therefore  have 


(1-3 -36) 


1  • 
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which  is  the  limiting  form  of  eq.  11-3-23. 

Thus,  the  near- independence  of  path  length  an  angle  near  the  normal 
incidence  regime  implies  a  slower  change  in  phase  thickness  with  angle  of 
incidence,  leading  to  the  large  predicted  by  eq.  11-3-23  or 

eq.  11-3-36. 

Within  the  normal  incidence  regime  itself,  i.e.  when  0  <  JaP  so  that 
eq.  11-3-33  no  longer  holds,  we  can  set 


X 


S  (cos  6) 
cos  & 


(SB)' 


(1-3-37 ) 


so  that 


se 


(l  -3-3t) 


or  more  specifically 


se 
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Thus,  collection  angles  of  order  IS*  can  be  obtained  at  wavelengths 
around  100A.  Multilayers  may  therefore  sake  it  possible  to  achieve  the  high 
resolutions  that  are  potentially  possible  at  x-ray  wavelengths  given  large 
miner i cal  apertures;  for  this  purpose  Lwas  et  al.  (1982)  have  constructed  an 
x-ray  microscope  based  on  a  Sehwartsehild  configuration  that  is  presently 
undergoing  testing  at  the  Brookhaven  storage  ring. 

At  shorter  x-ray  wavelengths,  the  X*  factor  in  the  nmerator  of 
eq.  u-3-23  in  effect  precludes  single-substrate  reflecting  devices  from 
having  large  collection  angles.  Further,  the  necessity  of  having  layer 
thicknesses  of  at  least  a  few  angstroms  precludes  operation  in  the  normal 
incidence  regime;  in  order  to  obtain  fields  of  view  that  are  as  large  as 
possible  it  is  therefore  necessary  to  operate  at  quite  small  angles  to  the 
substrate  surfaces  (typically  of  order  a  few  degrees  or  so) .  A  device 
designed  along  these  lines  will  be  discussed  in  chapter  III. 


Section  IZ-4  -  Non-Periodic  Multilayer  Designs 


In  the  optical  regime,  the  familiar  quarterwave  stack  is  the  most 
cannon  reflective  coating,  but  more  ocnplieated  designs  are  often  used. 

In  addition  to  periodic  designs  based  on  a  thickness  ratio  other  than 
quarterwave,  such  non-quarterwave  multilayer  designs  include  structures  that 
are  not  periodic. 

Non-periodic  designs  are  resorted  to  for  a  variety  of  reasons,  such  as 
to  increase  ban&ridth  over  the  quarterwave  value,  or  to  introduce  relative 
phase  changes  between  different  polarisations.  He  will  restrict  ourselves 
to  the  question  of  whether  or  not  aperiodic  x-ray  multilayer  designs  can  be 
used  to  reduce  the  effect  of  absorption  and  so  increase  the  reflectivity. 

Other  design  objectives  will  be  left  for  future  research;  however  we 
note  that  in  the  soft  x-ray  regime  there  is  a  fixed  budget  of  layers  that 
are  permitted  in  a  design  due  to  absorption,  and  often  this  fixed  number  of 
layers  is  inadequate  to  produce  a  high  reflectivity  at  a  single  wavelength 
even  in  the  absence  of  other  design  constraints.  Techniques  for  modifying 
periodic  designs  that  require  large  increases  in  the  rwnber  of  layers  used 
might  not  be  appropriate  in  the  x-ray  regime. 


[ 
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Carniglia  and  Apfel  (1980)  have  developed  aperiodic  multilayer  designs 
that  yield  higher  reflectivity  in  the  presence  of  absorption  than  can 
periodic  multilayers  containing  the  sane  nimtaer  of  layers.  We  will  consider 
whether  such  gains  are  possible  in  the  x-ray  regime. 

Our  first  result  in  this  connection  is  obtained  in  Appendix  8.  There 
we  apply  our  difference  equation  in  the  case  of  a  periodic  reflector  that 
has  been  optimized  according  to  the  formulas  of  sec.  I1-2-B,  and  prove  that 
if  the  reflector  contains  a  sufficient  number  of  layer  pairs,  its 
reflectivity  will  be  an  extremum  with  respect  to  an  arbitrary  variation  in 
structure  (the  variation  used  need  not  leave  the  structure  periodic) .  A 
periodic  multilayer  coating  having  an  optimized  structure  therefore  provides 
at  least  a  local  extremum  in  the  reflectivity. 

Limited  numerical  studies  suggest  that  this  local  extremum  reflectivity 
is  in  fact  a  global  maxummt.  Our  most  detailed  set  of  nonerical  results  is 
shown  graphically  in  fig.  11-4-1.  Curve  1  shows  attainable  reflectivity 
versus  timber  of  layer  pairs  when  all  23  layer  thicknesses  are  treated  as 
independent  degrees  of  freedom.  Curve  2  shows  reflectivity  versus  J  when 
the  reflector  is  periodic;  in  this  curve  the  t%rc  degrees  of  freedom  in  the 
basic  period  are  reoptimized  at  each  J  value.  Curve  3  shews  the 
reflectivity  of  a  periodic  multilayer  that  has  been  optimized  as  discussed 
in  sec.  II-2-B  to  have  maximisn  reflectivity  in  the  steady-state  limit.  All 
curves  are  obtained  using  the  nor  x-ray  characteristic  matrix  solution  of 
•q.  11-2-22  with  a  canned  optimization  routine. 
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COMPARISON  OF  OPTIMIZATION  SCHEMES 


1.  R(  J)  of  aperiodic  designs  optimized  at  each  J 

2.  R(J)  of  periodic  designs  optimized  at  each  J 

3.  R(  J)  of  a  periodic  design  optimized  at  J  -»  * 


In  the  ateady-state  limit  *11  three  curve*  converge.  Even  when  the 
number  of  layer  pairs  is  somewhat  less  than  1/Im( 2 ) ,  the  highest 
reflectivity  that  can  be  obtained  frcn  an  aperiodic  design  is  not  very  such 
larger  than  that  obtained  from  an  optimized  periodic  design.  However, 
though  the  difference  in  R  between  the  periodic  and  aperiodic  designs  is 
small,  the  set  of  parameters  ^ in  the  multilayers  of  curve  1  are  in 
general  quite  different  frcn  the  constant  values  for  these  parameters  in 
curve  3,  even  in  the  limit  of  J  <x>  (but  not  in  the  joint  limit  J  •*, 
K  -^*o). 

One  would  like  to  test  these  conclusions  as  ccnprehensively  as 
possible,  but  aulti-dimensional  optimization  is  very  expensive. 

Carniglia  and  Apfel  have  solved  a  simplified  version  of  the 
2J-dimensicnal  optimization  problem  for  visible  HR  reflectors.  They  show 
bow  to  calculate  the  particular  thicknesses  for  the  top  pair  of  layers  in  a 
bilayer  stack  that  will  maximize  the  total  reflectivity  of  the  multilayer, 
given  that  the  preceeding  layer  pairs  have  already  been  optimized  in  this 
way. 

They  hypothesize  that  such  a  two- by- two  optimization  is  fully  general; 
in  other  words  that  no  increase  in  total  reflectivity  can  be  obtained  by 
simultaneously  changing  the  thicknesses  of  any  preceeding  layers  when  a 
given  pair  of  layers  is  being  optimized  using  their  algorithm. 

Carniglia  and  Apfel's  mathenatical  calculations  are  not  applicable  in 
the  x-ray  regime  since  they  assizne  that  in  the  preceeding  stack  of  layer 
pairs,  1. 


Nonetheless,  their  basic  concept  is  of  Interest  in  the  x-ray  regime 
because  mirrors  made  by  the  in-situ  reflectance  monitoring  (XSFK)  technique 
discussed  in  Chapter  Z  will  have  a  related  property.  If,  during  a 
deposition  session  like  that  sham  in  fig.  9  of  Spiller  (1962b) ,  one  cuts  off 
the  depositions  of  successive  high  index  layers  at  successive  reflectance 
maxima,  the  ^  value  for  each  layer  pair  will  automatically  be  optimized 
under  an  optimization  schone  like  Carniglia  and  Apfel's.  The  thickness 
ratio  parameters  may  be  chosen  by  a  different  criterion,  but  the  total 
thickness  with  which  each  pair  of  layers  is  terminated  will  be  the  thickness 
which  maximizes  the  total  reflectivity  of  the  entire  stack. 

In  Appendix  9  it  is  sham  that  in  the  x-ray  regime  the  optimum  , 
given  an  arplitude  reflectance  £  fran  the  stack  of  preoeeding  layers, 
satisfies  the  condition: 


Im[  a.[i  +  t(p.  e'*  n  ♦ 

*  0 

(l-4- 

In  Appendix  9  we  shew  that  the  solution  for 

q>  is 
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Sere,  as  in  nc.  IX-2-B,  w  will  asaim  that  the  bean  can  be  treated  as 
linearly  (S)  polarised.  The  qualitative  conclusions  drawn  in  this  section 
apply  also  in  the  oase  of  unpolarised  radiation,  ft*  key  results  in 
eqs.  11-4-1,2,8  and  9  apply  to  both  polarisations. 

The  optical  constants  used  in  eqs.  11-4-1,2  are  those  of  the  XSIW 
monitoring  wavelength,  which  night  be  different  fran  those  at  which  the 
coating  is  used,  for  exanple  if  the  coating  is  want  to  operate  near  normal 
incidence. 

As  K  becomes  large,  the  multilayer  produced  by  ISHK  becomes  periodic, 
with  ^  approaching  the  Darvin-Prins  solution;  in  lowest  order  eq.  XI-4-1 
becomes  in  S  polarization 


In 
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CVA) 


_ t _ 

(  AJ  tanp. 


0  (I-4-J) 


If  we  now  make  the  rough  approximation  that 


A  *  0,  A  -  A  *  A  (1-4-4) 

t  'HI  H 


we  get 
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This  yields  the  correct  Bragg  condition 


S  -  0 


(1-4-4) 


given  that  the  optimization  condition  on  ^  of  sec.  1I-2-B  is  also 
satisfied,  so  that  the  sun  of  the  first  and  last  terms  in  eq.  11-4-5  vanish. 

However,  A^  E  0  is  only  a  rough  approximation,  so  that  in  general  the 
Bragg  condition  will  not  be  satisfied  under  the  two-by-two  optimization 
produced  with  ISRM. 

Thus,  in  the  x-ray  regime  the  hypothesis  of  Carniglia  and  Apfel  does 
not  hold,  and  a  2x2  optimization  scheme  based  on  an  L-H  unit  cell  does  not 
yield  a  fully  optimized  reflector;  as  a  non-optimality  in  the  Bragg 

condition  the  error  may  be  considered  to  be  equivalent  to  a  small  shift  in 

X  or  6  .  As  a  small  X  or  0  shift,  this  non-optimality  is  not  in  itself 

crucial  since  it  tan  be  compensated  for  by  an  intentional  mall  shift 

between  the  monitoring  X  or  0  and  the  X  or  0  of  the  coating's  intended 
application. 

Bowever,  a  similar  problem  arises  when  we  consider  a  two- by-two 
optimization  scheme  for  the 

In  general  the  in  an  ISRM  multilayer  could  be  chosen  in  a  variety 
of  ways.  For  example,  all  oould  be  set  equal  to  the  fi  of  an 

optimized  periodic  multilayer. 


xx-4-e 
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Alternatively,  one  oould  terminate  deposition  of  t he  low  index  layers 
at  reflectance  minima.  This  turns  out  to  prodooe  values  that  are 
miner ically  quite  close  to  cptixun. 

Finally,  one  oould  choose  each  p*  in  accordance  with  a  2x2 
optimization  scheme,  choosing  each  pm  successively  to  maximize  the 
reflectance  of  the  stack  of  cells  1  ...  K,  as  suggested  by  the  discussion  in 
(Spiller,  1976). 

As  shown  in  Appendix  9,  each  under  such  an  optimization  scheme  is 
given  by 


1  r  r  _ l/\H*)\ _ 
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where  A  ■  sin  V^))  , 
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which  is  the  solution  to  the  condition 


In  { V  *?(».♦  j- )-  f  ( V  \)  (e.  e,ip*  *  ~tk  ) 


(1-4-1) 

Xn  the  steady-state  regime,  eq.  XX-4-9  beccnves  in  lowest  order  (for  S 
polarization) 

1,n{v  ZTa  ( is *  »  0 

( 1-4-10 ) 


As  before,  this  will  not  be  consistent  with  the  results  of  sec.  II-2-B, 
unless  ve  make  the  crude  approximation  A  S  Q  .  Thus,  except  to  a  rough 
approximation,  the  2x2  optimization  scheme  converges  to  a  nan-optimal  value 
of  |3  as  well  as  a  nan-optimal  value  of  ^  . 

We  can  see  an  analogy  with  the  periodic  case  by  adding  eq.  II-4-3  to 
eq.  II-4-10.  We  obtain  in  lowest  order: 
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Sere  it  has  not  been  assumed  that  £>  *  0  . 

i. 

Me  showed  above  that  under  an  ISAM  scheme,  non- negligible  absorption 
and  refraction  in  the  low  index  layer  cause  the  first  tarn  in  the  above 
condition  to  be  non-zero.  She  sum  of  the  second  and  third  terms  in  the 
above  condition  cannot  then  be  zero;  and  since  such  a  cancellation  is  also 
the  requirement  for  an  optimized  fi ,  we  can  consider  the  non-satisfaction  of 
the  Bragg  aonditian  to  have  in  turn  prevented  the  2x2  optimization  sehotie 
from  converging  to  the  oorrect  thickness  ratio. 

this  is  similar  to  the  situation  tie  found  when  treating  the 
optimization  of  periodic  multilayers  in  sec.  II-2-B,  where  we  noted  that  our 
condition  for  optimizing  ^  in  a  periodic  multilayer  applies  only  if  the 
generalized  Bragg  aonditian  Re(S)«0  is  satisfied. 

In  the  usual  case  that  absorption  and  dispersion  in  the  low  index  layer 
is  snail,  the  runerical  consequences  of  these  non-optimalities  is  only 
moderate,  since  the  accuracy  of  eq.  11-4-4  increases  as  the  index  decrement 
in  the  low  index  layer  goes  to  zero.  K  $7.6h,  0*0*  reflector  with  J*200 
would  have  a  reflectance  of  .19  if  designed  by  the  2x2  procedure,  ctmpared 
with  a  .22  reflectance  from  the  optimum  periodic  design,  for  a  A R/R  of  15%. 

Further,  the  £  values  and  the  monitoring  conditions  can  be  chosen  in 
such  a  way  as  to  aonpensate  for  this  reflectivity  loss.  For  exanple,  if  the 
monitoring  and  application  wavelengths  were  the  same  one  eould  introduce  a 
shift  between  the  monitoring  angle  and  the  application  angle  that  would 
cancel  the  difference  between  aq.  11-4-2  and  eq.  II-2-30;  one  could  then 
choose  the  ^  under  one  of  the  alternative  schemes  described  above. 
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The  main  advantage  of  the  XSRM  technique  lie*  in  its  potential  for 
strongly  reducing  the  effect  of  randan  thickness  errors,  rather  than  in  its 
use  as  a  method  for  selecting  the  target  thiefoesses. 
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Section  II-5  ~  Effect  of  Aocunulating  Randan  Thickness  Errors  an 
Multilayer  Reflectivity 

Part  A)  Introduction 

X-ray  multilayers  nut  contain  a  large  number  of  layer  pairs  since 
the  reflectivity  of  individual  layers  is  very  small.  The  magnitude  of 
uncorrected  thickness  errors  in  the  multilayers  nut  therefore  be  kept  to 
a  small  fraction  of  a  layer  thickness  in  order  to  prevent  a  large 
cumulative  dephasing  from  developing  between  the  top  and  bottom  of  the 
stack.  Because  the  layers  themselves  are  only  of  order  ten  angstroms  in 
thickness,  the  tolerance  an  random  errors  in  the  layers  is  extrmnely 
tight. 

The  two  fabrication  methods  described  in  Chapter  I  can  be  regarded  as 
alternative  approaches  to  solving  this  problem. 

Barbee's  vacuum  sputtering  technique  (Barbee,  1982)  achieves  precise 
control  of  the  layer  thicknesses  through  the  use  of  sophisticated 
deposition  procedures  and  apparatus;  these  include  use  of  a  magnetron 
sputtering  sourae  to  prevent  high  energy  ion  bcmbar&sent  of  the  substrate 
surfaoe,  precise  control  of  the  deposition  voltage,  and  use  of  a 
multi-angle  preparatory  pre-deposition  onto  the  substrate  surfaoe. 

The  second  fabrication  method  is  Spiller's  in-situ  reflectance 
monitoring  technique  (ZSHM),  which  serves  to  eliminate  the  accumulation  of 
dephasing  fresr  any  thickness  errors  that  may  be  present.  The  intensity  of 
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the  XKSM  probe  shows  an  interference  oscillation  as  the  thickness  of  the 
top  layer  increases;  this  oscillation  is  essentially  determined  by  the 
absolute  difference  in  phase  between  the  tf£)er  surface  of  the  bop  layer 
and  the  stack  of  preceeding  layers,  the  maxima  in  the  peak  intensities 
occur  when  the  double-pass  phase  thickness  of  the  «$per  pair  of  layers  is 
approximately  one  wavelength;  the  precise  single-pass  phase  thickness  that 
is  required  is  given  by  eq.  11-4-2. 

8q.  11-4-2  obtains  even  when  the  reflectivity  ^  of  the  preceeding 
stack  is  nan-optimal  due  to  thickness  errors.  Under  IS  PM,  the  target 
thickness  for  the  qpper  layer  of  the  final  unit  cell  (aimed  at  by 
monitoring  the  1SRM  probe  intensity) ,  is  that  thickness  given  by 
eq.  II-4-2  which  puts  the  top  interface  of  the  layer  into  phase  with  the 
preoeeding  layers,  thereby  compensating  in  large  part  for  any  errors  in 
those  layers. 

Multilayers  of  the  first  kind  are  subject  to  what  we  will  call 
accumulating  thickness  errors,  while  those  fabricated  by  1SHM  are  subject 
to  what  we  will  call  non-accumulating  thickness  errors.  He  will  analyze 
the  effect  of  accumulating  errors  in  part  B,  and  will  disoiss  our  solution 
in  part  C.  Hen-accumulating  errors  will  be  discussed  in  eec.  II-6. 
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Part  B)  Analysis  of  Accumulating  Thickness  Errors 

Thickness  errors  of  an  accumulating  type  have  been  treated 
theoretically  by  Shellan  et  al.  (1978) ,  and  have  also  been  discussed  briefly  by 
Bauraeister  (1981).  Shellan  presents  solicit  closed-form  solutions  for 
the  reflectivity,  but  acme  of  his  approximations  are  not  appropriate  to 
the  x-ray  case. 

Shellan  treats  non-absorbing  multilayers  using  coupled  wave  equations 
that  contain  perturbation  terms  to  represent  the  effect  of  thickness 
errors.  Shellan' s  coupled  wave  equations  are  of  a  form  that  apply  to 
dielectric  structures  operating  at  the  dielectric  Bragg  aonditicn.  The 
advantage  here  is  more  than  the  avoidance  of  complex  indicies  of 
refraction;  if  the  nominal  phase  thickness  of  the  unit  cell  is  both  pure 
real  and  equal  to  X  ,  the  aoupled  wave  solutions  simplify  considerably. 

Shellan 's  analysis  assumes  a  small  coupling  constant  per  cell,  as  is 
necessary  when  a  large  number  of  layers  must  participate  in  the 
reflection;  in  this  feature  his  analysis  is  appropriate  to  the  x-ray 
regime.  In  the  limit  that  the  perturbation  terms  are  moall,  Shellan 
solves  rigorously  for  the  expectation  value  of  the  reflectance.  Be  also 
obtains  an  expression  for  the  reflectivity  in  the  presence  of  larger 
errors,  in  the  ease  of  structures  which  contain  sufficiently  few  layers 
that  their  overall  reflectivity  la  nail  oanpared  to  one;  to  do  this  he 
uses  what  he  calls  an  'undepleted  incident  wave'  approximation.  In  the 
context  of  x-ray  diffraction,  auch  structures  might  be  said  to  be 
operating  in  the  'Fourier  transform  regime'. 
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X-ray  reflector*  necessarily  have  a  small  coupling  constant  because 
the  x-ray  indicias  of  refraction  of  all  materials  are  cloae  to  one. 

However,  in  the  x-ray  oase  a  significant  fraction  of  the  coupling  between 
the  structure  and  the  radiation  field  is  due  to  absorption.  Also,  while 
the  reflectivity  of  x-ray  multilayers  can  be  quite  severely  limited  by 
absorption  (leading  to  a  situation  where  a  relatively  small  nutber  of 
layers  participate  in  the  reflection  prooess) ,  it  is  nonetheless  usually  a 
poor  approximation  to  apply  the  vsidepleted  incident  wave  approximation  to 
x-ray  multilayers. 

It  is  also  of  particular  interest  in  the  x-ray  case  to  be  able  to 
consider  operation  off  the  Bragg  condition,  since  x-ray  multilayers  have 
very  narrow  baivVidths. 

Another  consideration  in  the  x-ray  case  is  that  while  present 
technology  permits  the  fabrication  of  multilayers  having  2d  spacings  of 
100A  and  structures  of  fairly  good  quality,  the  soft  x-ray  spectrum  can  be 
considered  to  extend  to  normal  incidence  wavelengths  that  are  shorter  than 
this  by  an  order  of  magnitude;  in  addition  aone  experiments  have  been 
carried  out  in  higher  orders,  which  are  in  a  sense  equivalent  to  still 
shorter  2d- spacings.  We  will  see  below  that  the  tolerance  on  accumulating 
thickness  errors  has  an  approximately  quadratic  scaling  with  2d- spacing. 

We  therefore  consider  it  desirable  to  be  able  to  analyze  the  effect 
of  thickness  errors  that  are  large  enough  in  ccnpariscn  with  the  layer 
thicknesses  to  substantially  degrade  the  reflectivity;  in  this  case  a 
perturbation  treatment  is  inappropriate. 


-  ,-i  -  . 


II-5-5 


For  these  rtuoni  we  have  developed  a  new  approach  to  treat  the 
affect  of  randan  thickness  errors  an  Multilayer  reflection.  Our  approach 
renders  the  problan  of  thickness  errors  tractable  by  using  an 
approximation  that  applies  rigorously  in  both  the  limit  of  large  errors 
and  the  limit  of  snail  errors.  approximation  also  tends  to  be  quite 
accurate  in  the  intermediate  region. 

Our  analysis  is  based  an  a  deuanposition  of  the  wplitude 
reflectivity  into  what  are  essentially  coherent  and  incoherent  parts  (see 
aq.  II-5-10  below) . 

Our  key  approximation  is  to  neglect  cubic  and  higher  pesters  of  the 
incoherent  reflectivity.  It  is  essentially  this  approximation  that  allows 
steady-state  solutions  for  the  overall  reflectivity  to  be  found 
analytically  (see  eqs.  II-5-35  and  41). 

With  further  approximations  we  find  solutions  for  the  reflectivity 
outside  the  steady-state  regime  (this  analysis  is  carried  out  in  Appendix 
10;  the  results  are  shown  in  aqs.  11-5-42  and  45) . 

While  we  use  a  rumber  of  approximations  in  the  derivation,  it  is  the 
neglect  of  higher  order  powers  in  the  incoherent  reflectivity  that  is  the 
key  approximation  in  two  senses.  First,  this  is  the  approximation  to 
which  our  results  are  nsaerically  most  sensitive.  Second,  the  other 
approximations  that  we  use  are  akin  to  our  key  approximation  in  that  they 
allow  our  results  to  be  correct  in  both  the  limits  of  large  and  small 
thickness  errors.  However,  the  other  approximations  generally  relate  to 
specific  algebraic  equations  arising  in  the  derivation. 
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We  will  use  a  particular  r-ray  Multilayer  aa  an  axatple  while 
presenting  our  analysis,  namely  a  twgsten-carbon  multilayer  reflecting 
67.6A  radiation  (B-K^ )  at  normal  incidence,  the  layer  thicknesses  in  the 
absence  of  thickness  errors,  d  w  •  7.6A,  d  »  26. SA,  are  chosen  in 
accordance  with  the  optimisation  conditions  of  sec.  ZI-2-B.  Zn  cur 
exanples  we  will  generally  take  J  •  200,  which  is  a  sufficiently  large 
nimter  of  layer  pairs  to  plaoe  the  structure  in  the  steady-state  regime. 

Following  Shellan,  we  have  tested  the  accuracy  of  our  method  using  a 
Monte  Carlo  simulation,  in  which  a  oorputer  progrwn  is  used  to  generate  a 
large  matter  of  simulated  multilayer  stacks  containing  randan  errors,  and 
to  oonpare  their  mean  reflectivity  with  the  analytical  prediction. 

We  have  found  that  the  analysis  gives  accurate  results  over  s  full 
range  of  WE  error  heights,  so  long  as  the  mdegraded  reflectance  of  the 
structure  aider  consideration  does  not  approach  too  closely  to  one. 
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We  will  calculate  the  aean  reflectivity  of  multilayers  that  contain 
uncorrelated  errors  in  the  thicknesses  d(  of  the  different  cells.  We 
will  assune  that  the  nominal  Multilayer  structure  is  periodic,  and  that 
the  statistical  properties  of  the  errors  are  the  same  in  all  layers.  If 
d  is  the  nan  thickness,  then  let 
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(1-5-1) 


where 


<Ad(>  *  0 

«Ad/>  e  cr* 


(1-5*2) 


We  will  treat  the  variance  O'  as  a  known  quantity.  If  each  cell  consists 
of  an  L  layer  and  an  H  layer,  each  of  which  has  independent  errors  with 
variances  <AL*  >  and  <AH*>  ,  then  due  to  the  independence, 
a‘  «  <  AL\>  ♦  <AH*>  ,  and 
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^c«*e(<AL*>  ♦  <AH*>) 


4a‘cr* 

A4 


c os*  9 


(l-5-d) 


XI-5-8 


where  from  eg.  II-l-ll 


1x6 


9  *  >  ♦  Af  *  K  ~  — r —  cord  +  Af 


ZxAd* 

<,>  +  —r-‘*'>  „.s,4> 


With  these  definitions  eq.  n-1-20  becomes  in  first  order 


k«*i 


-r.Af 

e  e  *  ^  -  f«'r- p)- 


We  note  that  in  eq.  II -5-5,  r,  p,  end  >a  are  treated  as  independent  of 
<p  ;  this  approximation  is  accurate  to  within  terms  of  order  ^A  .  As 
discussed  in  sec.  II-l-B,  we  can  neglect  terms  of  order  f  ■  A  and  A  in 
analyzing  stochastic  problems  like  the  present  one,  but  we  must  retain 
high  orders  in  the  parameter  7  . 

The  fact  that  we  can  directly  incorporate  aq.  II-5-4  into  aq.  II-5-5 
ia  die  to  an  advantage  of  the  difference  aquation  formalist;  nsnely  that 
the  formalism  directly  incorporates  the  physical  pre^erties  of  the 
individual  cells. 

The  fact  that  our  equations  use  the  properties  of  the  different  cells 
explicitly  has  another  useful  consequence.  Because  is  determined  by 
the  properties  of  the  cells  K-l,  K-2,  ...  ,  but  not  by  those  of  cell  K 
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(»•«  mc>  II-l-B),  we  can  writ®  for  arbitrary  functions  f  and  g, 


and  in  particular, 


»2i  A| 

<«  -  <  e  *'«>  <?  > 


(*-5-7) 


The  first  expectation  value  in  «q.  II -5-7  is  the  characteristic  function 
of  the  variable  (Papoulis,  1965,  p.159) .  Any  probability 

distribution  with  known  characteristic  function  can  be  used  in  the 
analysis  to  folio;.  For  simplicity  we  will  assume  Gaussian  statistics;  as 
we  will  see,  the  final  results  are  very  insensitive  to  the  exact  choice  of 
distribution  function. 

Then  using 


•  Si  At  ^ 

<  e  T*> 


.Ka^> 


which  holds  for  a  Gaussian  randan  variable  with  tero  Bean,  we  take  the 


; 
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expectation  value  of  eq.  II -5-5,  and  get 


<f  > 


C  *  <  ^  >  -  f  ♦  p)  -(iMf)  <  f  > 


(1-5-1) 


Now,  we  define  the  variable  ^  through  the  relation 


^  ■  <  £  >  ♦  £  (1-J-JO) 


p  therefore  represents  the  random  or  nan-deterministic  pert  of  the 
reflectivity.  Since  the  phase  of  ^  is  ccnpletely  randan,  we  Bay  also 
regard  it  as  the  incoherent  oomponent  of  the  reflected  bean. 

Since 


<?*>  -  <<*>  ♦  *<?.>  ?,♦  **  > 
-  <•  >  ♦  <f‘> 

'l  » 
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eq.  H-5-9  becomes 
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-  (ir+p)<^>  -  (;r*p)<  ^*> 


In  order  to  solve  this  equation,  we  mat  obtain  <  S*  >  . 

*  K 

Subtracting  aq.  II-5-12  from  eq.  II-S-5, 


e  e 


e  e 


<*  > 

'  K 


-Z.<t>  -i.Af 

♦  «  c  “o 

'  ft 
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Next  we  square  eq.  II-5-13 ,  neglecting  terms  of  order  -  A  and  A 
He  find 
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Defining 


?  *  -  <  t  >  -  (T-  ip)  >  (1-5-15) 


we  can  combine  the  second  and  the  second  to  last  terns  cn  the  right  of 
eq.  11-5-14  to  obtain 


.  f  -4;<t>  -4i<t>  -4<a**>  .  -4;<«>  -*<*♦*>  •  *•'*♦*  .  / 
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the  error  in  this  substitution  is  of  order  f  •  A . 


Talcing  expectation  values,  and  using  <  |  >  s  0  ,  we  get: 
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He  will  first  find  steady-state  solutions  to  sqs.  11-5-12,17.  In 
order  to  operate  in  the  steady-state  limit,  a  multilayer  Bust  contain  a 
mmber  off  layers  sufficient  to  extinguish  the  incident  beam  before  it 
reaches  the  substrate.  In  addition,  the  statistical  properties  of  the 
thickness  errors  must  be  independent  of  K.  In  general,  extinction  may 
result  from  a  combination  of  absorption,  dephasing  due  to  detuning  from 
Bragg  resonance,  dephasing  due  to  thickness  errors,  and  depletion  by 
back-reflection. 

In  the  soft  x-ray  regime,  absorption  is  a  somewhat  more  important 
extinction  mechanism  than  depletion  by  back-reflection,  so  that,  following 
the  discussion  off  sec.  II -3,  we  may  consider  a  multilayer  to  operate  in 
the  steady-state  regime  if 
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Following  the  discussion  of  sec.  II-l-B,  we  can  consider  the  total 
phase  error  contained  within  the  structure  of  s  multilayer  operating  in 
the  steady-state  limit  to  be  off  order  ,  assuming  that 
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the  thickness  errors  accumulate  in  »  random  walk  fashion. 

When  ^  >u*  ,  the  errors  may  be  treated  as  a  mall 

perturbation  cn  the  ideal  periodic  structure,  and  nay  be  considered  to 
cause  mly  a  mall  drop  in  reflectivity.  (This  will  be  shown  aore 
rigorously  below). 

On  the  other  hand,  if  <A^  m  ,  the  errors  cause  a  substantial 
drop  in  reflectivity,  and  can  no  longer  simply  be  treated  as  a 
perturbation. 

2.  m 

In  the  soft  x-ray  regime,  the  condition  <Af  >  «  >u  is 
considerably  more  restrictive  than  the  condition  <Af*>  ^  1  .  When 

<  A  g*>  approaches  gne,  the  different  layer  pairs  can  no  longer  interact 
coherently  and  the  reflectivity  becomes  very  small,  assuming  that  the 
structure  has  a  mall  coupling  constant.  (In  our  67.6A  example  with 

<  Ag*>  set  to  1,  the  Monte  Carlo  simulation  program,  described  below  shows 

that  R  •  0.0027).  The  assmpticn  <Aj*>  1  is  therefore  an 

assimptian  that  the  structure  maintains  seme  resonant  reflectivity;  it  is 
a  generalization  of  a  standard  assumption  in  crystal  diffraction  theory, 
which  in  our  notation  states  that  g  <<C  i  . 

In  the  analysis  that  follows,  we  will  not  need  to  assume  that 
<Af*>  4^  M*  r  but  will  assume  that  <Ag*>  ^  1  ,  or  sore 
specifically  that  terms  of  order  A  g*  are  negligable.  If  in  cases  of 
greatest  interest,  <Ag  >  —  ju  »  then  terns  of  order  A  f  are 
numerically  comparable  to  terms  of  order  g  'A  . 
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The  distinction  should  be  born  in  nind  when  following  the  discussion 
in  (Shellan,  1978) ,  whose  perturbation  treatment  in  essence  bakes  the 
assunpticr  that  <Af*>4C  \u  |. 

We  should  note  that  the  above  arguments  are  somewhat  oversimplified 
in  that  they  only  consider  extinction  of  the  beam  due  to  error-induced 
dephasing  and  absorption.  We  will  see  below  how  these  criteria  are 
generalized. 

Neglecting  terms  of  order  A  9 *  ,  we  can  aet 
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independent  of  the  probability  distribution  satisfied  by  the  random 
errors. 

,  i 

Bq.  11-5-17  then  becomes  (neglecting  terms  of  order  g-  A  and  A  ) 


<o*  >  *  “4  <  A  <p  >  ♦  (l»4i$- J  <  Ag>*>)  <  > 


- Z(ir+  p)  <?’> 


(1-5-20) 
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The  key  approximation  used  in  the  present  analysis  is  based  on  the 
double  inequality 


<1^1  >  <  <1^|1>  <  1  (l-J-21) 


The  right  hand  inequality  is  fairly  strong  in  the  soft  x-ray  regime,  where 
the  absorption-limited  reflectivity  is  typically  of  order  201.  The  right 

,  a 

hand  inequality  is  particularly  strong  when  <Af  >  ju  ,  since  the 
multilayer  reflectivity  is  then  severely  degraded  by  the  randan  thickness 
errors. 

On  the  other  hand,  the  left  hand  inequality  becomes  very  strong  when 
<Aga>  4s.  **  •  *ince  in  the  limit  that  the  errors  are  snail,  <  ^  > 

approaches  p  (see  eq.  11-5-10) . 

Since  both  inequalities  are  always  in  foroe  to  a  certain  extent,  one 
can  generally  assisne  that  <  I  p  I  *>  is  amall  compared  to  one  (see 
fig.  11-5-1  which  shows  <(f  l*>  plotted  as  a  function  of  r).  We  will 
therefore  asstsne  that  terms  which  are  cubic  or  higher  in  ^  are 
completely  negligable. 

It  is  this  neglect  of  cubic  terms  in  ^  that  renders  the  problem  of 
tendon  thickness  errors  tractable  ewer  a  broad  range  of  error  magnitudes. 

These  cubic  terms  are  particularly  mall  in  the  limits 
<Af*>  4C  m"  «  <  Af*>  &  M*  •  Our  analysis  is  therefore  most 
sensitive  to  the  approximations  we  have  made  at  values  of  CT  in  the 
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•  •  *  • 

intermediate  region.  In  our  67.SA  axanple,  <  A  f  >  •  M  at  c  •  1.22A;  tat 
will  use  O’  ■  1A  when  presenting  mmerical  exanples. 

Our  key  approximation  has  been  found  to  be  quite  aocurate  for 
multilayers  operating  in  the  soft  x-ray  regime,  as  the  mmerical  examples 
presented  below  trill  show. 

If  we  neglect  terms  of  order  in  eg.  II-5-20,  we  get 


<?*  >  -  «5  l*>  t(4i?-f<Ag*>Kl6  |1>-4<Ag1>  <?  > 

'Di  ’«  '«  » 

(At  c  ■  1A  in  our  67. 6A  example,  the  neglected  cubic  terms  have  about 
0.04  of  the  magnitude  of  the  third  term  an  the  right  in  eq.  11-5-22). 

When  radiation  no  longer  reaches  the  substrate, 

<?Z  >  =  <  ?2  >  *  <  oZ>  so  that  eq.  XI-5-22  has  the  solution 


<?*> 


<  A  g  .>  <  ^  > 
“iS  ♦  2  <  A  ^  > 


(I -J*X3) 


A 
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In  the  steaty-state  this  solution  can  be  substituted  into  the  last 
term  of  sq.  XI-5-12,  which  then  becomes 


<*> «  e  e  <*>  -  (ir-p) 


•  (<r  ♦  p)  <p>  (l - — -  — - \ 

\  i<t>  *2<A^>  ♦  (ir  +  p)  ^ 


or,  neglecting  terms  of  order  q.  A  and  , 


<*>  Z(i<t>  ♦  <Afa>)  *  -  (ir-p) 

-(ir*p><*>  ^1"  ^ 


(l-J-2*; 


This  equation  reduces  to  a  cubic  polynomial  in  <  ^  >  . 

We  note  that  the  substitution  of  eq.  11-5-23  into  the  last  term  of 
eq.  Z 1-5-12  tends  to  be  quite  accurate  in  the  soft  x-ray  regime,  for  three 
reasons.  First,  |(ir  ♦  p)  |  is  generally  somewhat  less  than 
|2i<t>|  ,  so  that  our  approximation  is  made  in  one  of  the  less 
sensitive  terms  of  eq.  XI-5-12.  Second,  the  factor  <'$1>  in  the  fourth 
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X 

term  is  snail  oonpared  to  the  factor  <  ^  >  in  the  third  tarn,  except  when 

<  A  9*>  2>  m'  i  in  which  case  both  the  third  and  fourth  tens  of  the 

equation  become  nail  oaipared  to  the  second  tern.  In  our  67. 6A  txsrple 

the  last  tern  of  eq.  II-5-12  is  0.22  of  the  magnitude  of  the  second  tern 

on  the  right-hand  side  at  c  *  1A.  Finally,  the  cubic  terms  in  "p  that 

are  neglected  in  eq.  II-5-23  tend  to  be  quite  snail  as  discussed  above. 

Using  similar  approximations,  eq.  XX-S-25  can  now  be  converted  from  a 

cubic  equation  to  a  quadratic  aquation  with  very  little  loss  in  accuracy. 

I 

In  the  soft  x-ray  regime  <  p  >  is  fairly  snail  even  when  the  errors  are 
snail,  making  the  aeoond  to  last  tern  of  eq.  11-5-12  snail.  The  last  tern 
of  eq.  II-5-12  tends  to  be  acnewhat  snail  for  reasons  discussed  in  the 
preoeeding  paragraph.  The  steady-state  solution  of  eq.  II-5-12  when  these 
two  terns  are  neglected  is 


<?> 


[cr*ip)/2<t>  ] 

1  ♦ '  <A'Vi<t> 


U-5-Z6) 


(Prom  this  one  can  derive  the  expression 


!<?>!* 


iy+;piV* _ 

(<f  >  “  ju)  *  (<  m")* 
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which  is  quite  useful  foe  rough  estimations) .  The  argument  used  in 
deriving  sq.  II-2-26  can  be  considered  to  be  somewhat  analogous  to 
Shellan's  "undepleted  incident  wave"  approximation. 

Eq.  11-5-26  is  least  accurate  when  <At*>  «  jf ,  since  the 
neglected  second  to  last  tern  of  sq.  XI-S-12  can  then  approach  the  sane 
order  of  magnitude  as  the  second  tern. 

Bq.  II-5-26  can  be  normalised  to  give  the  correct  &  •  0  value  by 
re-writing  it  as 


<?> 


(t-S-lt) 


where  the  numerator  is  the  solution  to  eq.  11-5-5  when  O’  ■  0,  namely 
eq.  II-2-11. 

The  following  considerations  permit  a  further  improvement  in  the 
accuracy  of  sq.  11-5-26  when  O’  is  mnall.  If  O’  is  mall,  the  factor! 
in  the  denominator  of  eq.  11-5-23  can  be  set  equal  to  i  (given  by 
•q.  11-2-27) ,  since  the  left  hand  inequality  in  eq.  11-5-21  beocmes  very 
strong,  permitting  u  gleet  of  terms  of  order  (<  A />/,*')  .  Bq.  II-5-24 
then  beams 

<(>  {!.'<«>  *Z<  Afl>  «  -  (ir-f)  -  (irtp)<(>* 

i  <a/>  * 

-  (iT*p)  >  -  (  iT  ♦  p)  - — -  <?> 

(z-r-zf) 
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where  we  have  set  -  i  J  *  l  <  Af*>  «  -  it  in  the  lest  ten, 

since  we  are  foe  the  ament  considering  ( <  A  f*>  /  m  ")  to  be 
anall. 

If  we  define  if  B  < 9 >  ~  p  ,  and  neglect  terms  of  order 
eq.  II-5-29  beoomes 


Ap2i<t>+2<Af>p  *  “2^irtp)jAo 

•  •  ' 


<*t*>  2 

hr  +  p)  — —  9 
li  0 


(t'S-30) 


which  with  the  use  of  eq.  II-2-11  reduces  to 


<?>  -  (1-  •' <*f>  (- )) 

(This  may  be  regarded  as  a  second-order  perturbation  solution  to 
eq.  II-5-9.)  Writing 


I I -5-23 


we  can  rewrite  eq.  n-5-31  as 


<*> 


3-<t> 

2<x 


)  <&?*> 


(l-s-ii) 


bo  that  it  now  has  the  correct  limiting  behavior  when  XT  it  mall. 

Bqs.  11-5-26,33  are  plotted  in  fig.  11-5-2  for  the  case  of  our  67. 6A 
example. 

Substituting  eq.  11-5-33  into  the  last  term  of  eq.  11-5-25  (note  tha 
this  last  term  is  most  significant  at  mall  XT,  where  eq.  II-5-33  is 
designed  to  be  most  accurate) ,  we  obtain  the  quadratic 


<9>  (r-  ip) 


1- 


•  <t>  ♦  Z.<  A  j*>  - 


i<S  ♦<»> 


\ 


4  2  <p  >  (<t>  -  '  <  Aj*>)  ♦  (T*  ip)  *  0 


(1-5-34) 
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whose  solution,  after  sanipulaticn,  is  found  to  be 


<? >  *  - 


('-•pi  <At*>^i-^lS.j[<»>.;<AT1>  y-,t  j 


*  < 


(<t>  -  \<L)>)  4 


(<t>-i<Af*>)Z 


^  .  J.  (. 

r*-<*>ir  *' 

A  ,1 

<  A  f  >  1  2 

H 

<t>4i<Af  > 


(I-S-3S ) 

This  is  our  final  solution  for  <  p  >  in  the  steady-state  regime;  it  is 
plotted  in  fig.  II-5-2. 


K>(s» 


<R>  - 


<  I  <*>!*♦  2 Re  (<p> ♦  lpl*> 


-  !<?>!*  +  <lpl*> 


we  need  to  find  <1  p  l  ,>  in  the  steady-state  in  order  to  find  <  R.  >  . 

Taking  the  magnitude  squared  of  eq.  I I -5-13,  and  neglecting  terms  of 
order  ^  •  A  and  A  gives 

Ip  I  ®lpl  -A^A<fZ  *  )}|p  |  ♦4A^1|p|  +4A^2|<p>l 

*tA  t‘ge  ?B  )  4  w»  ((<4  f*>  -  A  f*-  i  <f>  t * J 

-2*4ir4rHl»/?.- ?*<?*>)) 
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We  have  neglected  terms  of  order  (A*1),  and  have  used  eq.  11-5-15. 

The  final  term  an  the  right  of  eq.  11-5-37  is  of  order  ^  ,  and  we 
will  neglect  it  for  the  reasons  discussed  above. 

Taking  expectation  values: 


<1*  |*>  *=  <|g  \t>-4l"<\Z  |a>  +  4<A^>|<0  >|* 

«♦*  «  *  '«  'x 

(I-5-if) 


In  the  steady-state 


<|f  |  >  *  *  <lf  I  > 


(1*5-31) 


so  that 


<•?!*>  -  ^  K»>l* 


(1-5-40) 
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The  neglect  of  cubic  tens  is  s  poorer  approximation  hers  then  it  was 
in  the  steps  between  eqs.  II-5-20  and  23,  since  in  the  latter  esse  the 
cubic  terns  were  neglected  only  in  less  sensitive  terns  of  eq.  II-5-12. 

In  our  67. 6A  example  at  O'  ■  1A,  the  terns  neglected  in 

eq.  11-5-37  total  about  0.06  of  the  aagnitude  of  the  terns  kept. 

Pram  eg.  XI-5-36  our  steadi-state  solution  for  <  R.  >  is  thus 


<R>  * 


(i-s-41) 


with  I  <^>  I*  obtained  from  eg.  11-5-35. 
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He  now  wish  bo  determine  <  R  >  outside  the  steady-state  regime, 
i.e.  for  the  case  of  a  multilayer  containing  sufficiently  few  layers  that 
a  significant  portion  of  the  incident  bean  is  transmitted  into  the 
substrate. 

The  analysis  for  the  finite-K  regime  is  presented  in  Appendix  10. 
Using  approximations  that  are  similar  in  spirit  to  those  used  in  the 
steady-state  analysis,  we  find  that  the  reflectivity  of  a  multilayer 
containing  J-l  oells  is 


<9  > 


*  ) 


l  T4,p 


Z< 


*  »;p(t-d 
>  e 


fl-5-42) 


where 


Here  <  g  >  is  the  steady-state  solution  of  eq.  II-5-35.  The  subscript 
on  I,  indicates  that  eq.  II-5-15  is  to  be  evaluated  in  the  steady-state. 
In  the  special  case  of  a  bilayer  reflector 
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<*_> 


i  -  “**/  <A*  >  \  ^  *  tit(3-l) 

1  e  /  — - 

*2<Afa>  j  x- 


(a-x-44) 


Our  solution  for  < I  I Z  >  outside  the  steady-state  regime  is 


2<A**>|<p>|4 

do  l*>  =  - — -T 

'j  *.(J-l)D|a 

1-  e  I 


i-  e 


-4(3-1)  ij" 


-  Re 


Zi  (3-1)0  -4(3- 1)1 

e  -  e  J 

i*+  iD/i 


- 413-1)0  -4(3-1) 

e  -  e 

*(‘V-  n 
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Here  is  given  by  eq.  21-5-15. 

<R^>  is  now  given  by  eq.  II-5-36. 
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In  order  to  test  the  o>erall  accuracy  of  the  different  approximations 
used  in  deriving  air  solution,  we  have  written  a  Monte-Carlo  progrss  to 
directly  simulate  the  effect  of  randan  thickness  errors. 

The  reflectivities  of  the  simulated  stacks  are  ocnputed  using  the 
characteristic  matrix  method.  The  exact  (non  x-ray)  characteristic  matrix 
solution  of  eq.  II-1-22  is  used,  rather  than  the  x-ray  approximation  of 
eq.  11-1-14.  A  Gaussian  randan  number  generator  is  used  to  obtain  the 
thickness  errors. 

Pig.  11-5-3  oorpares  our  solution  to  the  Monte  Carlo  results.  The 
tested  values  of  cr  include  the  parabolic  portions  of  the  curves  in  the 
small  -  cr  perturbation  regime,  where  AR/R  ^  1  ,  and  extend  into  the 

large  -  cr  regime  where  the  drop  in  R  fran  its  undegraded  value  is  quite 
substantial. 

By  varying  the  sample  size  and  studying  the  observed  scatter  in  the 
Monte  Carlo  results,  we  estimate  the  accuracy  of  the  results  due  to  the 
finite  sample  size  (usually  2000  simulated  stacks)  to  be  about  1.0005. 

The  theoretical  curves  for  A  A  /  X  in  fig.  II-5-3  were  obtained  by 
numerically  finding  the  half-max  points  of  ^  R,^  .  This  procedure 
implicitly  makes  the  assumption  that  the  half-width  (in  wavelength)  of  the 
expected  reflectivity  is  equal  to  the  expection  value  of  the  half-width  of 
reflectivity. 

This  is  a  reasonable  assumption  to  make  in  order  to  obtain  a  rough 
idea  of  the  behavior  of  AX/^  I  further  one  can  make  argoaents  showing 
explicitly  that  the  assimption  should  be  accurate  if  cr  is  snail.  However 
it  would  be  very  expensive  to  test  this  assimption  using  a  Monte  Carlo 
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Figure  II-5-3 


-  --  (Dashed  Curses) 


simulation;  we  have  not  done  so,  so  the  A  k/  X  curvet  in  fig.  11-5-3 
should  only  be  regarded  as  approximate. 

One  %ould  expect  air  solution  for  <  R,  >  to  becane  less  accurate  as 

the  undegraded  reflectivity  of  the  structure  under  consideration  becomes 

large,  because  in  that  case  neither  inequality  in  eq.  II-5-21  will  hold 

strongly  at  intermediate  values  of  O’  . 

•  • 

The  6A  to  125A  range  covered  by  Henke’s  compilation  of  scattering 
factors  (Henke,  et  al.  1982)  might  be  considered  to  represent  the  full 
extent  of  the  soft  x-ray  regime.  According  to  fig.  II-2-5,  the  largest 
soft  x-ray  reflectivities  that  can  be  attained  with  tungsten-carbon 
multilayers  are  at  short  wavelengths  and  at  wavelengths  just  above  the 
car bon- K  edge. 

Figs.  11-5-4  and  5  compare  our  solution  with  the  results  of  Monte 
Carlo  similations  in  each  of  these  cases.  The  agreement  is  quite  close; 
therefore  our  analysis  can  be  expected  to  accurately  model  the  effect  of 
randan  thickness  errors  in  tungsten-carbon  multilayers  operating  at  any 
wavelength  in  the  soft  x-ray  region,  since  reflectivities  in  the  region 
will  not  be  significantly  higher  than  in  the  two  cases  considered  in 
figs.  11-5-4,5. 

The  6A  case  cannot  easily  be  simulated  at  normal  incidence  because 
the  nsnber  of  layer  pairs  required  would  be  inpractically  large,  making 
the  cost  in  computer  time  prohibitive.  On  the  other  hand,  the  Monte  Carlo 
program  does  not  implement  algorithms  for  off-normal  incidence 
calculations. 
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reflecting  44A  radiation  at  normal  incidence. 
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Figure  11-5-4 
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According  to  eq.  II-2-15,  the  multilayer  parameters  r  ,  p  ,  and  m 
coale  as  see  1  6  .  Thus,  the  effect  of  a  non- normal  angle  of  incidence 

can  be  simulated  by  multiplying  A  nd  ^  by  this  factor.  In  this  way 
fig.  II-5-4  simulates  reflection  of  6A  radiation  from  a  2d«75A  multilayer 
at  an  angle  6  »  arcos (  6A  /  7SA  )  ■  65*. 

The  errors  introduced  with  this  technique  are  of  order  A*  and  f  A, 
and  ao  belong  to  the  class  of  approximations  made  in  sec.  1 1-1.  Since 
these  errors  include  none  of  the  statistical  approximations  made  in  the 
present  section,  the  Monte  Carlo  results  can  be  considered  to  be  a 
reasonable  test  of  the  statistical  assumptions  at  A  »  6A. 

(Of  course,  one  can  simply  regard  fig.  I 1-5-4  as  being  based  on  an 
artificial  example  designed  to  test  our  result  uider  conditions  of  high 
reflectivity) . 
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Part  C)  Discussion  of  Solution 

We  now  discuss  the  general  behavior  of  our  result. 

Frcm  the  approximate  result  of  eq.  11-5-27,  we  aee  that  (to  the 
extent  that  sultiple  reflections  can  be  neglected) ,  the  dephasing  <  A 
has  an  effect  cn  the  reflectivity  that  is  roughly  similar  to  the  effect  of 
/a" ,  i.e.  that  is  similar  to  the  effect  of  absorption. 

In  sec.  XI-2-C  we  saw  that  the  bandwidth  of  an  undegraded  multilayer 
* 

was  of  order  Zm  ,  The  similarity  between  the  effect  of  thickness  errors 
and  an  increase  in  absorption  therefore  causes  the  dependence  of  A  X/x 
on  O’  to  be  approximately  parabolic  (i.e.  linear  in<Af*>),  as  shown  in 
figs.  11-5-3,5. 

In  the  large  C  region  where  extinction  is  due  primarily  to  dephasing 
resulting  from  thickness  errors  (i.e.  in  the  region  <Ag*>  ^  m"  ), 
the  effective  lumber  of  layers  participating  in  the  coherent  reflection 
process  is  of  order  ,  because  absorption  imposes  only  a  minor 

constraint  cn  the  nunber  of  layers  involved.  The  coherent  reflectivity  in 
this  region  should  therefore  have  the  approximate  form: 


!<?>!  —  <A**>  (r  ♦ » p) 


I r+  '  P  I* 


<  A®*> 


(1  -5-46) 


which  is  indeed  the  limiting  behavior  of  eg.  11-5-35. 

Thus,  the  coherent  part  of  the  reflectivity  |<9  ^  I  falls  off 
approximately  as  CT  **  in  the  limit  of  large  errors. 


— T  - 


II-S-38 


In  this  limit  the  incoherent  ocnponent  < |  9  |*>  fells  off  more 
slowly,  as  V  ' 1  (see  fig.  1 1 -5-1) .  in  the  limit  of  very  lerge  errors 
(i.e.  in  the  region  <  A  f*>  >  1  ,  where  our  formalin  is  no  longer 
valid) ,  the  incoherent  reflectivity  can  be  expected  to  enter  a  regime  in 
which  the  reflections  from  the  individual  cells  add  purely  inooherently; 
in  this  case  the  dependence  on  cr  will  be  determined  by  the  effect  of  the 
thickness  errors  on  the  reflectivity  of  individual  unit  cells. 

However,  in  the  regime  of  interest  to  us  (where  <  A  Cl),  the 
incoherent  beam  can  be  considered  to  be  generated  by  "packets"  of  unit 
cells;  the  packets  are  groups  of  unit  cells  which  have  been  displaced  from 
their  expected  phase  positions  by  the  aocvnulation  of  thickness  errors, 
with  the  number  of  unit  cells  within  one  packet  (this  nxnber  is 
— ■  ( <  A  f*>  )  1  )  being  such  that  the  cells  in  the  packet  still  maintain 
strong  coherent  interaction  with  one  another.  (Thus,  by  definition  the 
reflectivity  of  a  packet  is  <  £  > ) . 

Since  (  ?  "  )  ”4  is  the  total  nissber  of  layers  in  the  structure, 

«q.  n-5-40  represents  the  incoherent  sixn  of  the  reflectivities  of  the 
packets;  i.e.  <  |  p  |*>  is  given  by  the  number  of  packets  in  the 
structure  multiplied  by  |  <  9  >  |  *  . 

The  incoherent  reflectivity  will  be  largest  when  there  is  an 
accumulation  of  about  one  radian  of  dephasing  across  the  effective 
thickness  of  the  multilayer,  ao  that  the  aixe  of  a  packet  is  comparable  to 
the  effective  multilayer  thickness;  this  occurs  when 


(1-5-47) 
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If  O’  increases  beyond  this  point,  the  incoherent  reflectivity  begins 
•lowly  to  decrease,  and  the  coherent  reflectivity  begins  a  very  rapid 
decrease. 

Therefore,  using  eq.  II-3-15,  we  can  establish  the  following 
tolerance  on  thickness  errors  that  the  multilayer  »ist  satisfy  in  order  to 
have  high  overall  reflectivity: 

<&$*>  <  =  4  *  30  *  (2d  .  )*  (l-S-4f) 


or 


°(A)  ~  3  '  *  (2dU))X  (t-S-41) 

Thus,  there  is  a  variation  in  sensitivity  to  thickness  errors 
amounting  to  approximately  two  orders  of  magnitude  across  the  soft  x-ray 
region. 

From  eq.  11-5-43  above,  we  see  that  this  same  criterion  can  be  used 
to  characterize  the  effect  of  thickness  errors  on  multilayers  operating 
outside  the  steady-state  regime,  except  that  the  fractional  fall-off  of 
the  reflectivity  with  (T  is  slower  when  J  <  (  M  )  i  this  is  because 

the  steady-state  regime  is  reached  more  rapidly  with  large  CT  due  to  an 
increase  in  the  imaginary  part  of  the  permneter  D  defined  in  eq.  11-5-43. 
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D*  may  be  regarded  as  a  generalised  extinction  par  an*  ter  that 
incorporates  the  effect  of  absorption,  dephasing  due  to  detuning  fren  the 
Bragg  condition,  and  dephasing  due  to  thickness  errors. 
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Section  II-6:  Non-Aocunulating  Random  Thickness  Errors 

Part  A)  Introduction 

In  section  II-5  we  have  analyzed  reflection  from  nominally  periodic 
multilayers  that  contain  accunulating  thickness  errors.  In  such 
multilayers,  the  target  thicknesses  for  the  layers  are  constant, 
pre-determined  quantities,  so  any  errors  made  in  meeting  these  target 
thicknesses  will  accumulate  randomly,  resulting  in  a  degradation  of 
reflectivity  through  the  accunulation  of  dephasing. 

On  the  other  hand,  when  x-ray  multilayers  are  fabricated  with  ISRM, 
the  target  thickness  for  each  high  index  layer  becomes  that  thickness 
which  will  bring  reflection  from  the  layer  into  phase  with  reflection  frcm 
the  preceeding  stack;  thereby  eliminating  the  accunulation  of  dephasing. 

Spiller  et  al.  (I960)  have  verified  experimentally  that  ISRM 
substantially  increases  the  nunber  rtf  layers  that  can  successfully  be 
fabricated  in  an  x-ray  multilayer.  While  the  tolerance  an  thickness 
errors  is  considerably  less  severe  when  ISRM  is  introduced,  it  is  still 
quite  stringent  (on  an  absolute  scale)  compared  to  the  tolerances  which 
govern  optical  sultilayers. 

In  this  section  Me  analyze  the  effect  of  non-accunulating  thickness 
errors.  We  will  find  that  with  this  kind  of  error,  the  reflectivity  loss 
is  due  to  a  degradation  in  the  structure  of  individual  unit  cells. 

Because  of  this  degradation  the  single-oell  reflectivity  is  reduced  by  a 
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Debye-Waller  factor;  the  effect  of  the  errors  will  also  turn  out  to  be 
similar  to  that  of  interlayer  diffusion. 

In  many  ways  the  analysis  of  non-acomilating  errors  is  a  more 
difficult  problem  than  that  of  accumulating  thickness  errors. 

First,  multilayers  fabricated  with  ISRM  are  not  periodic.  Second, 
with  ISRM  there  is  a  complicated  interaction  between  a  hunan  operator  and 
the  multilayer's  structure  and  reflectivity.  Third,  the  wavelength  and 
angle  at  which  the  coating  is  used  may  not  be  the  sene  as  the  wavelength 
and  angle  at  which  the  ISRM  is  performed  (see  sec.  II— 4) . 

Fourth,  the  error  magnitudes  of  greatest  interest  in  the  present  case 
are  larger  than  those  of  interest  in  the  case  of  accumulating  errors.  The 
a  priori  asserptior  of  a  Gaussian  distribution  is  therefore  less 
reasonable  with  non-acamulating  errors,  since  the  errors  may  not  be  small 
ccnpared  to  the  layer  thicknesses.  (Tor  example,  when  the  errors  are 
large  a  Gaussian  distribution  will  imply  the  existence  of  layers  with 
negative  thickness.)  Further,  in  contrast  to  the  accumulating  case,  in  the 
present  oase  the  reflectivity  loss  is  sensitive  to  the  detailed  shape  of 
the  probability  distribution. 

Finally,  when  ISRM  is  used  we  cannot  make  the  usual  kind  of 
siperpositian  to  treat  reflection  of  unpolarised  radiation;  the  radiation 
can  only  be  treated  as  polarised  under  certain  special  conditions. 

It  will  only  be  possible  to  arrive  at  fully  analytic  solutions  to  the 
problem  of  ncn-aecunulating  errors  under  certain  special  conditions. 
Nonetheless,  the  complexity  of  the  nwerical  computations  that  must  be 
made  in  the  general  case  will  be  reduoed  to  a  level  similar  to  that 
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Involved  in  carrying  out  the  two- ty- two  optimization  procedure  described 
in  sec.  11-4  (eqs.  II-1-20,  4-2,  and  4-8) . 

In  particular,  the  analysis  will  eliminate  the  need  for  expensive 
Monte  Carlo  simulations,  which  in  the  present  case  are  particularly 
carplicated  for  the  reasons  listed  above.  In  the  steady-state,  the 
required  computation  reduces  to  the  solution  of  at  most  two  coupled 
algebraic  equations,  which  can  (to  a  fairly  good  approximation)  be  solved 
analytically. 

He  will  use  two  methods  to  analyze  the  problem  of  non-accumulating 
errors.  The  first  method  will  be  mathematically  quite  rigorous  as  long  as 
certain  special  conditions  obtain.  The  second  method  will  be  more 
phenomenological,  but  can  be  expected  to  be  applicable  under  a  broader 
range  of  circumstances  than  the  first  method. 

Despite  the  oanplexity  of  the  problem,  we  will  in  each  case  try  to 
keep  the  assumptions  as  sinple  as  possible.  The  purpose  of  our  analysis 
is  to  explore  the  kinds  of  behavior  that  may  be  encountered  in  practice, 
rather  than  to  make  strongly  quantitative  predictions;  detailed  modeling 
will  be  warranted  as  experimental  work  in  the  field  becomes  more 
extensive. 

Bach  of  our  two  methods  will  be  based  on  a  stronger  version  of  the 
key  ^proximation  used  in  the  accumulating  ease,  namely  that  the  parameter 
^  defined  in  eq.  11-5-10  is  mall  compared  to  one.  In  the 
non-accumulating  case  we  will  show  below  that  it  is  permissible  to  neglect 
quadratic  as  well  as  cubic  powers  of  £  . 
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Since  ^  <  <  1  ,  it  will  be  possible  to  set 


1  *  <?> 


(1-4-1) 


(Clearly  this  is  not  a  good  approximation  for  the  first  few  layers.)  We 
will  find  that  in  this  one  feature  the  problem  of  non-accumulating 
thickness  errors  1b  simpler  to  analyte  than  that  of  accumulating  thickness 
errors. 

Both  methods  of  analysis  will  assure  for  simplicity  that  the 
thickness  errors  obey  Gaussian  statistics  that  are  independent  of  K.  The 
analysis  can  without  significant  modification  be  applied  to  sore 
carpi ica ted  distributions.  Such  distributions  might,  for  example,  keep 
the  layer  thicknesses  always  positive,  or  make  the  FMS  errors  inversely 
proportional  to  the  ISRM  signal  modulation. 

Both  methods  of  analysis  will  assume  that  a  set  of  target 
parameters  have  been  obtained  in  same  way  (for  example,  by  one  of  the 
methods  discussed  in  sec.  1 1-4) .  We  note  that  even  though  p^  is 
nominally  the  phase  thickness  of  the  high  index  layer,  pK  can,  as  shown 
in  eg.  ZI-2-43,  be  regarded  merely  as  a  division  parmeter.  As  such,  it 
can  be  considered  to  determine  the  desired  thickness  of  the  low  index 
layer,  rather  than  that  of  the  high  index  layer,  via  the  relation 
p  *  *  it  -  P*  *  *  it  -  p*  .  The  reflectivity  will  be  much  less 
sensitive  to  the  thickness  ratio  than  to  the  total  cell 
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thickness  Jt  -  •  . 

the  final  value  foe  each  cell  will,  in  the  absence  of  errors  in 

the  Kth  high  index  layer,  automatically  satisfy  eq.  11-4-2  due  to  the  use 
of  the  1SRM  technique;  this  will  be  true  independent  of  any  errors  in  the 
preceeding  layers. 

The  first,  more  rigorous  analysis  will  make  the  following  two 
assumptions  in  addition  to  those  made  above. 

First,  it  trill  assune  that  errors  occur  only  when  fabricating  the  low 
index  layers  (L  layers) .  This  is  reasonable  because  in  many  cases  the 
oscillations  that  the  ISHM  signal  will  undergo  when  the  L  layers  are  being 
deposited  will  be  considerably  weaker  than  the  oscillations  that  occur 
when  the  high  index  layers  {H  layers)  are  being  deposited;  this  will  in 
turn  lead  to  smaller  errors  during  the  B  depositions.  The  L  layer  signal 
oscillations  will  be  particularly  weak  when  multilayers  having  2d-spacings 
considerably  shorter  than  the  carbon-K  edge  are  fabricated,  and  it  is  in 
such  short  spacing  multilayers  that  thickness  errors  have  their  greatest 
effect.  However,  we  should  note  that  in  many  cases  one  might  choose  the 
ISRM  wavelength  to  be  one  which  causes  a  strong  contrast  during  L 
depositions  (Spiller ,  1982b) . 

Another  factor  that  will  tend  to  reduce  the  magnitude  of  the  errors 
made  in  fabricatirig  R  layers  is  that  their  depositions  will  be  truncated 
at  reflectance  maxima;  these  are  ocmparatively  well-defined  positions  on 
the  XSBM  trace. 
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Depending  on  the  algorithm  used  to  detemine  the  p*  .  the  desired  L 
truncstions  nay  not  be  u  easy  for  the  operator  to  relate  to  the  XSIH 
trace,  resulting  in  larger  errors  daring  L  depositions. 

In  addition  to  assuming  that  the  high  index  errors  are  negligible, 
the  first  Method  of  analysis  will  assvsne  that  the  ISFM  probe  bean  can  be 
treated  as  polarized.  This  will  be  strictly  valid  only  if  the  probe  angle 
is  near  grazing  incidence,  near  45*  incidence,  or  near  normal  incidence. 
(As  a  practical  Matter  it  is  doubtful  that  the  probe  6  can  be  made  less 
than  15*  or  so) . 

The  second  Method  of  analysis  will  not  need  to  assume  either  a 
polarized  probe  beam  or  to  assume  negligable  errors  in  the  R  layers;  it 
will,  however,  be  less  accurate  and  less  rigorous  than  the  first  Method. 
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Part  B)  First  Method  of  analysis 

If  there  are  no  errors  in  the  B  layers,  we  can  write  the  phase 
thickness  of  the  Kth  low  index  layer  as 


where  £  £  Zk/X  ccsB  f  ,  and  the  thickness  error  •f  is  a  render, 
variable  with  mean  CT  that  is  assumed  to  obey  Gaussian  statistics  that  are 
independent  of  K.  @L  „  0  *s  *he  desired  phase  thickness  of  the  Kth  low 
index  layer. 

We  note  that  eg.  11-6-2  will  not  obtain  if  there  are  errors  in  the  H 
layers,  since  an  error  in  the  K-lst  B  layer  will  cause  a  t,  term  with 
subscript  K-l  to  be  present  in  *  ,  due  to  the  automatic  error 

condensation  achieved  with  ISHM. 

In  the  remainder  of  this  section,  we  will  for  simplicity's  sake 
sur press  the  K  subscript  cn  all  multilayer  parameters. 

At  this  point  we  also  introduce  the  following  quantities: 

A  *  A  stc  9  A  s  A  sec2© 

MM  _  »•  l 

q  *  A  -  A 

_  “  ml  _ 

n  *  i  ♦  A  n  6  i  ♦  A 

N  *  t  l 
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On  the  grounds  that  £  ie  anall,  we  now  eet  ^  ■  <^>  in 
eq.  II -4-2  (and  ao  we  alao  aet 

as  # 

<  *  >  S  <Re**>*  !<*>!*  <e"*M>  andVf*V<<>  ). 
Under  the  assimption  of  no  errors  in  the  H  layers,  will  still  be 
given  by  eg.  II-4-2  despite  the  presence  of  randcn  errors  in  the  L  layers. 
Bq.  II-4-2  will  now  depend  on  these  errors  only  through  the  first  m' 
term  and  through  the  u  term  in  the  pa  ram  ter  t L  • 

Since  the  distinction  between  R  and  is  a  higher  order  effect,  the 
substitution 


71“  |  <  p  >  |*<  e'ZM  >  =  |<f  >  I*  e‘Z<“  >  (1-6-4) 

where  <  ju>  £  tZ  *  <[  ,  will  not  effect  our  result  except  to 

order  f  •  A  •  %  • 

Then  from  eg.  II-4-2  we  can  aet 


M  f 

?  -  <f>  ♦  *  0(f  A-5)  (2-6-5) 


where 


<*>  -  -7 


f 


< u>  ♦  V  ”  Arctan 

\l-|<»>l  e 


l*|<p>|l«'2“" 

2  - Z<M’>  ttLnV* 


(cont.rtutd  oh  neat  fojt) 
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k?>i  e _ vj, 

4  I  ^  _  ,*  *♦<*">  ,»  '10t'> 

e  -2|<f>l  e  eci2V 

(*- i- 1 ,  ton tmu*d) 

We  note  that  if  f  is  the  the  defect-free  f  value  for  the  Kth  cell, 

<9>  *  90  • 

<  ^  >  must  be  scaled  from  the  wvelength  and  angle  of  the  coating's 
fabrication  to  the  wvelength  and  angle  of  the  coating's  application,  if 
the  two  are  different  (i.e.  one  must  set 
<  s  <  ^  cos  ex  )  /  (  COS  6X  )  ) .  It  will 

then  also  be  necessary  to  propagate  <  9  >  using  eq.  11-1-20  in  two 
separate  sequences;  the  first  under  the  fabrication  conditions  and  the 
second  under  the  application  conditions,  with  the  scaled  values  of  <  9  > 
used  in  the  second  sequence.  Eq.  11-6-5  will  then  still  hold  for  the 
second  sequence,  but  the  error  term  of  order  •  A • £  will  be  different. 
Next  tie  determine  <^>and<t>.  We  have 

■  *  -  »  -  A. 

■  (*-<»>*  P,  .)-  i  -  i  f 

*  <\>  -i  -  i  7 


/ 


-  arcsin 


lP(e)l 


n-i-  7) 
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Similarly 


(  *  "  *  A 


-  n 


<  A 


■s,  x  _  i.  - '  ■ , 


8  <t  > 


T  3 


fa'*-*) 


where  again  both  <  t  >  and  <  ^  >  are  $  dependent. 

In  hf*wndix  9  we  show  that  the  propagation  equation  for  ^ 
(eq.  11-1-20)  has  the  form 


-lit. 


«♦* 


-•e  "  P(B)(  e  ^  ♦  A  sinp^e'^*) 

”»e  *  P  (0 )  {  A^ain/5^  e  4  n^e  '*")  f 

(JL-t-V 

in  the  bilayer  case.  To  derive  eq.  11-6-9  we  have  used  the  definitions  in 
eq.  II-1-15  to  recalculate  the  quantities  in  eq.  11-2-15  under  a 
non-oentroeynmetric  geometry  (L  layer  an  substrate  side,  B  layer  on 
incident  side) ,  and  have  retained  terms  of  order  g  •  A  ,  which  are 
neglected  in  eq.  11-2-15. 
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Assvsning  Gaussian  statistics,  we  have  from  eqs.  II-6-6  and  5-6 


=  e 


-2i <t> 


+  oU‘) 


(1-i-iO) 


(Here  we  have  used  the  same  identity  as  is  used  in  eq.  Il-S-8,  which  can 
be  shewn  to  hold  even  when  the  randan  variable  in  the  exponent  is 
multiplied  by  a  complex  coefficient). 

We  can  new  set  the  expectation  value  of  the  first  term  on  the  right 
of  eq.  II-6-9  equal  to 


-Zi  <t  > 
C 


z 

<?> 


In  order  to  calculate  the  expectation  value  of  the  remaining  terms  in 
the  equation,  we  must  determine  the  quantities: 

<sin£  •  e  '*> 

<  sinfi  eiA*  e'S 

Ksinfi*  e  ‘  e  > 

- . 
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Beginning  with  the  £ir»t  quantity  in  eq.  XI-6-12,  we  expand  ain 
into  expo nentials,  and  first  calculate 


(1-4-13) 


Using  eqs.  11-6-2,7,  and  8,  we  can  aet 


K-Pr'  *  <f>.>  -  -  <t>  -  zs  -  {;?■ 


so  that  under  Gaussian  statistics 


(J-6-14) 


-;<t> 
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C  + 


old) 
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Through  similar  steps  the  remaining  quantities  in  eq.  11-6-12  are 
readily  obtained.  Neglecting  quadratic  and  higher  terms  in  ^  ,  and 
neglecting  terms  of  order  A*  ,  we  obtain  for  the  expectation  value  of 
eq.  II-6-9 
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I 


<f  > 


e  <o  > 


P(8) 

2 


♦  A, 

-  ^r2 « ,i<’>  <♦,>  [a.«‘^  (e:<fc>- e',<<‘>) 


(1-4*16) 

Me  can  now  propagate  <p  >  through  the  stack  using  eqs.  11-6-6,7,8 
and  16.  Since  p  1  ,  we  can  set  <R>  *  l<$>  I  *  to  obtain  the 

intensity  reflectance.  This  is  as  far  as  we  can  proceed  towards  a  fully 
analytic  eolation  outside  the  steady-state  regime.  In  the  steady-state 
regime,  eqs.  XI-6-6  and  16  reduoe  to  two  ooypled  algebraic  aquations,  which 
can  be  rapidly  solved  by  an  iterative  procedure  in  which  <  f  >  ,  <  t  >  , 
and  <  ^  >  are  initially  set  equal  to  their  defect-free  values. 


ka 
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Part  C)  Second  Method  of  analysis 

A  comparison  of  eq.  II-6-16  with  eq.  II-6-9  shows  three  ways  in  which 
non-accunulating  errors  Manifest  thwnselves:  first,  in  the  difference  between 
the  expectation  values  of  the  parameters  j ,  p*,  and  t,  and  the  defect-free 
values  cf  these  parameters;  second,  in  the  difference  between  p  and  its 
expectation  value;  and  third,  in  a  direct  way  through  the  presence  of  the 
Debye-Waller  factors  exp{-2<^*>  ). 

In  our  second  analysis  we  will  treat  the  first  of  these  effects  as 
negligible.  In  essence  we  will  assume  that  the  IS PM  technique  keeps  the 
total  phase  thickness  of  each  cell  a  constant;  in  other  words,  we  assure 
that  each  H  layer  thickness  is  made  by  the  IS  KM  technique  to  ocrpensate 
for  any  OPD  error  present  in  the  preceeding  L  layer.  (To  within  order 
A  ,  we  can  also  set  the  optical  path  difference  equal  to  the  physical 
path  difference.  After  the  expectation  value  iB  taken,  the  error  in  such 
a  substitution  will  always  be  found  to  be  of  order  A*  ). 

With  the  second  method  of  analysis  we  can  treat  the  possibility  of 
errors  in  the  B  layers  as  well.  In  order  to  do  ao,  we  go  to  a 
quasi-oentrosynmetric  decomposition  in  which  the  cell  interfaces  lie 
within  the  low  index  layers.  (The  decomposition  reduces  to  that  shown  in 
fig.  11-2-1  in  the  absence  of  errors).  We  define  these  fictitious 
interfaces  to  occur  at  the  same  physical  heights  above  the  substrate  as 
would  the  Midpoints  of  the  low  index  layers  in  the  absence  of  errors;  thus 
the  interfaces  that  divide  the  cells  db  not  change  positions  in  the 
presence  of  errors. 
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The  error*  db  cause  shifts  in  the  the  positions  of  the  physical 
interfaces  that  separate  the  H  and  L  layers.  We  will  seasure  these  shifts 
relative  to  the  defect-free  positions  of  the  physical  interfaces. 

Finally,  we  will  assure  that  these  shifts  of  the  physical  layer 
interfaces  can  be  treated  as  unoorrelated  Gaussian  randan  variables.  This 
is  suggested  by  the  Debye-Waller  factors  in  sq.  IX-6-16,  which,  as  we  will 
see,  represent  degradations  in  the  reflectivity  of  individual  interfaces 
due  to  displacements  whose  magnitudes  have  a  Gaussian  distribution. 
However,  we  should  note  that  this  assumption  is  not  strictly  consistent 
with  the  results  of  part  B.  For  simplicity  we  will  also  neglect  terms  of 
order  9  •  A  in  our  seoond  analysis.  The  algebraic  details  of  the 
derivation  are  given  in  Appendix  11. 

There  we  show  that 


<p  >  = 


■  «♦» 


~z>t 

€  <9  > 


where 


*  Twec. 

<$*>  E  ~  uiB  cr 

,  i  At 
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Her*  a  and  CT  arc  the  standard  deviations  of  the  errors  in  the  L  and 
B  layers,  and  t  ,  A  »  and  d  are  the  defect-free  values  of  the 
parameters  t,  p*  ,  and  .  Except  for  the  non-oentroeynmetrie  geometry 
and  the  non-zero  errors  in  the  B  layers,  our  result  is  equivalent  to 
eq.  XX-6-16  with  the  expectation  values  of  the  parameters  g  ,  t,  and  p 
set  equal  to  their  defect-free  values,  and  with  terns  of  order  g •  A 
neglected. 

we  also  find  that,  to  an  approximation,  <1  p  I*  >  in  the 
steady-state  is  given  by 
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(This  as s uses,  for  purposes  of  illustration,  that  £  m  t  m  t  .  gee 

•>  % 

below).  Aus, 


<  1 9 1*>  —  0  -  A  1  (t-6-zi) 


as  assumed  earlier.  In  our  67.6A  standard  exanple  (see  sec.  XI-5) ,  with 
no  errors  in  the  B  layers,  <  I  9  |*>  reaches  its  largest  value  in  the 
steady-state  regime  with  (7  >  5A;  in  this  regime  <19  I X >  *  0.004. 

The  physical  reason  for  the  small  size  of  <  I  9  |*>  is  discussed  below. 

In  the  second  analysis,  the  polarization  of  the  ISHM  probe  enters 
indirectly  via  the  defect- independent  parameters  f  .  These  parameters 
are  obtained  by  making  a  defect-free  propagation  of  <  9  >  using 
eq.  II-6-9,  with  eq.  II-4-2  (in  S  or  P  polarization)  used  to  determine 

V 

The  ^  parameters  are  polarization  dependent,  and  so  eq.  11-4-2  will 

not  apply  to  an  unpolarized  XSFM  beam  except  wider  special  conditions. 

However,  since  only  the  defect-free  case  need  be  considered  at  this  point, 

it  would  be  computationally  quite  rapid  to  determine  the  f  values  by 

0 

msnerical  optimization. 

As  an  alternative,  since  our  approach  in  this  second  analysis  is 
somewhat  phenomenological,  we  will  oonsider  it  adequate  to  determine  each 
f'  value  using  a  weighted  average, 


W*p»l,0,« 
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with  •  and  •  given  fay  «q.  II-4-2  in  each  of  the  two 
polarizations.  Xn  the  defect-free  calculation  g  and  "ust  be 
propagated  using  eq.  11-6-17  in  separate  sequences,  but  in  each  sequence 
the  9  value  used  bus',  be  that  given  by  eq.  11-6-22. 

If  necessary  the  resulting  m  values  can  then  be  scaled  to  the 
wavelength  and  angle  of  the  coating's  application;  <  f  >  can  then  be 
re-propagated  in  each  polarization,  and  the  final  intensity  reflectances 
averaged. 
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Part  D)  Discussion 

As  in  MCtion  II-5  we  have  used  a  Monte  Carlo  simulation  to  test  the 
accuracy  of  our  theoretical  treatment.  In  order  to  be  able  to  aaopare 
both  methods  of  analysis,  we  have  only  considered  cases  that  satisfy  the 
restricted  conditions  under  which  our  first  method  of  analysis  can  be 
applied  (see  part  A  above) . 

Our  Monte  Carlo  simulation  program  has  no  off-normal  Incidence 
capability,  so  we  have  used  a  normal  incidence  XSKM  geometry  in  our 
simulations. 

As  discussed  in  sec.  II-5-C,  the  non-normal  incidence  case  can  be 
reduced  to  the  normal  incidence  case  fay  using  a  sec*  6  scaling  of  A^ 
and  A  ,  as  long  as  the  ISRM  beam  can  be  treated  as  polarized.  (See 
part  A  above.  The  error  involved  in  assusing  a  sec*  6  scaling  should  not 
involve  the  statistical  kinds  of  approximations  made  in  this  section.) 

However,  we  note  that  in  the  non-accumulating  case,  the  unpolarized 
reflectivity  can  only  be  treated  as  the  mean  of  the  reflectivities  of  the 
orthogonal  polarizations  wider  the  siaple  model  of  part  C. 

We  have  also  chosen  to  have  the  Monte  Carlo  program  prohibit  layers 
with  negative  thickness;  when  the  Gaussian  randan  number  generator  forces 
such  s  thickness,  the  thickness  is  set  to  zero  instead. 

Our  theory  could  model  this  restriction  by  using  the  appropriate 
non-Gauss i an  probability  distribution  for  the  errors,  but  we  have 
intentionally  left  this  discrepancy  unresolved  since  the  discrepancy  is 
somewhat  representative  of  the  kind  of  oversimplification  that  modeling  of 
this  complex  problem  aust  necessarily  have. 
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Pig.  II-6-1  shore  the  expectation  of  R  plotted  over  a  range  of 
values  from  sero  bo  ~  <dj>/  3  .  (Even  at  cr  ■  ,  the 

fraction  of  layers  isvHer  a  Gaussian  distribution  that  will  have  negative 
thickness  is  fairly  mall.)  the  solid  curve  is  calculated  by  the  method  of 
part  B,  the  dashed  curve  by  the  method  of  part  C. 

The  Monte  Carlo  results  in  fig.  II-6-1  have  been  obtained  using  the 
exact,  non  x-ray  characteristic  matrix  solution  of  aq.  11-1-22.  the  B 
layer  truncations  were  determined  by  numerical  solution  for  the 
reflectance  maximum  at  each  K  value,  in  each  simulated  stack. 

Prior  to  constructing  the  simulated  stacks  that  contain  thickness 
errors,  the  program  constructs  a  defect-free  stack  that  is  optimized  in  a 
two-by-two  fashion  (see  sec.  11-4),  with  the  optimization  performed 
numerically.  The  purpose  of  this  initial  optimization  is  to  obtain  the 
defect-free  p*  division  parameters,  which  determine  the  target  L  layer 
truncations.  As  discussed  above,  our  equations  can  be  used  with  a  set  of 
division  parameters  obtained  by  any  method. 

As  discussed  in  part  A,  one  could  easily  justify  adopting  a  more 
complicated  model  for  both  the  Monte  Carlo  simulation  and  the  theoretical 
analysis,  but  in  this  preliminary  w>rk  we  have  used  only  the  minimal  set 
of  assunptions  discussed  above. 

As  the  figure  shows,  the  agreement  between  the  two  methods  of 
analyses  and  the  Monte  Carlo  simulations  is  quite  good,  with  the  first 
analysis  being  senewhat  more  accurate  than  the  second.  Me  estimate  the 
error  due  to  the  finite  sample  size  to  be  about  t  .001.  The  tw> 
theoretical  calculations  disagree  at  cr  *  0  due  to  neglect  of  terms  of 
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EFFECT  OF  NON-ACCUMULATING 
THICKNESS  ERRORS 

(Tungsten-Carbon  Multilayer  made  with  ISRM) 


J  =  75 


- Theoretical  Curvet 

Monte  Carlo  Simulations 

•  J  =  200 
O  J  =  75 
o  4*30 


J  *  30 


RMS  ERROR  IN  LOW  INDEX  LAYER  (A) 

Plot*  are  for  •  W/C  multilayer  reflecting  67.6  A  radiation 
at  normal  Incidence.  d„  =  7.6  A,  de  *  26.5  A 


Figure  II-6-1 
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order  f  •  A  in  the  Method  of  pert  C. 

Both  analyses  depart  somewhat  from  the  simulations  at  large  cr 
values,  perhaps  because  of  the  positive  layer  thickness  constraint  that  we 
have  placed  in  the  simulation  program,  tie  should  mention  at  this  point 
that  at  large  CT  values,  the  second  method  of  analysis  will  not  always 
track  the  more  rigorous  solution  of  eq.  II-€-16  as  closely  as  it  does  in 
our  standard  67.6A  exaitple;  in  fact  it  is  even  possible  in  some  cases  for 
the  second  solution  to  predict  an  increase  in  reflectivity  as  cr 
increases.  This  can  occur  only  when  the  nominal  B  layer  thickness  is 
sufficiently  small  and  cr  sufficiently  large  that,  under  the  single  model 
of  part  C,  the  H  thickness  can  frequently  become  negative.  The  6imple 
model  can  then  also  predict  an  increase  in  single-cell  reflectivity  as  o~ 
increases,  without  an  increase  in  absorption,  leading  to  an  increased 
reflectivity  for  the  stack  as  a  whole. 

In  contrast  to  the  similar  situation  that  occurs  with  the  L  layers, 
it  would  not  be  valid  to  aorrect  for  this  effect  by  going  to  a 
non-Gaussian  distribution,  sinoe  the  effect  is  an  artifact  intrinsic  to 
the  system  of  cell  decomposition  that  we  use  in  our  in  our  second  method 
of  analysis. 

Bowever,  we  have  always  found  the  second  solution  to  be  accurate  in 
the  region  of  snail  and  moderate  O"  values. 


/' 


II-6-23 


We  now  pc  exseed  to  discuss  our  solution  in  physical  tans.  Zn  order 
to  do  so  we  take  cur  less  rigorous  solution  in  sq.  11-6-17,  and  1st  the 
total  number  of  layers  J  (and  therefore  the  relevant  values  of  K)  became 
sufficiently  large  that  the  multilayer  operates  In  the  steady-state  regime 
and  becomes  periodic. 

We  also  assume  somewhat  arbitrarily  for  purposes  of  discussion  that 
the  errors  in  the  high  and  low  index  layers  have  equal  RMS  magnitudes. 

We  then  obtain  from  sq.  11-6-17 


R  S  l<9>|X 


(This  nay  be  regarded  as  an  approximate  analytic  solution  to  the  two 
ooupled  algebraic  equations  that  our  first  method  of  analysis  reduces  to 
in  the  steady-state.  We  also  note  that  it  is  not  necessary  to  assune 
equal  L  and  H  errors  in  order  to  arrive  at  an  analytic  steady-state 
solution.) 

Bj.  n-6-23  is  sinply  a  standard  Darwin-Prins  solution  (see 
sq.  11-2-11)  with  the  single  cell  reflectivity  irff  degraded  by  a 
Debye-Waller  factor. 
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To  an  extent  the  pretence  of  a  Debye-Waller  factor  la  not  surprising. 
In  x-ray  diffraction  such  factors  result  when  the  phase  position  of  a 
scatterer  is  randomly  displaced  according  to  a  Gaussian  distribution 
(usually  due  to  thermal  notion,  see  Janes,  1965,  p.193) .  The  Debye-Waller 
factor  may  be  regarded  either  as  the  Fourier  transform  of  a  Gaussian 
distribution,  or  as  the  expectation  of  a  oonplex  exponential  whose 
amplitude  obeys  Gaussian  statistics. 

The  central  principle  of  x-ray  diffraction  states  that  when 
extinction  can  be  ignored,  the  amplitude  reflected  from  a  saiple  will  be 
the  Fourier  transform  of  the  diffracting  structure.  (This  is  of  course 
the  basis  by  which  x-ray  diffraction  is  used  to  determine  the  6tructjre  of 
quasi -per iodic  objects.) 

In  the  present  case  extinction  and  multiple  reflections  cannot  be 
neglected  excsept  within  individual  unit  cells?  thus  only  the  unit  cell 
reflectivity  in  eq.  II-6-23  is  multiplied  by  a  Debye-Waller  factor. 

However,  there  is  a  more  subtle  feature  present  in  our  result,  in 
that  aq.  II-6-23  has  the  standard  Darwin-Prins  form.  This  implies  that 
the  total  interaction  of  the  different  cells  in  the  multilayer  is  no 
different  from  the  interaction  that  the  corresponding  set  of  "averaged” 
cells  would  have  in  the  absence  of  errors;  this  is  true  despite  the 
presence  of  extinction  and  multiple  reflections,  which  must  he  considered 
in  the  structure  as  a  whole. 

This  point  can  be  clarified  by  considering  the  similarity  between 
non-accumulating  thickness  errors  and  interlayer  diffusion.  We  will 
continue  to  consider  only  the  sisple  ease  where  the  multilayer  has  equal  L 
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and  B  layer  errors  (or  diffusion  profiles) ,  is  periodic,  end  operates  in 
the  steady-state.  (We  should  also  note  that  interlayer  diffusion  has  been 
analyzed  by  Underwood  and  Barbee  (1982)}  their  treatment  is  discussed 
briefly  in  sec.  IV-1). 

According  to  eq.  II-1-15,  the  parameters  r  and  m  are  essentially 
(truncated)  Fourier  transform  components  of  the  structural  profile 
A  (*). 

If  we  model  the  diffusion  process  as  a  convolution  of  the 
sharp- interface  A  ( z )  function  with  same  smoothening  function,  then  r 
and  m  will  be  multiplied  by  the  appropriate  Fourier  transform  components 
of  the  smoothening  function.  (Since  r  and  M  are  actually  truncated 
Fourier  components,  the  proceeding  argument  is  somewhat  oversimplif ied, 
but  in  Appendix  12  we  show  that  the  result  is  still  correct) .  Since  M  is 
the  DC  component  of  the  transform,  it  is  unchanged  by  diffusion  if  the 
diffusion  profile  is  correctly  normalized.  Thus,  neither  interlayer 
diffusion  nor  non-accumulating  errors  will  change  the  (average) 
transmittance  of  the  cell. 

However,  if  the  diffusion  profile  is  the  same  normalized  Gaussian 
function  as  the  Gaussian  probability  distribution  used  in  our  analysis  of 
non-accumulating  errors,  the  diffusion  will  cause  r  to  be  multiplied  by 
the  easw  Debye-Waller  factor  that  appears  in  eq.  XI-6-23. 

Application  of  the  Darwin-Prins  solution  of  eq.  11-2-11  to  the 
diffusion  problem  then  yields  a  solution  identical  to  that  in  eq.  II -6-2 3 
for  win-accumulating  errors.  The  two  solutions  remain  identical  when 
generalized  to  allow  for  unequal  J»G  errors  or  diffusion  profiles  in  the 
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Is-B  and  B~L  interfaces. 

this  identity  of  results  implies  that  the  randomly  varying  unit  cells 
in  the  stack  interact  together  in  the  sene  way  as  do  a  corresponding  set 
of  cells  having  an  averaged  (or  diffused)  structure. 

This  behavior  is  a  consequence  of  the  small  magnitude  of  9  .  This 
can  be  aeen  by  considering  the  step  in  the  analysis  where  we  take  the 
expectation  value  of  eq.  11-6-9  while  neglecting  |T  t  sinoe  ^  is 
neglected  <  $  >  is  propagated  only  by  the  averaged  or  expectation  values 
of  the  cell  structure.  Under  the  simple  model  of  part  C,  the  averaging  is 
represented  by  the  Debye-Waller  factors. 

Thus,  as  long  as  the  multilayer  contains  at  least  a  few  layers,  it 
need  not  really  be  treated  as  a  randan  structure  except  to  the  trivial 
extent  of  diffusing  the  cell  structure.  (This  is  only  true  within  the 
sinple  model  of  part  C,  of  course) .  In  a  sense  the  ISRM  technique  may  be 
said  to  surpress  the  randomness  by  eliminating  the  accunalation  of 
dephasing. 

Physically,  9  is  snail  in  the  non-aocvxnulating  case  because  the 
ocxiponents  reflected  from  the  different  interfaces  are  unoorrelated  with 
one  another,  and  so,  in  fact,  add  together  only  incoherently. 

Since  incoherent  reflection  is  inherently  *  weak  process  in  the  x-ray 
regime,  the  overall  9  component  is  rnnall;  thus  in  eq.  11-6-19  <  1 9 1  > 
is  of  the  order  of  the  intensity  reflectance  of  each  oell  miltiplied  by 
the  total  nutter  of  cells  (rather  than  by  the  square  of  the  total  nunber 
of  cells  as  would  be  the  ease  with  coherent  reflection) . 
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In  contrast ,  in  the  case  of  aocimulating  errors  the  displaoanents  are 
well-correlated  from  cne  interface  to  the  next,  with  the  correlation 
between  displacements  vanishing  only  between  one  'packet*  and  the  next 
(see  sec.  11-5).  The  packet  as  a  whole  will  radiate  in  a  primarily 
incoherent  way;  however  within  the  packet  the  cells  interact  coherently, 
staking  the  overall  incoherent  reflectivity  considerably  larger 
than  in  the  non-accimilating  case. 

In  the  accumulating  ease  <  ^  >  and  ^  are  therefore  strongly 
inter connected,  and  ao  there  is  no  simple,  non-stochastic  analogue  for 
accumulating  randan  thickness  errors;  in  contrast,  as  we  have  seen  there 
is  a  straightforward  analogy  between  non-accumulating  thickness  errors  and 
interlayer  diffusion. 

The  scaling  of  accumulating  and  non- accumulating  errors  with 
2d- spacing  is  also  different.  In  order  for  the  Debye-Waller  factor  to 
depart  significantly  frcm  one 


2 


cose  c 


In  Z 


(H-b-2*) 


or 


C T  <  0.1  t  (2d) 


(a-b-15) 


This  aay  be  compared  to  the  corresponding  result  in  eq.  II-5-62  for 
accunulating  errors.  In  each  formula  a  should  be  regarded  as  the  total 
RMS  error  per  cell,  rather  than  per  layer. 


Chapter  III  Applications  of  X-Ray  Multilayer  Coatings 


Section  212-1  -  General  Considerations 

Part  A)  Introduction 

In  this  chapter  we  discuss  possible  applications  for  x-ray  multilayers 
in  terms  of  the  theory  developed  in  chapter  II.  We  will  concentrate  on  two 
applications  that  are  of  current  interest  at  the  Laboratory  for  Laser 
Energetics  of  the  University  of  Rochester;  however,  we  will  first  discuss 
some  points  of  general  relevance  concerning  the  performance  of  multilayers 
in  the  x-ray  region. 

In  sec.  IIX-2  we  will  then  discuss  the  possibility  of  using 
multilayers  to  form  cavity  mirrors  for  projected  x-ray  lasers  (wavelength 
of  order  100A).  In  sec.  III-3  we  will  discuss  the  use  of  multilayers  in 
short  wavelength  x-ray  microscopy  (wavelength  of  order  1A) . 
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Part  B)  General  Constraints  on  Multilayer  Performance 

The  results  of  secs.  II-2  and  3  ah an  that  the  reflecting  performance 
of  x-ray  multilayers  is  strongly  limited  by  absorption  and  by  the  weak 
interaction  of  the  layer  materials  with  the  incident  bean. 

According  to  figs.  II-2-4  and  5,  absorption  will  limit  the 
reflectivities  of  multilayers  to  at  most  0.8,  and  much  poorer 
reflectivities  will  oe  typical  in  most  of  the  soft  x-ray  spectrvsn  unless 
new  materials  combinations  prove  feasible. 

In  sec.  II-3  we  have  shewn  that  x-ray  multilayers  have  a  very  narrow 
bandwidth  as  a  result  of  the  weak  coupling  between  the  structure  and  the 
radiation.  This  makes  very  difficult  the  development  of  x-ray  optical 
components  that  are  achromatic  or  wi de-field.  (However,  the  spectral 
selectivity  of  x-ray  multilayers  may  be  useful  in  certain  applications,  see 
below) . 

Eq.  II-3-39  suggests  that  centered  optical  systems  of  moderately  high 
aperture  can  be  devised  for  the  long  wavelength  region  around  100A.  The 
strongest  drawback  to  such  systems  may  be  the  limitation  that  absorption 
will  inpose  an  the  total  ninber  of  reflecting  surfaces. 

As  the  wavelength  falls  below  about  100A,  either  the  angle  of 
incidence  to  the  surface  or  the  2d- spacing  spacing  must  be  reduced; 
according  to  eq.  11-3-21  the  aperture  and  field-of-view  of  a  coated  x-ray 
optical  element  must  then  decrease. 

As  the  wavelength  is  decreased  still  further,  the  necessity  of 
maintaining  layer  thicknesses  of  at  least  a  few  angstroms  will  require  that 
the  multilayers  operate  very  far  from  normal  incidence;  according  to 
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eq.  II-3-23  the  largest  possible  fields  of  view  will  then  be  obtained  at 
angles  as  close  as  possible  to  the  surface. 

The  geometrical  aberrations  in  such  highly  decenter ed  systems  will 


then  injpose  additional  constraints  on  field  of  view. 
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So  far  te  have  cnly  diicussed  the  limitations  on  performance  that  are 
inposed  by  materials  properties.  Experimental  observation  shows  that  the 
reflectivities  of  oultilayers  fabricated  with  present  technology  are 
usually  somewhat  lmer  than  those  of  an  ideal  structure,  although 
uncertainty  in  the  optical  oonstants  makes  such  a  conclusion  difficult  to 
establish  unambiguously.  Fairly  good  agreement  with  theory  can  be  obtained 
when  the  theory  shakes  provision  for  defects  in  the  multilayer  structure 
(Haelbich,  Segmuller,  and  Spiller,  1979;  Barbee,  19B2). 

We  now  discuss  published  experimental  results  in  order  to  gauge  their 
implications  for  practical  devices. 

Multilayer  reflectivity  generally  drops  off  rapidly  when  the 
2d-spacing  is  reduced  below  40  or  50A.  The  measured  performance  also 
generally  acmes  closer  to  the  defect-free  value  under  short  wavelength 
1.5A)  illumination  at  glancing  angles,  than  otherwise. 

To  date  there  have  only  been  tentative  attempts  to  compare  the 
observed  degradation  to  theoretical  models,  and  cnly  the  simplest 
degradation  model  has  been  considered.  Haelbich,  Segmuller,  and 
Spiller  (1979) ;  and  Barbee  (1982) ,  have  each  attempted  to  model  their 
experimental  results  with  an  expression  of  the  form 
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Multiplication  of  the  entire  aultilayer  reflectivity  by  such  a 
De by e-Waller  factor  (rather  than  the  individual  interfacial  reflectivity  as 
in  sec.  II-6) ,  has  been  shewn  by  Eastman  (1978)  to  correspond  to  what  he 
calls  the  identical  film  model  of  roughness,  in  which  same  particular 
roughness  profile  (generally  that  of  the  substrate)  is  replicated 
identically  in  all  succeeding  layers. 

The  problem  of  multilayer  roughness  is  potentially  a  very  complicated 
one,  since  there  is  an  enormous  range  of  possible  statistical  magnitudes 
and  correlations  for  the  roughness,  with  varied  correlations  being  possible 
both  within  each  layer,  and  across  the  different  layers.  Further,  the 
analysis  of  roughness  will  in  general  require  the  solution  of  a 
three-dimensional  electromagnetic  problem. 

We  will  describe  our  own  initial  investigation  of  aultilayer  roughness 
in  the  section  of  chapter  XV  that  deals  with  future  wort;  the  single, 
limiting  case  models  that  we  analyze  there  nay  prove  useful  in  planning  the 
kind  of  detailed  experimental  investigation  of  the  reflecting  performance 
of  irtiltilayers  that  mil  be  neoessary  to  characterize  their  structures  in 
detail.  In  chapter  XV  we  briefly  describe  certain  characteristic 
signatures  in  multilayer  reflecting  properties  that  defects  of  different 
kinds  (including  various  kinds  of  roughness)  will  produce.  However,  the 
experimental  results  published  to  date  appear  to  permit  only  limited 
interpretation. 

Barbee  (1982)  has  fit  the  Debye-Waller  factor  of  eq.  III-l-l  to 
several  measuraments  of  peak  reflectivity  made  in  first  order  at  1.54A  with 
multilayers  having  a  variety  of  2d-spacings.  (See  hia  fig.  7.  Note  that 
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in  this  figure,  the  solid  curve  is  an  interpolation  between  experimental 
points,  rather  than  a  theoretical  calculation.)  Barbee  has  used  four  data 
points  in  his  fit  to  eg.  III-l-l,  and  has  obtained  a  final  WB  error  of 
i5%.  (The  best-fit  value  of  O’  was  3. ISA).  Ihe  experimentally  observed 
fall-off  is  slightly  sharper  than  that  predicted  by  eq.  III-l-l.  (We 
should  note  that  Bar  bee's  multilayers  probably  do  not  contain  enough  layers 
to  operate  in  steady-state). 

Unfortunately,  the  sensitivity  of  multilayer  reflectivity  to  any  of 
the  various  degradation  mechanisms  considered  in  chapter  II  will  have  a 
dependence  on  2d- spacing  that  is  similar  to  that  of  eq.  III-l-l. 

For  exanple,  multilayer  reflectivity  In  the  presence  of  either 
interlayer  diffusion,  non-accunulating  thickness  errors,  or  the  kind  of 
roughness  that  we  term  "smoothening  films*  (see  chapter  IV)  will,  in  the 
eteady-state,  be  approximately  given  by  (here  reproduced  from  eq.  11-6-23) 


ft  = 


-*<r> 


T0e 


Jto-r.e  '  *, 


where  /<$*>  “  Zncr  cos  q/\  . 

jjq.  II 1-1-2  has  the  same  similarity  in  shape  to  eq.  III-l-l  that  the 
Darvin-pr ins  solution  of  eq.  II-2-11  bears  to  our  approximate  solution  in 
#q.  xi-3-1.  thus,  it  will  usually  be  difficult  to  experimentally 
distinguish  the  two  kinds  of  structural  defects  by  varying  the  2d-spacing, 
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particularly  with  a  mall  lumber  of  data  points,  and  with  the  Debye-Waller 
coefficient  serving  as  an  adjustable  parameter. 

Other  kinds  of  nultilayer  defects  will  also  resemble  the  dependence  of 
eq.  II-l-l,  at  least  qualitatively. 

For  example,  we  find  from  eqs.  II-5-27  and  3-15  that  accumulating 
errors  will  have  a  dependence  on  2d- spacing  of  the  approximate  form: 


f U6) 


^Ideal 


where  A  is  independent  of  the  spacing. 

This  function  will  have  a  slightly  sharper  fall-off  than  eq.  111-1-1 
in  the  soft  x-ray  regime,  but  will  have  a  slightly  more-rounded  toe. 
Accumulating  errors  will  produce  a  greater  change  in  bandwidth  than  will 
non-accimulating  errors. 

Eqs.  111-1-1,2  and  3  are  compared  in  fig.  III-l-l;  in  each  equation 
the  (7  value  has  been  chosen  to  make  the  degradations  equal  to  0.5  at  a 
cannon  4-spacing  (d„  ).  The  optical  constants  used  in  eq.  II 1-1-2  are 

those  of  timgsten  and  carbon  at  X  *  1.66A  (comparable  to  Barbee's 
wavelength  of  1.54A). 

The  near  coincidence  of  these  curves  shows  the  difficulty  in  trying  to 
use  the  dependence  of  reflectivity  on  2d-spacing  to  determine  the  nature  of 
the  structural  defects  in  suit! layers. 
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SCALING  OF  REFLECTIVITY  DEGRADATION 
WITH  2d-SPACING 
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To  date  there  has  been  one  direct  experimental  indication  that 
roughness  has  •  detectable  effect  on  multilayer  reflectivity;  Spiller  et 
al.  (1980)  have  found  that  the  measured  reflectivity  of  a  55  layer 
gold-palladiim/carbon  multilayer  is  increased  if  the  detector  aperture  is 
increased,  indicating  that  the  reflected  beam  contains  scattered  radiation. 

Me  might  mention  that  at  short  x-ray  wavelengths,  the  onset  of 
reflectivity  I06S  predicted  by  eq.  XII-1-2  will  became  somewhat  sharper 
than  the  onset  at  longer  wavelengths.  The  presence  of  structural  defects 
satisfying  this  equation  might  then  explain  why  moderate  2d- spacing 
multilayers  tend  to  have  short  wavelength  reflectivities  that  are  closer  to 
the  defect-free  values  than  are  the  reflectivities  at  longer  wavelengths. 

Haelbich,  Segmuller,  and  Spiller  (1979)  have  made  use  of  a 
multidimensional  optimization  code  that  numerically  determines  the 
particular  thicknesses  for  each  of  the  layers  in  a  multilayer,  and  the 
particular  coiplex  indicies  of  refraction,  that  will  produce  the  best 
agreement  with  a  measurement  of  reflectivity  versus  0  under  monochromatic 
illumination.  A  Debye-Waller  coefficient  then  forms  a  final  degree  of 
freedom.  The  agreement  acM  *ved  is  moderately  good  (see  their  fig.  3) ,  but 
with  so  many  free  parameters,  it  is  difficult  to  judge  hew  reliable  a  test 
this  is  of  the  assumed  degradation  model. 

In  chapter  TV  we  will  discuss  briefly  the  kinds  of  measurements  with 
which  one  might  characterize  the  structural  defects  in  multilayers.  For 
the  monent  we  draw  one  important  conclusion  from  existing  experimental 
work,  namely  that  x-ray  multilayers  can  be  expected  to  show  a  rapid  loss  in 
reflectivity  when  the  2d- spacing  drops  be lew  about  40A. 
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Sue *  a  degradation  is  not  surprising.  According  to  sq.  II-6-25,  an 
interfacial  gradient  with  a  characteristic  depth  of  order  4A  will  cause  a 
factor  of  two  reflectivity  loss  in  a  multilayer  having  a  2d-spaclng  of  40A. 

One  might  consider  such  a  fall-off  in  reflectivity  to  be  essentially 
inevitable,  since  an  interface  cannot  be  sharp  on  a  sub-atonic  scale. 

While  this  may  prove  to  be  correct  as  a  practical  matter,  we  might 
note  that  there  are  x-ray  crystallographic  argunents  which  suggest  that  the 
ultimate  limit  imposed  by  atonic  structure  on  2d- spacing  is  actually  much 
smaller.  Under  a  classical  atonic  model,  the  gradient  in  the  atonic 
electron  density  distribution  causes  a  fall-off  in  the  atonic  scattering 
factor;  this  fall-off  is  usually  measured  in  terms  of  the  quantity 
sin  6  /  X  (see  James,  1965,  chapter  III,  for  a  discussion  of  these 
concepts) .  The  fall-off  typically  beoones  pronounced  when 
sin  6/X  ~  0.5A'1  (i.e.  2d  -  2A). 
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Section  XIX-2  -  Cavity  Mirror*  for  X-Ray  Laser* 

Part  A)  Multilayer-Coated  Cavity  Mirrors 

There  has  been  considerable  interest  in  recent  years  in  the 
possibility  of  achieving  optical  vilification  at  x-ray  wavelengths 
(Forsyth,  et  al.,  1976;  Maynant  and  Elton,  1976).  Zn  cost  schemes,  the 
anplifying  region  is  a  plama  formed  at  or  near  the  focus  of  a  high  power 
visible  pup  laser  (where  characteristic  tvperatures  are  in  the  range  of  1 
kev).  Such  an  x-ray  amplifier  could  form  the  basis  for  a  true  x-ray  laser 
if  a  cavity  of  adequate  quality  could  be  constructed. 

Much  of  the  x-ray  laser  research  conducted  to  date  at  the  Laboratory 
for  Laser  Energetics  (U£)  has  concentrated  on  a  scheme  based  on  three-body 
recombination  (Bhagavatula  and  Yaakobi,  1978,  Conturie,  1982) . 

In  these  experiments,  a  pulse  from  the  pup  laser  (duration  typically 
1  nanosecond)  is  focussed  through  a  slit  of  width  —  200 Mm  (or  past  a 
half-plane)  onto  a  target  having  a  2  of  order  10  (the  target  is  typically 
located  ~  200 mm  below  the  slit).  The  laser-produced  plasma  undergoes  a 
jet-like  expansion  back  into  the  direction  of  the  incident  beam,  where  it 
is  cooled  through  oollisian  with  a  cold  plasma  formed  in  the  slit. 

It  has  been  established  experimentally  that  population  inversions 
occur  in  the  rapidly  cooled  plasma.  Most  probably  through  the  mechanism  of 
three-body  recombination  (Conturie,  1982) . 

The  most  extensively  studied  inversion  is  one  that  has  been  observed 
on  the  iw4  to  hf3  transition  in  heliisn-like  alroimsn  (wavelength  ~  130A) , 
but  inversion  has  also  been  observed  from  a  nunber  of  lof-to-moderate  Z 
targets,  and  on  a  nunber  of  x-ray  transitions  (including  in  particular  the 
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80A  3-2  transition  in  hydrogen- like  flour ine) . 

For  •  frequency-tripled  Nd:Yag  pump  laser  (X  ■  .35mm)*  Ccnturie 
(1982)  has  tentatively  identified  9  <  Z  ^  13  as  the  region  of  best 
performance,  which  will  tend  to  limit  transition  wavelengths  to  the  region 
above  50a. 

For  this  reason  our  discussion  will  initially  concentrate  on  the 
possibility  of  using  multilayers  to  form  cavity  mirrors  for  the  100A 
regime. 

According  to  fig.  II-2-5,  the  reflectivity  of  a  tungsten-carbon 
multilayer  at  such  wavelengths  will  be  of  order  0.1.  As  shown  in 
fig.  II-2-4,  we  might  hope  to  achieve  higher  reflectivities  with  new 
materials  combinations,  but  it  is  clear  that  in  the  soft  x-ray  region, 
multilayers  cannot  achieve  the  reflectivities  of  virtually  unity  that  are 
possible  in  the  visible. 

However,  even  with  a  reflectivity  of  0.1,  it  is  possible  to  get  above 
the  threshold  for  true  lasing,  and  so  improve  the  coherence  of  the  vnitted 
radiation.  If  the  gain  in  the  ariplifier  can  be  raised  above  10,  so  that 
each  90%  loss  an  reflection  will  be  compensated  for  by  a  pass  through  the 
amplifier  (fig.  111-2-1). 

the  key  to  the  two-mirror  configuration  shown  in  the  figure  is  that 
the  multilayer  structures  have  been  tuned  to  reflect  precisely  at  normal 
incidence  (i.e.  their  structures  are  optimized  according  to  the  results  of 
sec.  11-2-B) ,  so  that  each  single  lossy  reflection  is  ocnpensated  for  by  a 
pass  through  the  wplifier. 


PROPOSED  CAVITY  CONFIGURATION 


F i gure  III- 
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According  to  iq».  II-3-20  and  1-11,  the  absolute  tolerance  on 
systematic  error  in  the  wit  cell  thickness  d  that  ia  iiqposed  by  the  need 
to  have  adequate  reflectivity  at  normal  Incidence  Is 


d,;> 

T  * 
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If  these  tolerances  are  not  set,  so  that  the  nultilayer  resonance 
angle  is  not  sufficiently  near  6  ■  0*,  one  night  turn  to  a  cavity 
configuration  like  that  shown  in  fig.  111-2-2*,  where  a  grazing  reflection 
is  employed  to  return  the  bean  to  the  second  mirror. 

However,  there  are  now  two  lossy,  near  normal- incidence  reflections 
required  for  every  pass  through  the  mrplifier,  so  that  a  higher  gain  is 
required  than  is  needed  with  the  two  mirror  cavity. 

Further,  because  of  the  synmetry  that  exists  between  the  grazing 
regime  and  the  normal  incidence  regime  (see  sec.  XX-3)  r  w  can  show  that 
when  the  aultilayer's  detuning  error  is  sufficient  to  rule  out  the 
configuration  of  fig.  IIX-2-1,  the  grazing  reflection  will  then  likely 
OCTTur  *°°  large  an  angle  for  efficient  reflection. 

According  to  eq.  XI-3-39,  the  half-angular  width  of  the  multilayer 
reflectenoe  profile  at  noraal  incidence  is 


III-2-5 


OTHER  CAVITY  CONFIGURATIONS 

- 


gain  «*  100 


A.  Cavity  formed  from  non-optimal  multilayers 


Ring  should  contai.  at  least 
M  mirrors,  where 


M  = 


ir*  lm  (n) 


2(1  -  Re(n))5/*  In2 


Refractive  Index  of  mirror 
material  Is  n. 

The  limiting  0  is  given  by 


Q  *  exp  [-4tt 


Re  ) 

I  n*  - 1 1 


] 


B.  Ring  cavity  formed  from  gracing  reflectors 


•waning  p  «  1,  «  £»*  .  This  is  also  the  approximate  half-width 

of  the  region  about  nonaal  incidence  within  which  S6,WNM  is 
approximately  constant.  The  critical  angle  (no  longer  a  sharply  defined 
cutoff  at  long  wavelengths)  is  given  by  (e  *  •  where  is  the 

index  of  the  bulk  reflector.  Xn  order  that  the  two-mirror  sehme  of 
fig.  II 1-2-1  be  ruled  out  while  the  scheme  of  fig.  XIX-2-2a  be  still 
worthwhile. 


(EL-2  -4) 

or 

A  L" 

Kl  £  -3T- 

Cven  in  the  soft  *-ray  regime  where  absorption  is  quite  significant, 
it  will  generally  be  possible  to  technically  satisfy  eq.  II 1-2-5. 
McnetheltM,  it  is  clear  that  in  addition  to  implying  a  greater  loss  per 
cavity  transit,  the  incorporation  of  s  grazing  reflection  as  in 
fig.  III-2-2s  will  not  greatly  reduce  the  fabrication  requirement  of 
•q.  II 1-2-1,  because  it  will  not  in  general  be  possible  to  satisfy  the 
inequality  of  eq.  II 1-2-5  in  a  strong  *my. 
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roc  th>  above  reasons  we  consider  the  aost  premising  aevity 
configuration  for  an  a-ray  laser  in  the  100&  regime  to  be  one  based  on 
nonsal  incidence  Multilayers. 


... 
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tart  B)  Xing  Cavity 

The  degradation  mechanists  diaeusaed  in  ate.  XII-1  isply  that  the 
cavity  configuration  of  fig.  XXX-2-1  will  beocne  inadequate  foe  tranaition 
wavelengths  shorter  than  about  40A. 

At  ahort  %«velengtha  one  night  attwnpt  to  oonfigure  the  multilayers  in 
a  polygonal  ring  in  order  to  increase  the  2d-spacing,  but  the  resulting 
exponentiation  of  the  reflection  lots  will  imply  a  low  overall  throughput 
(the  ideal  aingle-pass  resonant  reflectivity  from  multilayers  can  be  at 
most  0.8  in  the  soft  x-ray  regime) . 

Bovver,  as  the  angle  of  incidence  against  a  specular  reflecting 
surface  goes  to  sero,  the  reflectivity  can  approach  unity. 

On  this  basis,  Brener  and  Kaihola  (1980)  have  devised  a  ring  cavity  in 
which  a  large  number  of  grazing  reflections  fren  polished  substrates  are 
used  to  return  the  beam  to  the  wplifier  (fig.  III-2-2b) .  They  have  shown 
that  if  the  ring  contains  sufficiently  many  mirrors,  the  fraction  of  the 
beam  returned  to  the  snplifier  approaches  the  limiting  value 


where  n#  »  1  ♦  is  the  index  of  the  mirror  material. 

In  Appendix  13,  we  show  that  in  order  to  approach  the  limiting 
throughput,  the  ring  must  contain  at  least  M  reflectors,  where 
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(We  note  that  M  ^  ia  often  substantially  larger  than  the  nnber  of 
•irrors  that  teuld  be  required  sinply  to  put  the  incidenoe  angle  at  the 
critical  angle) .  The  fall-off  in  R  as  the  wisher  of  mirrors  ia  decreased 
below  M  ia  very  rapid. 

We  also  ahow  in  Appendix  13  that  a  general  formula  for  the  throughput 
in  5  polarization  tdien  M  )>>  M  ia 


Mi"/-*) 

|  n*-  1 1  m -Z-S) 


Eq.  111-2-8  can  differ  substantially  from  eq.  III-2-6  in  the  ultrasoft 
x-ray  regime. 

Using  eq.  111-2-6,  and  the  limited  set  of  x-ray  optical  constants  that 
were  then  available,  Brener  and  Kaihola  (I960)  identified  two 
wavelength-ref lector  combinations  where  moderately  high  performance  might 
be  expected;  these  are  polystyrene  at  A  *  67. 6A  (Q^  ■  0.04) ,  and 
magnesium  at  X  «  2.1A  (Qs  •  0.3). 

They  reported  that  even  performances  such  as  these  appeared  generally 
to  be  limited  to  isolated  wavelength-material  combinations.  Further,  a 
literature  review  reveals  that  the  particular  optical  constants 
measurements  which  form  the  basis  for  the  67, 6 A  polystyrene  reflector 
have  essentially  been  retracted  by  the  spectroscopist  who  made  the 
measurements  (Lukirskii  et  el.,  1965,  initial  results  reported  in  Lukirskii 
at  al.,  1964). 


In  addition,  application  of  aq.  111-2-7  to  the  2.1A  Mg  reflector 
indicates  that  the  ring  will  require  approximately  240  reflections.  Even 
if  arte  avoided  the  resulting  alignment  problem  by  using  a  cylinder  rather 
than  a  polygon,  the  minim un  allowable  size  of  each  facet  (dictated  either 
by  the  amplifier  length  or  by  damage  considerations)  would  probably  imply  a 
cavity  circumference  greater  than  one  meter.  Thus,  one  would  have  to  add 
to  the  formidable  difficulty  of  achieving  amplification  at  short 
wavelengths,  the  need  to  maintain  the  inversion  for  even  longer  than  the 
nanosecond  durations  that  are  oontenplated  in  the  100A  regime. 

For  these  reasons  we  initially  considered  the  ring  cavity  to  be 
considerably  less  promising  than  the  multilayer  configuration  of 
fig.  111-2-1. 

Bcwever,  we  have  recently  conducted  a  materials  search  using  the 
recent  aonpilatian  of  optical  constants  made  by  Henke  et  al.  (19B2) .  The 
results  are  shown  in  table  II 1-2-1 i  in  the  longer  part  of  the  soft  x-ray 
spectrisn  the  search  program  has  identified  cavity  materials  that  yield 
nominal  throughputs  of  up  to  0.4. 

While  these  are  not  as  high  as  the  reflectivities  that  are  obtained  in 
a  search  far  sultilayer  materials  (fig.  II-2-4),  they  are  certainly  high 
•nough  to  make  the  ring  cavity  look  quite  interesting  for  transition 
wavelengths  above  —  65A. 
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Optimum  Materials  for  Ring  Cavity 


0 


X  (A) 

(S  Pol) 

M(l/2) 

M(Crlt) 

Element 

124.000 

.3641025953573 

7.51 

7.29 

RHODIUM 

121.040 

.3921299229287 

7.47 

7.54 

RHODIUM 

11B.1S1 

.4182819668938 

7.44 

7.80 

RHODIUM 

115.331 

.3955354172620 

7.98 

8.11 

RHODIUM 

112.578 

.3S6659S935351 

8.71 

8.40 

RHODIUM 

109.891 

.3291879913218 

9.32 

8.67 

RHODIUM 

107.268 

.2850703515461 

10.24 

8.96 

RHODIUM 

104.707 

.2522392515925 

10.83 

9.05 

RUTHENIUM 

102.208 

.2264432242987 

11.57 

9.31 

RUTHENIUM 

99.768 

.2023748964901 

12.34 

9.58 

RUTHENIUM 

97.387 

.2029227829216 

12.62 

9.80 

SILVER 

95.062 

.2125785074919 

12.94 

10.22 

SILVER 

92.793 

.1775126130385 

14.24 

10.64 

SILVER 

90.578 

.1483135039153 

15.53 

11.05 

SILVER 

•8.416 

.1163074612112 

17.19 

11.51 

SILVER 

•6.306  . 

.0890465894562 

17.88 

11.28 

RUTHEN- 

84.246 

.0740252495011 

19.06 

11.59 

RUTHEN 

82.235 

.0629664087553 

20.16 

11.90 

RUTHEN 

80.272 

.0540046097159 

24.38 

14.03 

INDIUM 

78.356 

.0584580751483 

25.19 

14.71 

INDIUM 

76.485 

.0628288060968 

26.04 

15.42 

INDIUM 

74.660 

.0444836744491 

28.89 

16.13 

INDIUM 

72.878 

.0299965650159 

32.04 

16.85 

INDIUM 

71.138 

.0353809877597 

30.37 

16.36 

LANTHUNUM 

69.440 

.0449242525942 

30.39 

17.00 

LANTHUNUM 

67.782 

.0558698077888 

30.42 

17.66 

LANTH’JNUM 

66.165 

.0717400426828 

30.31 

18.42 

LANTHUNUM 

64.585 

,0750066570867 

31.28 

19.18 

LANTHUNUM 

63.044 

.0416450779567 

36.13 

19.99 

LANTHUNUM 

61.539 

.0219788948089 

41.04 

20.71 

LANTHUNUM 

60.070 

.0140820920987 

57.01 

27.26 

CARBON 

58.636 

.0133696653412 

59.19 

28.14 

CARBON- 

57.236 

.0124442741563 

61.66 

29.08 

CARBON 

55.870 

.0113654363679 

64.43 

30.08 

CARBON 

54.537 

.0101470675794 

67.71 

31.22 

CARBON 

53.23S 

.0087420769778 

71.40 

32.40 

CARBON 

51.964 

.0070910928000 

76.14 

33.81 

CARBON 

50.724 

.0053461688895 

•2.02 

35.42 

CARBON 

49.513 

.0035967330555 

•9.62 

37.32 

CARBON 

48.331 

.0018793132260 

100.91 

39.78 

CARBON 

47.177 

.0006449053022 

118.27 

43.09 

CARBON 

46.051 

.0000561201248 

154.81 

48.82 

CARBON 

44.952 

.0000249689810 

72.20 

21.49 

CHROMIUM 

43.879 

.0000234809186 

74.10 

22.01 

CHROMIUM 

42.832 

.0000226497750 

75.92 

22.53 

CHROMIUM 

41.809 

.0000213902624 

77.90 

23.08 

CHROMIUM 

40.811 

.0000199946934 

79.99 

23.64 

CHROMIUM 

39.837 

.0000185326460 

•2.25 

24.24 

CHROMIUM 

38.886 

.0000173651901 

•  4.41 

24.82 

CHROMIUM 

37.958 

.0000160923393 

•6.77 

25.44 

CHROMIUM 

37.052 

.0000147284590 

•  9.27 

26.08 

CHROMIUM 

36.168 

.0000133237546 

91.91 

26.75 

CHROMIUM 

35.304 

.0000121280221 

94.63 

27.44 

CHROMIUM 

Table  III-?-] 
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Optimum  Materials  for  Ring  Cavity 


Q 


x  (A) 

(S  Pol) 

Hll/2) 

H(Crlt) 

Element 

34.462 

.0000106690404 

97.59 

26.15 

CHROMIUM 

33.639 

.0000094204381 

100.67 

28.90 

CHROMIUM 

32.536 

.0000081844422 

103.95 

29.67 

CHROMIUM 

32.052 

.0000069709225 

107.46 

30.48 

CHROMIUM 

31.287 

.0000058342894 

111.31 

31.35 

CHROMIUM 

30.540 

.0000051698340 

114.91 

32.22 

CHROMIUM 

29.811 

.0000043060991 

119.01 

33.13 

CHROMIUM 

29.100 

.0000035484615 

123.49 

34.12 

CHROMIUM 

28.405 

.0000028021737 

128.37 

35.15 

CHROMIUM 

27.727 

.0000022382416 

111.57 

30.14 

COBALT 

27.065 

.0000021347739 

114.56 

30.92 

COBALT 

26.419 

.0000020114377 

117.72 

31.72 

COBALT 

25.789 

.0000018689708 

121.06 

32.55 

COBALT 

25.173 

.0000021463589 

209.34 

57.12 

BERYLLIUM 

24.572 

.0000026708443 

212.73 

58.56 

BERYLLIUM 

23.986 

.0000032707119 

216.28 

60.03 

BERYLLIUM 

23.413 

.0000041001999 

219.67 

61.54 

BERYLLIUM 

22.854 

.0000050253768 

223.28 

63.09 

BERYLLIUM 

22.309 

.0000060673091 

227.08 

64.68 

BERYLLIUM 

21.776 

.0000077398673 

230.37 

66.31 

BERYLLIUM 

21.257 

.0000096400645 

233.91 

67.98 

BERYLLIUM 

20.749 

.0000117805326 

237.66 

69.69 

BERYLLIUM 

20.254 

.0000146593254 

241.24 

71.45 

BERYLLIUM 

19.770 

.0000179611181 

245.01 

73.25 

BERYLLIUM 

19.299 

.0000220340768 

248.78 

75.09 

BERYLLIUM 

18.838 

.0000284884506 

251.93 

76.98 

BERYLLIUM 

18. 388 

.0000361553334 

255.28 

78.92 

BERYLLIUM 

17.949 

.0000418610502 

259.80 

80.91 

BERYLLIUM 

17.521 

.0000486134398 

264.33 

82.94 

BERYLLIUM 

17.103 

.0000606588480 

267.92 

85.03 

BERYLLIUM 

16.694 

.0000727212247 

272.05 

87.17 

BERYLLIUM 

16.296 

.0000881710569 

276.03 

69.36 

BERYLLIUM 

15.907 

.0001078347626 

279.88 

91.61 

BERYLLIUM 

15.527 

.0001306676291 

283.85 

93.91 

BERYLLIUM 

15.157 

.0001593946264 

287.71 

>6.26 

BERYLLIUM 

14.795 

.0001917968792 

291.78 

9J.6B 

BERYLLIUM 

14.442 

.0002272059188 

296.11 

101.15 

BERYLLIUM 

14.097 

.0002694053500 

300.41 

103.69 

BERYLLIUM 

13.760 

.0003208726122 

304.64 

106.29 

BERYLLIUM 

13.432 

.0003774263806 

309.10 

108.95 

BERYLLIUM 

13. Ill 

.0004432890668 

313.57 

111.68 

BERYLLIUM 

12.798 

.0005249066838 

317.87 

114.47 

BERYLLIUM 

12.493 

.0006232975746 

322.08 

117.34 

BERYLLIUM 

12.195 

.0007324482448 

326.50 

120.27 

BERYLLIUM 

11.904 

.0008701721502 

330.63 

123.28 

BERYLLIUM 

11.620 

.0010350574823 

334.68 

126.36 

BERYLLIUM 

11.342 

.0012161348034 

338.98 

129.51 

BERYLLIUM 

11.071 

.0014132805211 

343.52 

132.74 

BERYLLIUM 

10.807 

.0016628845464 

347.69 

136.06 

BERYLLIUM 

10.549 

.0019342937161 

352.12 

139.45 

BERYLLIUM 

10.297 

.0022503197159 

356.49 

142.92 

BERYLLIUM 

10.052 

.0025908233891 

361  12 

146.48 

BERYLLIUM 

9.812 

.0029B04056776 

365.73 

150.13 

BERYLLIUM 

9.577 

.0034261007700 

370.30 

153.86 

BERYLLIUM 

Table  111-2-1 
(conti nued ) 
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Optimum  Materials  for  Ring  Cavity 


Q 


X  (») 

(S  9o1) 

«(l/2) 

K(CMt) 

IlcKCnt 

9.349 

.0039337953130 

374.87 

157.69 

BERYLLIUM 

9.136 

.0044735839306 

379.69 

161.61 

BERYLLIUM 

9. 908 

.0051194666344 

384.24 

165.63 

BERYLLIUM 

8.695 

.0058053048374 

389.06 

169.74 

BERYLLIUM 

8.488 

.0065369399343 

394.14 

173.95 

BERYLLIUM 

8.385 

.0073837789903 

398.95 

178.27 

BERYLLIUM 

8.087 

.0083597341518 

403.63 

182.69 

BERYLLIUM 

7.894 

.0093833485641 

408.61 

187.22 

BERYLLIUM 

7.706 

.0105731968466 

413.35 

191.86 

BERYLLIUM 

7.533 

.0118135474950 

418.38 

196.61 

BERYLLIUM 

7.343 

.013)696703343 

423.46 

201.48 

BERYLLIUM 

7.167 

.0147163189309 

428.35 

206.47 

BERYLLIUM 

6.996 

.0163986744697 

433.27 

211.57 

BERYLLIUM 

6.839 

.0183406833019 

438.19 

216.80 

BERYLLIUM 

6.666 

.0303318331051 

442.93 

222.16 

BERYLLIUM 

6.507 

.0334733434753 

447.97 

227.65 

BERYLLIUM 

6.353 

.0348780019589 

452.85 

233.28 

BERYLLIUM 

6.300 

.0374467481336 

457.83 

239.04 

BERYLLIUM 

Table  III-2-1 
(continued ) 
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Part  C)  Daaage 

Itie  high  tenperature  environment  of  the  plasma  amplifier  suggests  that 
damage  may  be  a  significant  problsn  foe  x-ray  laser  cavities. 

While  an  accurate  assessment  of  the  damage  threshold  can  only  be  made 
experimentally,  there  is  reason  to  believe  that  the  limitation  isposed  by 
damage  may  be  less  severe  than  one  might  expect,  at  leaat  in  the  context  of 
the  LLE  recombination  laser  experiments  described  in  part  A. 

In  the  first  place,  in  these  experiments  amplification  does  not  occur 
precisely  in  the  focal  region.  If  fig.  XXX-2-1  is  oriented  with  the 
mirrors  at  left  and  right,  and  if  the  pimp  beam  is  brought  from  the  top  to 
a  left-right  line  focus  an  e  target  at  bottom,  then  the  region  of  highest 
time-averaged  gain  will  be  formed  near  the  cooling  slit,  whose  axis  will  be 
parallel  to  the  optical  axis.  The  cooling  slit  will,  in  turn,  be 
positioned  —  200  above  the  target.  The  time-averaged  high  gain  region 

will  have  a  width  of  ~  200  and  a  length  approximately  equal  to  that  of 

the  line  focus.  (This  time-averaged  description  is  somewhat 
oversimplified;  Oonturie's  simulations  (1962)  show  that  the  region  of  high 
instantaneous  gain  travels  outward  with  the  velocity  of  the  expanding 
plasma. ) 

Thus,  the  mirrors  could  be  shielded  from  the  hot  plasma  formed  at  the 
focal  line,  and  need  aee  only  the  relatively  cold  plasma  in  the  amplifying 
region  near  the  slit.  Oonturie's  experiments  shew  that  the  particulate 
blow-off  would  not  flow  isotropically  through  the  slit,  but  vould  instead 
flew  preferentially  in  the  direction  of  the  pimp  laser  (Conturie,  1982) . 
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The  x-ray*  emitted  by  the  amplifier  art  then  likely  to  be  the 
principal  agent  of  damage  in  the  mirrors.  Even  though  Boat  of  the  x-ray 
radiation  emitted  by  the  plama  may  be  flourescent  (arising  from  a  number 
of  x-ray  transitions) ,  it  is  the  stimulated  •nission  which  is  most  likely 
to  cause  damage,  because  it  will  be  concentrated  on  a  narrow  area  of  the 
mirrors.  Stimulated  amission  into  the  coherent  spatial  modes  of  the 
flat-mirror  cavity  of  fig.  XI1-2-1  will  undergo  virtually  no  diffraction 
spreading  at  the  mirrors,  because  the  cavity  length  can  be  no  more  than  a 
few  centimeters  if  ooherent  oscillation  is  to  be  established  during  the 
—1  nsec  duration  of  the  gain,  (This  assimes ,  somewhat  arbitrarily  for 
purposes  of  illustration,  that  the  gain  duration  is  of  the  order  of  the 
purrp  pulse  duration) . 

Por  purposes  of  a  rough  estimation  we  will  consider  the  maximixn 
permissible  beam  energy  to  be  determined  by  thermal  loading  on  the  mirrors. 
Tungsten  has  a  thermal  conductivity  K  »  1.78  Matt/cm  -'K  and  a  heat 
capacity  C  «  0.132  Joule/gm  -*K  (Bechtel,  1975).  The  thermal  properties  of 
graphite  are  typically  K  ■  1.29  Matt/on  —  *K,  C  •  0.712  Joule/ -*K  (CRC 
Handbook  of  Tables  for  Applied  and  Engineering  Science,  1973). 

The  distanoe  that  heat  can  propagate  in  a  time  t  is  of  order 


5  — 


(n-i-  <?) 


which  (assaying  bulk  densities)  works  out  to  about  3000A  in  a  nanosecond 


t 
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foe  either  tungsten  or  carbon. 

Using  aq.  XI- 3-19,  we  have  that  the  approximate  total  thickness  of  the 
aiultilayer  stack  is 


D. 

M> 


N  d 


(A) 


Zd  • 

(A) 


(S-2-10) 


Thus,  the  heated  region  will  not  be  a  great  deal  thicker  than  the 
aiultilayer  stack,  so  that  we  can  use  the  thermal  properties  of  the 
aiultilayer  aiaterials  themselves  in  order  to  estimate  the  heat  rise. 

The  smiting  points  of  bulk  tungsten  and  carbon  are  each  about  3500  *K. 
Given  a  heated  voliane  3000A  deep  with  a  cross-sectional  area  200*  m  across 
(that  of  the  collimated  laser  beam) ,  we  can  estimate  that  the  maxinan 
energy  allowed  in  the  beam  is  ~  8-10  *  Joules. 

Conturie's  simulations  Indicate  that  inversion  occurs  at  electron 
densities  of  —  10**  cm  *  (Conturie,  1982) j  this  corresponds  to  an  ion 

if  *■> 

density  of  —  10  an  for  Z  ■  10.  Me  also  assume  that  recombination 
could  bring  as  many  as  1%  of  these  ions  into  inversion  across  the  4-3 
transition.  At  the  upper  limit  in  which  coherent  emission  accounts  for  all 
de-excitation  of  the  tpper  level,  the  total  output  from  a  cylindrical 
anplifier  200*m  in  dimeter  and  1  as  long  muld  then  be  ~  2*10*  Joules. 

Thus,  output  at  the  estimated  image  limit  of  8*10  Joules  could  be 
regarded  as  a  quite  respectable  performance  in  an  initial  demonstration 
experiment. 


A 
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Of  oourae,  even  for  refractory  Materials  like  tungsten  and  carbon,  the 
above  Method  of  damage  estimation  May  be  unrealistic.  For  example,  the 
separation  between  photon  absorption  sites  is  only  about  2QA  under  the 
above  conditions,  so  that  a  large  nunber  of  photoelectrons  will  be 
generated  in  the  aultilayer  stack.  (However,  based  on  a  crude  oatparison 
with  the  x-ray  exposure  characteristics  of  photoresist,  one  can  estimate 
that  the  lifetime  of  these  photoelectrons  will  likely  be  of  order 
femtoseconds,  so  that  the  instantaneous  density  of  secondary  electrons  will 
be  substantially  reduced) . 

It  may  be  possible  to  raise  the  damage  threshold  somewhat  by  using 
curved  mirrors  to  expand  the  beam.  Given  the  short  wavelength,  the  cavity 
would  have  to  have  a  borderline-unstable,  near-aoncentric  geometry,  in 
order  to  produoe  a  narrow  diffracting  waist  at  the  center  of  the  anplifier. 
Further,  since  only  a  fairly  small  nunber  of  reflections  are  likely  to  be 
possible  Airing  the  —  1  nsec  gain  duration,  the  intensity  distribution  on 
the  mirrors  trill  likely  resemble  the  initial  ASE  distribution  for  a 
significant  portion  of  the  shot;  therefore  the  effective  illuminated  area 
would  likely  be  little  larger  than  that  illuminated  by  ASE  despite  the 
curvature. 

According  to  Slaymaker  (1978) ,  the  far-f  ield  WHM  angular  spread  of 
the  ASE  will  be  approximately  equal  to  the  aspect  ratio  of  the  amplifier . 
We  might  estimate  that  if  the  mirrors  are  in  the  mid-field  at  —  1  an  fran 
the  vide  of  an  —  1  an  amplifier,  the  ASE  illuminated  area  at  the  mirrors 
Might  be  a  few  times  the  aperture  area.  A  curved  mirror  cavity  might  then 
succeed  in  increasing  the  damage  threshold  to  an  extent. 


6ino«  each  shot  will  only  Irradiate  a  nail  portion  o f  the  mirror 
■urfaee,  it  will  also  bs  poasibla,  if  neoessary,  to  use  the  multilayers  in 
a  quasi-one-shot  mode  in  *fiich  most  of  the  surface  is  masked  off  during 
each  shot. 

Me  might  mention  at  this  point  that  if  future  devioes  operating 
substantially  above  lasing  threshold  should  exceed  the  threshold  for  mirror 
damage,  it  may  be  possible  to  devise  target  configurations  that  distribute 
the  laser  radiation  over  a  larger  area  of  the  mirrors. 

For  example,  if  one  rotated  the  oooling  slit  90*  out  of  the  plane  of 
fig.  III-2-1,  and  spaoed  a  msnber  of  such  slits  beneath  the  cavity  axis  in 
a  grating-like  fashion,  then  with  an  appropriate  focussing  system  to 
generate  a  focal  line  beneath  each  slit,  one  would  form  an  amplifier  that 
was  segmented  along  the  cavity  axis,  but  expanded  in  the  perpendicular 
diractian. 

For  the  above  reasons,  we  feel  that  there  is  reason  to  be  optimistic 
about  the  question  of  damage  in  lulti layer  cavity  mirrors,  but  of  course 
the  wcertainties  are  such  that  the  issue  can  only  he  settled 
experimentally. 
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Section  XII-3  -  Multilayer  Coating*  for  Short  Wavelength  X-ray 
Microscopy 

Part  A)  Introduction 

In  this  section  we  oonsider  the  possibility  of  using  x-ray  multilayers 
to  form  spectrally  selective  coatings  for  a  short-wavelength  x-ray 
microscope.  To  our  knowledge  this  idea  was  first  aoncieved  at  Lawrence 
Livermore  Laboratory,  but  has  not  been  discussed  quantitatively  in  the 
literature. 

The  microscope  configuration  considered  is  that  invented  fay 
Kirkpatrick  and  Baez  (see  fig.  III-3-1,  after  fig.  8  of  Franks,  1977)  and 
is  one  of  the  first  x-ray  optical  systems  ever  developed  (Kirkpatrick  and 
Baez,  1948) .  This  configuration  is  used  at  HE  to  image  the  x-ray  amission 
from  laser  fusion  targets.  We  will  discuss  our  own  investigation  of 
possible  applications  for  multilayer  coatings  in  the  U£  system. 

At  present,  each  of  the  spherical  substrates  in  the  system  (radii  of 
curvature  about  3600  an)  is  coated  with  a  single  nickel  layer,  which  is 
illuminated  at  a  grazing  (suberitical)  angle  to  the  surface  i,  .  At  mall 
£  the  spherical  substrates  are  effectively  cylindrical  because  the 
curvature  transverse  to  the  beam  posse ses  virtually  no  power. 

In  the  longitudinal  meridian,  however,  the  substrate  can  be  shewn  to 
have  a  focal  length  R-  i,  /2  (Kirkpatrick  and  Baez,  1948).  Thus,  two 
spherical  substrates  with  crossed  orientations  can  provide  the  analogue  of 
s  thin  lens.  One  may  alto  regard  the  second  element  as  correcting  for  an 
extreme  asti^natim. 
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In  the  configuration  in  use  at  U£,  %  ie  about  0.8*,  the  working 
distance  is  about  20  cm,  and  the  Magnification  is  about  5. 

Since  the  individual  substrates  are  operating  very  far  off-axis 
(89.2*),  the  geometrical  performance  of  the  overall  system  is  very  poor  by 
conventional  standards.  Por  a  5  Mm  resolution  at  a  20an  working  distance, 
spherical  aberration  limits  the  aperture  to  about  60  Mm  when  £  ■  0.8*, 
and  ocma  limits  the  field  of  view  to  about  200  Mm  (Hopkins,  1981).  This 
field  of  view  is  very  small,  but  is  about  the  size  of  a  laser  fusion 
target. 

The  aberrations  become  less  severe  if  £  is  increased;  however  at  a 
given  wavelength,  a  single  layer  cannot  reflect  radiation  at  angles  larger 
than  critical.  At  an  angle  of  0.8*  to  the  surface,  the  shortest  wavelength 
that  the  present  nickle  ooating  can  reflect  is  about  2.5A. 

In  this  section  we  discuss  the  possibility  of  converting  the  single 
layer  coating  to  a  spectrally  selective  multilayer  coating  in  order  to 
bring  the  operating  wavelength  down  to  to  the  1.5A  regime,  while 
Maintaining  the  angle  of  incidence  at  about  0.8* . 

The  spectral  selectivity  can  be  very  useful  in  laser  fusion  diagnostic 
applications.  Por  exanple,  the  spectral  selectivity  permits  discrimination 
between  the  x-rays  emitted  in  the  oore  of  s  laser  fusion  target,  and  those 
•mitted  by  dements  in  the  target's  shell.  (At  short  wavelengths,  most  of 
the  x-rays  are  emitted  into  isolated  spectral  lines  having 
epeci es-dependen t  wavelengths,  rather  than  into  oontimnm) .  The  spectral 
selectivity  could  also  be  used  to  discriminate  betoken  x-rays  mitted  by 
the  primary  target  and  those  mitted  by  a  backlighting  source  (Richardson, 


l 
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1981} . 

With  Multilayer  coatings  and  abort  wavelength  illumination,  ana  can 
operate  the  microscope  at  angles  of  up  to  a  few  degrees  from  the  aurfaoe 
before  the  required  Id-spacing  beocnes  iapractically  nail. 

Hcwever,  eq.  II-3-23  shows  that  the  acoeptanoe  angle  or  field  of  view 
of  the  coating  will  decrease  as  the  angle  to  the  aurfaoe  is  increased.  We 
nay  regard  a  single  layer  eoating  as  a  Uniting  version  of  a  Multilayer,  so 
eq.  11-3-23  iaplies  that  the  aooeptanoe  angle  of  a  Multilayer  coating  Must 
necessarily  be  nailer  than  that  of  the  single  layer.  As  (  is  increased, 
the  limitation  imposed  on  field  of  view  by  the  coating  will  eventually 
beocane  equal  too  that  inpoeed  by  ana;  at  this  point  one  begins  to  trade  off 
field  of  view  for  increased  collection  aperture  and  resolution  (the 
spherical  aberration  continues  to  ii^rove  as  £  is  increased). 

In  the  present  application,  we  wish  to  nintain  the  field  of  view  at 
the  —  15041*  limit  irposed  by  geometrical  aberrations  in  the  single  layer 
systn;  this  requires  that  $  main  at  0.8* ,  which,  it  turns  out,  is  the 
angle  at  which  the  acoeptanoe  angle  of  a  tungsten-carbon  Multilayer  baser.*? 
approximately  equal  to  the  angular  field  alined  by  eons  (ace  below) . 

In  future  ays  toms,  it  might  be  worthwhile  to  exploit  the  trade-off 
between  coating  acceptance  angle  and  increased  collection  aperture,  by 
having  the  system  operate  at  an  angle  of  2*  or  bo;  such  a  system  could  be 
used  to  obtain  high  resolution,  narrow  field  images  of  the  snail  eempressed 
core  of  a  laser  fusion  target.  Further,  if  the  hard  ir-ray  mission  from,  a 
source  were  broad-band,  it  would  be  possible  to  use  a  larger  (  ay* ter 
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without  decreasing  the  field  of  view. 

The  alignment  and  fabrication  tolerances  in  the  larger  £  systems 
would  be  sore  stringent  than  those  in  the  J  •  O.B*  system  discussed  here, 
however. 

The  max  mm  angle  at  which  a  single  layer  coating  can  reflect  is  given 
by  (Bilderback,  1982) 


1.6  *  10~3 


which  at  ^  *20  gnv'an3  and  X  ■  1.66A  gives  (  »  0.7*.  Since  the  present 
system  is  designed  to  operate  at  an  angle  of  0.8*,  the  principle  advantage 
gained  with  multilayer  coatings  over  high  density  coatings  in  our 
application  is  that  multilayers  are  spectrally  selective. 
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Part  B)  Calculation  of  Performance 

Fig.  XII~3-2  shows  the  calculated  spectral  and  angular  selectivity  of 
tungsten/carbon  sultilayer  coatings  designed  to  convert  the  112 
Kirkpatrick-Baez  microscope  f ran  operation  at  a  minium  wavelength  of  2.5A 
to  operation  at  1.66 A  (Ni-K  H  ) .  The  calculations  use  X  *1.66X1  which  is 
the  wavelength  of  1 cm  temperature  Ni-K  M  amission.  The  actual  line  used 
would  be  the  Be-like  resonance  line,  which  would  have  a  slightly  different 
wavelength. 

The  thicknesses  chosen  for  the  coating  are  d  ■  d  »  4 0.0 A  in  the 

w  c 

front  mirror  (2d  *  160A),  and  dw  ■  dc  ■  41. 5A  for  the  rear  mirror 
(2d  *  166A) . 

It  might  ultimately  be  worthwhile  to  reoptimize  the  thickness  ratio 
slightly;  however,  large  changes  in  the  ratio  would  cause  a  decrease  in 
performance.  A  multilayer  with  a  substantially  thinner  W  layer  would  have 
higher  reflectivity,  but  only  at  the  expense  of  quite  eeverely  restricting 
the  field  of  view.  A  substantial  increase  in  the  W  thickness  could 
increase  the  field  of  view,  but  only  at  the  expense  of  rapidly  degrading 
the  double  mirror  reflectivity.  A  small  change  in  the  thickness  ratio  front 
d  ^  «d  c  would  cause  cnly  small  changes  in  the  performance  of  the  coating. 

Pig.  III-3-2  alao  assumes  that  different  coatings  are  used  on  the  two 
mirrors  of  the  Kirkpatrick-Baez  pair;  this  oorpensates  for  the  difference 
in  wgle  of  incidence  that  results  from  the  different  distances  of  the  two 
mirrors  from  the  target. 


Spectral  and  Angular  Sensitivity 
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Only  about  4  layer  pair*  are  required  to  approach  the 
absorption-limited  peak  reflectivity.  Comparison  of  this  value  with 
eq.  II-3-19  ahcvs  that  the  aoft  x-ray  approximations  of  sec.  IX-3  are 
becaning  inaccurate  at  the  present  short  wavelength. 

In  fact,  the  number  of  layers  participating  in  the  reflection  is 
becaning  sufficiently  small  that  inaccuracy  in  the  basic  x-ray  formalist! 
derived  from  eq.  XI-1-20  is  becaning  noticeable. 

We  have  found  that  our  formalisn  is  still  accurate  to  about  It  in  its 
calculation  of  peak  reflectivity,  but  beaames  increasingly  inaccurate  in 
the  wings;  in  addition  the  calculated  position  of  the  reflectance  peak  is 
off  by  about  15%.  For  these  reasons  we  have  performed  the  final 
calculations  in  this  section  numerically,  with  the  x-ray  formulas  used  only 
to  determine  trial  values  for  the  parameters.  The  numerical  calculations 
are  carried  out  with  the  non  x-ray  characteristic  matrix  solution  of 
eq.  II-1-22. 

The  calculations  use  densities  measured  by  Barbee  (1981)  of  15  gm/ar' 
and  2  gm/anJ  for  the  tungsten  and  carbon  layers,  respectively.  The  real 
index  of  refraction  used  for  tungsten  is  based  an  anomalous  dispersion  data 
obtained  by  Brentano  and  Baxter,  ms  quoted  in  Table  XV.  11  off  Janes  (1965) . 
The  real  index  for  carbon  was  calculated  with  the  free  electron  model.  The 
imaginary  indicies  for  tungsten  and  carbon  are  based  on  data  obtained  by 
Henke  et  al.  (1982). 

The  optical  constants  used  are 


■J  -* 

A  *  -4.JZ  *  JO  *  i  4.14  »  10 
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Ac  =  ~  7.44  »  JO  «  • 
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Despite  the  small  lumber  of  layer  pairs  required  st  the  center  of 
field,  we  have  chosen  to  use  35  pairs  in  the  system.  Even  with  35  pairs, 
there  is  still  a  faint  residual  ringing  in  the  ijoage  outside  the  ±  150 *.m 
field.  (This  is  a  consequence  of  extinction  being  due  to  depletion  by 
back-reflection,  rather  than  by  absorption) . 

The  worst  drawback  to  the  coating  is  probably  the  limitation  that  its 
finite  acceptance  angle  inposes  on  the  microscope's  field  of  view.  This 
limitation  is  probably  tolerable;  the  angular  IHHM  of  double  reflectance 
for  an  ideal  coating  is  0.08  degrees  at  X  -1.66A,  which,  as  indicated, 
leads  to  a  1 150  micron  field  of  view  at  a  working  distance  of  21as;  this 
represents  only  a  slight  decrease  from  the  field  of  view  allowed  by 
aberrations.  However,  the  narra#  field  leaves  little  margin  for  error  in 
designing  and  fabricating  the  coating. 

The  angular  boundaries  of  the  sene  of  high  double  reflectance  from  the 
multilayer  are  reasonably  sharp,  even  though  the  rone  is  somewhat 
asymmetric.  (Because  of  the  asynmetry,  the  angle  of  peak  reflectivity  has 
been  decentered  f ran  the  optical  axis) .  The  wavelength  response  of  the 
coating  is  fairly  broad.  Neglecting  dispersion,  the  nominal  wavelength 
range  (TWtt  of  the  double  reflectance  curve)  is  from  1.5A  to  l.BA.  (Kith 
dispersion  included  the  range  becomes  1.55A  to  1.75A.) 

The  tolerance  on  the  layer  thicknesses,  as  far  as  maintaining  high 
reflectivity  for  an  OT-axis  point  is  concerned,  is  cooperatively  mild, 
about  1  4A.  However,  despite  the  relatively  wide  range  of  layer 
thicknesses  that  can  provide  high  on-axis  reflectivity,  there  is  little 
margin  for  error  in  the  alignment  of  the  zone  of  high  reflectance  with  the 


optical  axis  (assisting  a  quasi-aonochronatic  aouroa) .  Therefore  it  would 
be  preferable  to  be  able  to  accurately  calibrate  the  deposition  system  in 
an  effort  to  obtain  layer  thicknesses  accurate  to  about  1A. 

The  9nm  aeparation  between  the  two  mirrors  implies  a  difference  in 
incidence  angle  of  0.03  degrees.  Fortunately  this  is  somewhat  smaller  than 
the  0.08  degree  FVHM  of  the  coating.  However (  because  of  the  need  to 
center  the  high  reflectance  cones  of  both  coatings  accurately,  one  would 
prefer  that  the  two  coatings  be  made  with  slightly  different  2d-spaeings. 

To  consider  the  effect  of  fabrication  error,  suppose  that  the  layer 
thicknesses  are  39A  instead  of  4QA,  and  that  the  sane  coating  is  used  on 
both  mirrors.  One  then  finds  that  the  optical  axis  intersects  the 
half-maximum  of  the  double  reflectance  curve,  so  that  roughly  speaking  only 
half  of  the  field  will  be  imaged  (i.e.  the  angular  half-width  of  the  field 
is  aero) .  If  the  second  coating  is  independently  optimized,  and  has  the 
desired  layer  thicknesses  of  41. SA,  the  angular  half-width  in  the 
previously  is>- imaged  direction  becomes  0.025  degrees  (given  the  ue  1A 
systematic  error  in  the  layers  of  the  front  mirror) .  If  the  intended 
thicknesses  of  the  layers  of  the  second  coating  are  41. 5A,  but  they  instead 
have  an  error  equal  to  that  in  the  first  coating  (so  that  their  thickness 
is  40. 5A),  the  angular  width  becomes  0.021  degrees.  If  the  layers  of  the 
second  coating  have  the  same  errors  in  the  apposite  direction  (i.e.  a 
thickness  of  42. SA),  the  half-width  becomes  0.03  degrees.  (As  amntioned 
above,  with  no  errors  in  the  independently  specified  coatings,  the 
half-width  is  0.04  degrees.) 
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Thus,  the  use  of  •  different  target  thickness  in  the  rear  coating  will 
definitely  reduoe  the  sensitivity  to  error. 

The  coatings  should  be  relatively  insensitive  to  other  kinds  of 
defects  in  structure,  due  to  the  large  2d-spaeings. 

Prom  the  36m  radius  of  curvature,  the  60  micron  aperture,  and  the 
working  distance  of  2 lan,  one  can  calculate  that  that  the  variation  in 
angle  of  incidence  across  the  aperture  due  to  curvature  of  the  mirrors  is 
0.007S  degrees.  (The  main  variation  is  due  to  displacement  in  the 
longitudinal  direction) . 

Thus,  the  narrow  0.08*  bandwidth  of  the  coatings  will  not  reduce  the 
aperture  of  the  system. 


